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Abstract 


This  document  reports  the  results  associated  with  a 
magnetostrict ive  fiber  optic  interf erometric  sensor  which  was 
fabricated  and  evaluated.  To  evaluate  the  sensor's  performance, 
a  Helmholtz  coil  was  designed  and  constructed  to  generate  a 
uniform  magnetic  field.  The  sensor's  response  was  determined  by 
examining  the  amplitude  of  the  interferometer's  output  light. 

DC  magnetic  fields  could  be  detected  due  to  excessive  noise 
contributed  by  the  fused  coupler  component  in  the  interferometer 
AC  measurements  were  successfully  accomplished  for  a  25  Hz 
magnetic  field  by  inspecting  the  modulation  of  the 
interferometer's  output  light  at  25  Hz.  Sensors  composed  of  RF 
sputtered  metal  coated  fibers  did  not  detect  the  AC  magnetic 
field  because  the  thin  metal  films  did  not  generate  enough  force 
to  change  the  lengch  of  the  optical  fiber  a  detectable  amount. 
The  sensor  did  detect  the  AC  magnetic  field  when  the  sensor  arm 
was  composed  of  a  bulk  nickel  rod.  Thermal  annealing  the  nickel 
rod  in  the  presence  of  an  external  magnetic  field  improved  the 
sensor's  response.  Thinning  the  optical  fiber  in  the  sensor  arm 
did  not  significantly  affect  the  sensor's  performance,  but 
increasing  the  length  of  the  optical  fiber  did  enhance  the 
sersor's  response.  The  sensor's  performance  also  improved  when 
the  AC  magnetic  field  was  superimposed  on  a  DC  bias  magnetic 
field.  In  the  sensor1''  optimum  configuration,  an  AC  magnetic 
field  flux  density  of  10  mG  was  detected. 
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Investigation  of  a  Magnetostrictive 
Fiber-Optic  Interferometric  Sensor 

I.  Introduction 


Background 

Magnetic  fields  generated  by  human  mental  activity  can 
provide  useful  physiological  information,  such  as  an  indication 
of  pilot  workload  or  loss  of  consciousness  (1:1).  However,  the 
flux  densities  of  these  weak  biomagnetic  fields  only  span  one 
decade  (10  9  to  10  8  Gauss)  (1:8)  .  (Appendix  A  provides  an 
explanation  of  magnetic  flux  density  and  magnetic  field  units) . 
The  intensity  of  these  fields  are  extremely  small  compared  to  the 
ambient  magnetic  field  which  has  a  flux  density  of  approximately 
10  1  Gauss  (2:360).  The  current  state-of-the-art  method  of 
detecting  small  magnetic  fields  involves  using  a  superconducting 
quantum  interference  device  (SQUID)  (1:8).  This  instrument, 
however,  is  expensive,  bulky,  and  must  be  cryogenically-cooled . 

An  overview  of  the  other  available  magnetic  field  sensing 
approaches  and  their  limitations  is  thoroughly  discussed  in 
Chapter  II  (Literature  Review).  An  alternative  technology,  known 
as  a  magnetostr ictive  fiber  optic  interferometric  sensor,  has 
been  proposed.  Theoretically,  this  technology  should  be 
sufficiently  sensitive  to  detect  the  weak  biomagnetic  fields 
without  the  disadvantages  associated  with  the  competing  existing 
technologies  (3:87). 
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The  magnetostrictive  fiber  optic  sensor  concept  is  based 
upon  the  classic  Mach-Zehnder  interferometer  (3:87).  To 
implement  this  technology,  one  fiber  is  coated  with  a 
magnetostrictive  material,  and  it  serves  as  the  sensor  arm.  A 
second  fiber  is  uncoated,  ana  it  is  used  as  the  reference  arm. 

In  practice,  a  magnetic  field  will  deform  the  magnetostrictive 
coating  covering  the  fiber  in  the  sensor  arm,  and  this  induced 
deformation  will  cause  perturbations  in  the  light  propagating  in 
that  arm.  These  perturbations  can  subsequently  be  detected  via 
the  constructive/destructive  interference  pattern  that  is 
generated  when  the  sensor  arm's  light  is  recombined  with  the 
light  propagating  in  the  reference  arm.  The  use  of  the 
interference  property  is  responsible  for  the  high  sensitivity  of 
the  sensor  since  changes  in  the  length  of  the  sensor  arm  (on  the 
order  of  the  wavelength  of  the  propagating  light)  can  be 
detected.  The  basic  arrangement  is  illustrated  in  Figure  1.1. 
More  detailed  illustrations  of  the  exact  interferometers  built 
are  included  in  Chapter  IV  (Experimental  Procedure)  . 

The  initial  theoretical  study  postulated  that  the  sensor 
would  be  sufficiently  sensitive  to  detect  a  magnetic  flux  density 
as  small  as  10  12  Gauss,  which  would  certainly  enable  the 
detection  of  biomagnetic  fields  (3:87).  However,  the  prototype 
sensors  fabricated  thus  far  have  not  achieved  the  theoretically 
predicted  minimum  nor  the  level  associated  with  human  biomagnetic 
activity  (4:47).  A  minimum  detectable  magnetic  flux  density  of 
1 . 4  x  10  4  G  at  2 5  Hz  has  been  achieved  by  the  American  Research 
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Signal  Processing 
Electronics 


Figure  1.1.  Interferometer  Arrangement 


1-3 


corporation  of  Virginia,  and  this  value  represents  the  progress 
made  thus  far  (5:36).  Chapter  II  (Literature  Review)  provides  a 
more  comprehensive  discussion  concerning  the  previous 
developments  of  the  magnetostrictive  fiber  optic  sensor. 

Problem  Statement 

This  thesis  project  attempted  to  realize  the  desired  minimum 
detectable  biomagnetic  field  flux  density  of  10’9  Gauss  by 
designing  and  fabricating  an  interferometer  with  optimized  sensor 
arms.  The  frequency  range  of  interest  spans  DC  to  30  Hz. 

The  sensor  arm  in  most  of  the  previously  fabricated 
interferometers  were  of  two  types  (4:48).  The  first  class  of 
sensor  arms  used  a  standard  fiber  optic  cable  epoxied  to  a  slab 
of  bulk  magnetostrictive  material.  This  type  of  sensor  arm  is 
difficult  to  accurately  characterize  and  reproduce  due  to  the 
epoxy  that  serves  as  an  interface  between  the  magnetostrictive 
material  and  optical  fiber.  The  second  type  of  sensor  arms  were 
conventional  fiber  optic  cables  that  were  coated  with  a 
magnetostrictive  material.  Many  methods  have  been  used  to  coat 
the  fiber,  including  electron  beam  evaporation  (6:408), 
electroplating  (7:86),  and  RF  sputtering  (8:32).  However,  no  one 
has  thoroughly  characterized  the  sensor's  performance  with 
respect  to  utilizing  RF  sputtered  materials  and  specific  coating 
geometries.  It  was  the  goal  of  this  research  project  to 
characterize  RF  sputtered  sensor  arms.  For  completeness,  several 
magnetostrictive  coating  materials  and  bulk  magnetostrictive 
materials  were  investigated  along  with  the  influence  of  thermally 
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annealing  the  thin  films  while  they  are  exposed  to  an  externally 
applied  magnetic  field. 

Scope 

The  method  of  depositing  the  magnetostrictive  materials  was 
limited  to  the  highly-rel iable ,  monolithic  silicon  integrated 
circuit  thin  film  RF  sputtering  technology.  The  magnetostrictive 
materials  used  included:  cobalt,  iron,  and  nickel.  These 
materials  all  exhibit  magnetostrictive  properties,  and  they  have 
been  utilized  in  similar  applications  (1:14).  The  wavelength  of 
light  used  to  excite  the  sensor  was  restricted  to  633  nm  to 
reduce  the  number  of  parameters  investigated.  The  influence  of 
wavelength  on  the  sensor's  performance  is  thoroughly  addressed  in 
the  theoretical  analysis  section. 

Assumptions 

The  ultimate  goal  of  this  project  was  to  design,  fabricate 
and  evaluate  the  performance  of  a  magnetostrictive  fiber  optic 
sensor  possessing  sufficient  sensitivity  to  detect  biomagnetic 
fields.  However,  to  test  the  sensor,  the  biomagnetic  fields  were 
simulated.  A  Helmholtz  coil  was  designed,  fabricated,  and 
calibrated  for  this  purpose.  During  calibration,  the  Helmholtz 
coil's  dynamic  magnetic  field  intensity  and  uniformity  were 
determined.  The  magnetic  field  intensity  measured  was  limited  by 
the  fabrication  tolerances  associated  with  the  Helmholtz  coil, 
the  regulation  of  the  applied  current,  and  the  ambient  magnetic 
flux  density,  and  not  by  the  sensitivity  of  the  Hall  effect 
magnetometer  which  was  used  for  calibration  (9:5). 
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sensor  arm 


Approach 

This  project  encompassed  three  critical  tasks: 
design  and  fabrication,  interferometer  instrumentation  design  and 
configuration,  and  the  sensor's  performance  evaluation.  This 
section  only  outlines  the  approach.  The  details  associated  with 
accomplishing  the  various  tasks  are  discussed  in  Chapter  IV 
(Experimental  Procedure) . 

Sensor  Arm  Design  and  Fabrication 

The  first  step  in  designing  and  fabricating  the  sensor  arm 
was  to  prepare  and  examine  the  ends  of  the  fiber  optic  cable. 

The  length  of  the  sensor  arm  was  limited  to  10  cm,  which  was  the 
extent  of  the  uniform  magnetic  field  generated.  Once  the  ends  of 
the  fiber  optic  cable  were  adequate  to  make  optical  connections, 
the  fiber's  outer  glass  cladding  was  cleaned  and  thinned  using 
chemical  etchants.  The  magnetostrictive  materials  were  then  RF 
sputtered.  Coating  variables  included:  thickness,  geometry,  and 
thermal  annealing  conditions,  which  included  annealing  in  the 
presence  of  an  external  magnetic  field.  Additional  sensor  arms 
were  fabricated  by  bonding  a  fiber  to  a  bulk  magnetostrictive 
material  sample. 

Interferometer  Configuration 

This  task  required  connecting  the  fibers  (coated  and 
uncoated)  and  components  to  from  an  interferometer.  All  of  the 
fiber  ends  were  carefully  prepared  to  make  high-quality  optical 
connections.  Once  the  fiber  ends  were  prepared,  the  fiber-to- 


fiber  connections  were  made  using  glass  ferrule  tubes.  All  of 


the  fibers,  couplers,  and  connections  were  evaluated  for  optical 
power  loss  before  the  interferometer  was  assembled,  as  well  as 
during  the  sensor's  fabrication  sequence. 

Sensor  Performance  Evaluation 

Prior  to  evaluating  the  sensor,  the  Helmholtz  coil  was 
designed,  fabricated,  and  calibrated.  The  sensor,  with  each  of 
the  various  sensor  arm  configurations,  was  evaluated  for  its 
response  to  DC  and  AC  magnetic  fields.  These  performance 
measurements  were  accomplished  by  varying  the  applied  magnetic 
field's  magnitude  and  examining  the  interferometer's  output  light 
amplitude  and  frequency  information. 

Plan  of  Development 

Following  this  Chapter,  which  serves  as  an  overview  of  the 
thesis  research,  a  literature  review  is  presented.  Chapter  II 
identifies  the  competing  magnetic  field  sensing  technologies  and 
highlights  the  evolution  and  the  current  status  of  research 
concerning  magnetostrict ive  fiber  optic  sensor  technology.  The 
applicable  theory  is  presented  in  Chapter  III,  which  includes  a 
discussion  of  Helmholtz  coil  design,  Mach-Zehnder  interferometer 
operation,  fiber  optic  cables,  electromagnetic  propagation  in 
cylindrical  waveguides,  magnetostrict ive  materials,  and  the 
sensor  arm  designs.  Chapter  IV  describes  the  experimental 
procedures  used  to  fabricate  and  evaluate  the  interferometric 
sensor  and  the  Helmholtz  coil.  The  results  of  the  experiments 
are  reported  in  Chapter  V.  The  last  Chapter  concludes  with  a 
summary  and  recommendat ions  for  subsequent  research. 
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II.  Literature  Review 

Introduction 

Since  the  initial  concept  for  detecting  magnetic  fields  with 
a  fiber  optic  interferometer  was  reported  in  1980,  dozens  of 
papers  concerning  this  sensor  technology  have  been  published. 
Research  regarding  this  sensor  concept  is  being  conducted  in 
government,  industry,  and  university  laboratories  in  Denmark, 
Great  Britain,  Japan,  and  the  United  States  (10:57).  After 
identifying  the  other  available  methods  of  detecting  magnetic 
fields  and  their  limitations,  this  section  highlights  the 
important  advances  in  the  development  of  the  magnetostrict ive 
fiber  optic  sensor  and  reports  on  the  sensor's  current 
performance  status. 

Competing  Technologies 

Many  methods  are  currently  available  to  detect  magnetic 
fields.  However,  as  this  section  will  illustrate,  none  of  the 
technologies  competing  with  the  magnetostrictive  fiber  optic 
sensor  have  the  same  potential  for  high  sensitivity,  compactness, 
and  low  cost.  The  other  technologies  addressed  include: 
induction  coils,  Hall  effect  magnetometers,  magneto-optic 
devices,  superconductor  quantum  interference  devices  (SQUIDs), 
and  magnetoresistive  devices. 

A  fluxmeter  is  an  example  of  an  induction  coil  that  has  a 
typical  sensitivity  of  10'5  G.  A  fluxmeter  detects  magnetic 
fields  through  the  movement  of  a  coil  placed  in  a  magnetic  field 
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(1:10).  The  coil  rotates  as  it  is  moved  through  an  area  of 
changing  magnetic  flux  density  (1:10). 

A  Hall  effect  magnetometer  detects  magnetic  fields  through 
the  voltage  that  appears  when  a  semiconductor  with  a  current 
passing  through  it  is  placed  in  a  magnetic  field.  The  Hall 
effect  is  often  used  in  laboratory  settings  to  detect  magnetic 
fields  as  small  as  10'5  G  (9:3). 

While  less  sensitive  than  Hall  effect  magnetometers, 

magneto-optic  devices  are  used  to  measure  power  line  current 

(10:60).  Magneto-optic  devices  operate  using  the  Faraday  effect, 

which  states  that  linearly  polarized  light  is  rotated  when  the 

light  propagates  in  special  media  that  are  exposed  to  magnetic 

fields.  Quartz,  special  glass,  ferrite  and  some  semiconductors 

have  been  used  as  the  Faraday  element.  The  magnetic  field 

sensitivity  of  the  magneto-optic  device  depends  on  the  Verdet 

constant  and  propagation  length  of  the  Faraday  element.  A 
.  .  .  -2  • 

sensitivity  of  10  G  is  representative  of  magneto-optic  devices 
(l: 11)  . 

SQUID  magnetometers  currently  are  the  most  sensitive 
magnetic  field  detectors  available,  and  they  possess  a  minimum 
detectable  magnetic  flux  density  of  10’9  Gm^Hz 'k.  There  are  two 
types  of  SQUID  magnetometers:  DC-SQUIDs  and  RF-SQUIDs.  DC-SQUIDs 
consist  of  a  DC  bias  current  flowing  through  a  loop  composed  of 
two  Josephson  junctions.  A  change  in  the  applied  magnetic  flux 
density  can  be  detected  from  the  change  in  current  flowing 
through  the  loop.  RF-SQUIDs  contain  only  one  Josephson  junction 
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iii  a  loop  that  is  excited  with  an  RF  bias  current.  While  SQUIDs 
are  sufficiently  sensitive  to  detect  biomagnetic  fields,  they 
require  cooling  to  4.2°  K,  which  makes  the  device  bulky  and 
expensive  to  operate  because  of  the  liquid  helium  dewar 
requirement.  Since  SQUIDs  cost  hundreds  of  thousands  of  dollars, 
they  are  too  expensive  for  many  research  applications  (1:12). 

In  October  1989,  IBM  and  the  Helsinki  University  reported  that 
they  had  fabricated  a  sensor  capable  of  detecting  magnetic  fields 
generated  by  the  brain  that  was  composed  of  four  SQUIDs  on  a 
silicon  chip  (11:5). 

The  final  technology  competitive  with  the  magnetostrict ive 
fiber  optic  sensor  is  the  magnetoresistive  sensor,  which  detects 
magnetic  fields  using  the  magnetoresistive  effect.  The  change  in 
resistance  of  a  current  carrying  magnetic  material  induced  by  an 
externally  applied  magnetic  field  is  referred  to  as  the 
magnetoresistive  effect  (12:2).  Philips  produces  a 
magnetoresistive  sensor  that  uses  a  Wheatstone  bridge  arrangement 
composed  of  magnetic  materials.  Philips  claims  their  detector  is 
one  and  one  half  times  more  sensitive  than  Hall  effect 
magnetometers  (12:6). 

Initial  Proposal 

Yariv  and  Winsor  first  proposed  the  magnetostrict ive  fiber 
optic  sensor  concept  in  March  1980  (3:87-89).  In  their  paper, 
they  described  the  sensor's  configuration,  the  principle  of 
magnetostriction,  and  calculated  the  minimum  magnetic  field  flux 
density  that  could  be  detected  using  their  theoretically  modeled 
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sensor.  Their  sensor  was  proposed  using  nickel  as  the 
magnetost.r  ict  ive  material.  By  assuming  that  the  magnetic  field 
to  be  detected  was  a  small  AC  field  superimposed  upon  a  larger  DC 
field,  they  stated  that  the  change  in  the  phase  of  light 
propagating  in  the  sensor  arm  would  be  directly  proportional  to 
the  length  of  the  sensor  arm  and  the  magnitude  of  the  AC  field. 

To  determine  the  minimum  detectable  magnetic  field,  they  assumed 
that  the  dominant  noise  source  would  be  shot-noise  in  the  light 
detector.  Shot-noise  is  the  result  of  random  electron  emissions 
in  the  electrical  components  (13:104).  Using  an  excitation  laser 
source  with  an  output  power  level  of  1  mW,  an  emission  wavelength 
of  1  micron,  a  detector  possessing  a  quantum  efficiency  of  0.5, 
and  a  detection  bandwidth  of  1  Hz,  they  calculated  a  minimum 
detectable  magnetic  field  flux  density  of  1.6  x  10 9  Gauss  for 
each  1  meter  length  of  the  coated  sensor  arm. 

Experimental  Verification 

One  month  after  Yariv  and  Winsor  published  their  proposal 
for  the  sensor  concept,  Dandridge  and  others  at  the  Naval 
Research  Laboratory  published  the  first  experimental  evidence 
that  such  a  sensor  would  work  (6:408-409).  They  examined  the 
response  of  two  classes  of  sensor  arms:  fiber  optic  cables 
epoxied  around  a  nickel  tube,  and  1.5  micron  thick  electron  beam 
evaporated  films  of  nickel  that  were  deposited  directly  on  the 
surface  of  the  fiber  and  then  annealed  at  900°C  in  hydrogen. 

Their  measurements  indicated  that  the  sensor's  sensitivity 
increased  with  the  frequency  of  the  AC  magnetic  field.  They 
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attributed  this  increase  in  sensitivity  to  the  interferometric 
detection  system  and  not  to  the  magnetostrictive  material.  Fiber 
optic  interferometric  systems  are  less  sensitive  at  lower 
frequencies  due  to  the  1/f  noise  in  the  light  detector  (14:141). 
They  also  found  that  the  sensors  using  bulk  magnetostrictive 
material  were  three  orders  of  magnitude  more  sensitive  than 
sensors  using  thin  film  magnetostrictive  material.  Their  bulk 
magnetostrictive  material  sensors  exhibited  a  linear  response 
with  magnetic  field  variations,  and  peak  sensitivities  of 
8  x  10  G  for  each  one  meter  length  of  the  coated  fiber  at 
10  kHz,  and  10'5  G  at  100  Hz. 

Elimination  of  Interferometer  Drift 

The  first  improvement  to  the  magnetostrictive  fiber  optic 
interferometric  sensor  was  adopted  from  a  technique  employed  for 
an  optical  fiber  hydrophone.  In  September  of  1980,  Jackson  and 
others  at  the  Naval  Research  Laboratory,  published  a  report 
concerning  the  utilization  of  a  piezoelectric  cylinder  to 
maintain  the  operation  of  a  fiber  optic  interferometer  at  the 
quadrature  point  (15:2926-2929).  An  interferometer  is  most 
sensitive  to  changes  in  the  relative  phase  of  the  propagating 
light  intensity  when  operated  at  the  quadrature  point;  this 
concept  is  illustrated  in  Figure  II. 1. 

In  their  paper,  they  described  attaching  the  reference  fiber 
and  the  sensor  fiber  around  piezoelectric  cylinders  whose 
diameter  changes  with  an  applied  voltage.  By  changing  the 
applied  voltage  to  the  piezoelectric  cylinder,  the  relative  phase 
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Figure  II. 1.  Quadrature  Operation 


of  the  light  in  one  arm  could  be  changed  until  the  quadrature 
point  was  obtained.  For  optimum  operation  of  the  sensor,  the 
ability  to  change  the  relative  phase  of  the  output  signal  is 
necessary  for  two  reasons.  First,  the  lengths  of  the  sensor  and 
reference  arms  cannot  be  precisely  cut  to  within  one  wavelength 
of  the  excitation  sourc*  .  Second,  when  the  sensor  arm  is  coated 
with  a  magnetost rict ive  material,  the  sensor  arm  will  change 
dimensions,  depending  upon  the  thermal  coefficient  of  expansion 
of  the  material  and  the  ambient  temperature  conditions.  In  the 
interferometer  they  fabricated,  the  relative  phase  of  the  light 
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propagating  in  the  two  arms  could  be  varied  from  10  s  to  1 
radian.  The  phase  change  was  linear  at  a  rate  of  2 it  per  volt 
applied  to  the  piezoelectric  cylinder. 

Materials  Investigation 

Once  Dandriage  and  others  had  proven  that  the  sensor  idea 
would  work  in  the  laboratory,  many  scientists  began  investigating 
improvements  to  the  sensor's  performance.  The  use  of  metallic 
glasses  as  the  magnetost rict ive  material  in  the  sensor  arm  has 
been  identified  as  one  of  the  major  achievements  in  the 
development  of  the  sensor  (10:57).  Koo  and  Sigel  of  the  Naval 
Research  Laboratory  first  reported  the  use  of  metallic  glasses  as 
the  magnetostrict ive  material  in  July  1982  (16:334-336).  In 
their  paper,  they  describe  sensor  arms  fabricated  by  bonding 
fiber  optic  cables  to  strips  of  commercially  available  metallic 
glasses,  such  as  2605  SC  produced  by  the  Allied  Chemical  Company. 
They  observed  that  annealing  the  metallic  glass  sensor  arms 
affected  the  performance  of  the  sensor.  The  reported  minimum 
detectable  magnetic  flux  density  was  5  x  10~9  G  for  each  one 
meter  length  of  coated  optical  fiber  (with  an  applied  1  kHz 
magnetic  field). 

Recent  Patent 

Although  fiber  optic  magnetic  field  sensors  are  close  to 
achieving  the  theoretical  minimum  detectable  magnetic  field,  many 
laboratories  are  continuing  to  investigate  methods  to  improve  the 
sensor's  performance.  The  SpecTran  Corporation  received  a  patent 
in  1987  for  a  unique  method  of  bonding  the  magnetostrict ive 
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material  to  the  fiber  optic  cable  in  the  sensor  arm  (17:97). 
SpecTran's  method  of  material-fiber  bonding  requires  a  new  type 
of  fiber  optic  cable  that  is  illustrated  in  Figure  II. 2.  The 
conventional  fiber  optic  cable  used  in  the  sensor  consists  of  a 
concentric  arrangement  of  a  light  carrying  glass  core,  glass 
outer  cladding,  and  the  magnetostr ict ive  jacket.  On  the  other 
hand,  SpecTran's  sensor  arm  fiber  optic  cable  consists  of  a 
central  core  of  magnetostrict ive  material,  an  inner  glass 
cladding,  a  middle  light  carrying  glass  cladding,  and  finally,  an 
outer  glass  cladding.  By  placing  the  magnetostrict ive  material 
inside  the  fiber  optic  cable  instead  of  outside  the  cable, 
SpecTran  predicts  that  the  fiber  optic  cable  will  be  more 
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sensitive  to  length  changes  induced  by  an  applied  magnetic  field. 
SpecTran  has  not  published  any  performance  results  to  date. 
Current  Work 

In  addition  to  new  fiber-material  bonding  techniques,  many 
laboratories  are  continuing  to  research  new  materials  and  new 
sensor  configurations.  An  example  of  the  current  research 
directed  at  detecting  biomagnetic  fields  is  the  work  being 
conducted  at  the  American  Research  Corporation  of  Virginia. 
Previously,  they  fabricated  a  sensor  using  fiber  optic  cables 
epoxied  to  bulk  metallic  glass  and  achieved  a  minimum  detectable 
magnetic  flux  density  of  1.4  x  lO'4  G  at  25  Hz  for  a  sensor  arm 
having  a  one  meter  coated  length  (5:36).  '"’heir  current  research 
involves  the  investigation  of  magnetostrictive  materials  and 
signal  processing  techniques  (18:2).  They  are  investigating 
vapor  deposition  and  annealing  of  metallic  glasses  and 
fabricating  new  sensor  arms  using  Terfenol ,  a  rare  earth  magnetic 
alloy  (Tb0  3Dy1  7Fe2  9)  ,  as  the  magnetostrictive  material.  The  use 
of  lock-in  amplifiers  is  being  studied  as  one  method  to  improve 
sensitivity.  The  lock-in  amplifier  is  used  to  amplify  only  those 
signals  from  the  light  detector  at  the  frequency  of  the  applied 
magnetic  field. 
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III.  Theoretical  Analysis 

This  chapter  contains  the  theory  that  is  pertinent  to  the 
sensor's  performance  and  evaluation.  The  chapter  begins  with  a 
discussion  of  the  concepts  used  to  design  the  Helmholtz  coil. 

Once  the  method  of  generating  the  magnetic  field  is  explained, 
the  magnetostrictive  materials,  which  sense  the  magnetic  field, 
are  described.  The  major  component  in  the  interferometer,  the 
fiber  optic  cable,  is  then  described  as  an  introduction  to  a 
discussion  of  light  propagation  in  a  cylindrical  waveguide. 

Next,  the  theory  concerning  an  interferometer  is  presented.  The 
results  of  the  chapter  are  summarized  in  a  section  listing  the 
critical  performance  parameters  of  the  sensor,  which  is  followed 
by  the  intended  sensor  arm  designs  for  this  research  project. 
Helmholtz  Coil  Design  (19:152-158) 

In  order  to  evaluate  the  performance  of  the  magnetostrictive 
fiber  optic  sensor,  the  biomagnetic  flux  density  will  be 
simulated  using  a  Helmholtz  coil,  which  consists  of  two  coils  of 
equal  size  arranged  perpendicular  to  a  common  axis  as  illustrated 
in  Figure  III.l.  The  Helmholtz  coil  is  often  used  in  scientific 
research  to  produce  a  relatively  uniform  magnetic  field  whose 
magnitude  is  controlled  by  the  amount  of  current  flowing  through 
the  two  coils,  the  number  of  turns  in  each  coil,  the  radius  of 
each  coil,  and  the  distance  separating  them. 

The  magnetic  field  generated  by  a  Helmholtz  coil  can  be 
explained  using  Figure  III. 2  and  the  Biot-Savart  law  which 
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N  -turns 


N- turns 


Figure  ITI.l.  Helmholtz  Coil  (19:156) 


states  (bold  indicates  vector  quantity) : 

b,.2>  -  y-  L  - 


r2  '  r. 


where , 


B(r2)  =  magnetic  flux  density  vector  at  a  point  of 


interest  on  coil  2 


H  =  permeability  of  free  space 


I1  =  magnitude  of  t*e  current  in  coil  1 

r1  =  vector  measured  from  an  arbitrary  reference  point  to  a 
point  on  coil  1 

r2  =  vector  measured  from  an  arbitrary  reference  point  to  a 
point  of  interest  on  coil  2 


dlrX  (rrr,) 


Figure  III. 2.  Biot-Savart  Law  Geometry  (19:153) 


dl1  =  vector  line  integral  increment. 

In  order  to  describe  the  magnetic  flux  density  that  can  be 
generated  in  a  multiple-turn  Helmholtz  coil,  the  Biot-Savart  law 
is  first  applied  to  determine  the  on-axis  magnetic  flux  density 
of  a  single-turn  of  wire  with  the  geometry  shown  in  Figure  III. 3 
Since  the  magnetic  flux  density  for  off-axis  points  cannot  be 
calculated  in  closed  form,  numerical  techniques  are  used,  and 
they  will  be  explained  following  the  on-axis  calculation.  Based 
on  the  geometry  of  the  problem  (Figure  III. 3),  the  following 
substitutions  into  equation  1  can  be  made: 
dl,  =  ad# (-isintf  +  jcosp) 
r2  -  r1  =  -iacosi?  -jasin«  +  kz 
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Figure  III. 3.  Axial  Magnetic  Flux  Density  for  a  Single-Turn 

of  Wire  (19:155) 

I  r2  -  r1 1  =  (a2  +  z2)\ 

Thus,  with  the  substitutions,  the  Biot-Savart  law  for  the  singl 
turn  of  wire  becomes: 


rT  (izacosf  +  izasing  +  ka1-' )  d* 
=  -  I  -2  ,  3/2 


47 t  J  ,,  (z  +  a") 


(2) 


Using  the  fact  that: 


J7  sin^dfl  -  J2t  cos?d?  =  0, 


equation  2  simplifies  to: 


J±rM- 


B(Z)  =  ■  7  2  3/? 

2 ( z "  +  a2)5" 


(3) 


Consequently,  this  equation  can  be  used  to  predict  the  magnetic 
field  produced  by  two  circular  on-axis  coils  with  the  same 


radius.  Using  the  geometry  shown  in  Figure  III.l,  the  magnetic 
flux  density  at  a  point  on-axis  between  the  two  N-turn  coils 
contributed  from  the  N-turn  lower  coil  is  given  by: 

Nu.Ia2 

B(2)  =  2  ( z2  +  a2) 3/2  *  (4)  • 

To  determine  the  magnetic  flux  density  contribution  from  the 
upper  N-turn  coil,  an  alternate  reference  system  is  used.  The 
new  z-coordinate  origin,  z',  is  located  2b  above  z  as  shown  in 
Figure  III.l.  The  magnetic  flux  density  at  a  point  on-axis 
between  the  two  N-turn  coils  that  is  contributed  from  the  upper 
N-turn  coil  is  given  by: 


B(z)  = 


_ MJJl _ 

2  (z  '2  +  a2)  3/2 


Since  z'  =  2b-z,  equation  (5)  can  be  rewritten  as: 

_ NiiJLaf _ 

B(Z)  "  2(<2b-z)2  +  aV'2  *  (6>’ 

Combining  the  two  expressions  for  the  magnetic  flux  density 
contributions  from  the  upper  and  lower  N-turn  coils  (equations  4 
and  6)  yields  the  magnitude  of  the  total  magnetic  flux  density  at 
a  point  on-axis  which  is  given  by: 


B  (  z )  = 


Nu  la* 


(z  +  a  ) 


z  sir?  + 


[  ( 2b-z )  2  +  a2] 3/2 


The  purpose  of  the  Helmholtz  coil  is  to  produce  a  uniform 
magnetic  field.  By  determining  the  necessary  conditions  for  the 
first  and  second  derivatives  (with  respect  to  z)  of  equation  (7) 
to  vanish,  the  values  of  coil  separation  (2b)  and  coil  radius  (a) 
can  be  optimized  for  the  most  uniform  region  of  the  magnetic 
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field.  The  first  derivative  of  the  above  expression  with  respect 
to  z  is  given  by: 


Nu  la 

dB  = 


(z2  +  a2) 5/2 


3  C  z-2b ) 

[  ( 2b-z )  2  +  a2]5/2 


This  expression  is  zero  when  z  =  b.  Differentiating  equation  8 
with  respect  to  z  yields: 


3Nu  la 


(z2  +  a2) 5/2  ( z2+a2)  7/2 


[  ( 2b-z)  2  +  a2)  ]5/2 


5  fz-2br 
[  ( 2b-z) 2  +  a2] 7/2 


When  z=b,  equation  9  simplifies  to: 


2a2  -  8b2 
(b2  +  a2)772 


(10)  . 


It  is  noted  that  the  second  derivative  evaluated  with  z=b 
vanishes  when  a=2b.  These  results  indicate  that  the  coil  radius 
should  be  equal  to  the  coil  separation.  Under  these  conditions, 
the  magnetic  flux  density  at  the  midpoint  between  the  coils  is 
given  by: 


B(z=Ha)  = 


8g.NI 


A  Taylor  series  approximation  for  the  magnetic  flux  density 
indicates  how  slowly  it  changes  near  the  midpoint.  The  Taylor 
series  approximation  using  the  first  five  terms  is  given  by: 


B  ( z )  =  B  ( z  =  -j  a )  •<  1  - 


144  (  z  -  Hal  4 


(12)  , 


where  it  is  noted  that  the  middle  three  terms  of  the  expansion 


vanish.  The  above  expression  shows  that  for  a  position  equal  to 
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a  distance  (a/10)  measured  relative  to  the  midpoint  along  the 
center-line,  the  magnetic  flux  density  is  0.99988  of  its  midpoint 
value. 

To  illustrate  how  slowly  the  magnetic  flux  density  changes 
in  the  axial  direction,  the  magnetic  flux  density  was  calculated 
for  the  specific  Helmholtz  coil  fabricated  for  this  experiment 
(an  average  coil  radius  of  10.5  cm  with  250  turns  of  wire).  For 
a  current  flow  of  200  mA,  Figure  III. 4  shows  the  axial  variation 
of  the  magnetic  flux  density  between  the  two  coils.  The  axial 
magnetic  flux  density  at  the  midpoint  between  the  coils  was 
calculated  to  be  4.28  Gauss  (G) . 

To  examine  the  change  in  the  magnetic  flux  density  in  the 
off-axis  direction,  numerical  methods  were  used  since  a  closed- 
form  expression  does  not  exist.  The  starting  point  for  the 
calculation  is  the  differential  form  of  the  Bict-Savart  law;  that 
is  : 


which 
coil 
made : 


dB(r?)  = 


Lr2_-.r,j 


47T 


r2  *  ri 


utilizes  the  geometry  of  Figure  III. 2. 
geometry  of  Figure  III. 3,  the  following 


(13)  , 

For  the  single-turn 
substitutions  can  be 


r,(fr)  =  iacosr  jasinr 

r2(<?)  =  jy  + 

dl ( e )  =  -iacosj  +  jasinc 


where  e ,  y,  and  z  are  specified  for  a  single  point  on  the  coil. 
As  e,  y,  and  z  are  varied  for  points  around  the  coil,  the 


magnetic  flux  density  at  the  observation  point  (P)  is  calculated 
by  summing  the  contribut ions  for  all  the  points  on  the  circle's 
circumference.  This  calculation  was  accomplished  using  the 
numerical  values  discussed  above  for  the  Helmholtz  coil  to  be 
used  in  the  sensor's  experimental  performance  evaluation.  The 
results,  produced  using  MATHCAD  software,  are  illustrated  in 
Figure  III. 5.  The  computer  program,  attached  as  Appendix  B, 
provides  vector  expressions  for  the  magnetic  flux  density  for  the 


off-axis  points,  and  it  is  noted  that  the  vector  quantities  are 
no  longer  only  in  the  z-direction. 
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Figure  III. 5.  Off-Axis  Magnetic  Flux  Density  from  a  Helmholtz 

Coil 


Figure  III. 5  depicts  a  plot  of  the  magnitude  of  the  magnetic 
flux  density  vector  for  the  off-axis  points.  The  on-axis  value 
of  the  magnetic  flux  density  that  was  numerically  calculated 
agrees  with  the  value  calculated  in  closed  form,  to  s:>:  decimal 
places.  This  calculation  illustrates  that  any  small  off-axis 
errors  that  might  be  involved  concerning  the  positioning  of  the 
fiber  optic  cable  along  the  center-line  of  the  Helmholtz  coil 
will  not  drastically  affect  the  magnetic  flux  density  that  the 
fiber  optic  cable  will  experience. 
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This  research  effort  required  a  Helmholtz  coil  that  could 
not  only  generate  a  small  magnetic  flux  density  to  evaluate  the 
performance  of  the  sensor,  but  also  a  large  flux  density  to 
permit  its  calibration  using  a  state  of  the  art  laboratory 
magnetometer.  The  expression  for  the  magnetic  flux  density  at 
the  center  point  between  the  coils  (equation  11)  illustrates  the 
design  trade-offs  that  exists  between  the  following  coil 
parameters : 

i)  number  of  turns  (N)  , 

ii)  current  (I) ,  and 

iii)  coil  radius  (a) . 

Equation  11  is  repeated  here: 


8uNI 

B  (z=ija)  = 


(11)  . 


Other  constraints  imposed  on  the  design  included  limiting: 

i)  the  current  to  550  mA, 

ii)  power  dissipation  due  to  the  wire's  resistance  to  2.5  K , 

and 

iii)  the  maximum  number  of  turns  of  wire  to  250. 

Based  on  these  constraints  and  trade-offs,  the  final  design 
consisted  of  a  10.5  cm  radius,  multiply-tapped  coil  that  was 
wound  with  18  gauge  magnet  wire.  The  taps  at  250- ,  100- ,  10-, 

and  2-turns  facilitated  flexibility  in  generating  the  magnetic 
flux  densities  of  interest.  For  calibration  purposes,  the 
Helmholtz  coil  was  configured  with  250-turns,  and  AC  and  DC 
currents  ranging  from  550  to  50  mA  were  applied  to  produce 
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magnetic  field  densities  spanning  11.78  to  1.07  G. 

For  a  practical  Helmholtz  coil  with  many  turns  of  wire,  the 
coil  separation  and  coil  radius  are  not  constant  for  all  the 
turns  of  wire  due  to  the  finite  size  of  the  wire  in  each  coil. 
Consequently,  the  values  of  coil  separation  and  coil  radius  in 
the  above  calculations  are  average  values.  Figures  III. 6  and 
III. 7  show  the  effects  of  the  extreme  values  of  coil  separation 
(8  to  13  cm)  and  radius  (10  to  11  cm)  for  the  axial  magnetic  flux 
density  relative  to  the  center  position  between  the  coils.  It 
was  assumed  that  the  average  values  contributed  a  negligible 
variation . 

Magnetostrict ive  Materials 

Now  that  the  mechanism  for  generating  the  magnetic  fields 
has  been  explained,  the  material  that  reacts  to  the  magnetic 
field  is  described.  Magnetostriction  refers  to  the  small  change 
in  the  dimensions  of  ferromagnetic  materials  that  are  exposed  to 
an  external  magnetic  field  (15:627).  For  this  research  project,, 
the  change  in  the  length  of  the  material  applied  to  the  outer 
surface  of  a  fiber  optic  cable  is  the  important  parameter.  This 
change  in  length  is  typically  on  the  order  of  parts  per  million. 
This  section  begins  by  describing  the  magnetic  nature  of 
ferromagnetic  materials.  The  origin  of  magnetostriction  is  then 
explained,  which  is  followed  by  a  discussion  of  enhancing  the 
magnetostrictive  property  of  materials. 

Ferromagnetic  Materials.  To  understand  ferromagnetic 
materials,  it  is  necessary  to  understand  the  source  of  magnetic 
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fields  in  materials.  When  an  electron  orbits  the  nucleus  of  an 
atom,  it  produces  a  magnetic  field  just  as  the  current  flowing 
through  a  conductor  produces  a  magnetic  field.  The  electron 
itself  also  produces  a  magnetic  field  from  its  inherent  spin.  In 
nonmagnetic  materials,  the  magnetic  field  produced  by  the  orbit 
and  spin  of  one  electron  is  cancelled  by  the  magnetic  field  of  an 
oppositely  oriented  electron  to  produce  a  net  zero  magnetic  field 
in  the  bulk  material.  However,  in  certain  materials  like  cobalt, 
nickel,  and  iron,  the  magnetic  field  contributions  from  each  atom 
contribute  to  a  net  average  external  magnetic  field  (2:210).  A 
quantum  effect  referred  to  as  the  exchange  coupling  between  atoms 
in  a  ferromagnetic  material,  locks  their  magnetic  moments  in  a 
parallel  configuration,  and  this  regions  is  called  a  domain 
(2:214).  The  size  of  the  domains  vary,  but  average  0.1  to  0.01 
cm  across  in  iron.  Domains  are  separated  by  boundaries  called 
Block  walls  which  are  approximately  fifty  nm  thick  in  iron 
(21:66).  Block  walls  are  composed  of  atoms  whose  magnetic  moment 
orientation  transitions  between  those  of  neighboring  domains 
(20:815).  Above  a  cerrain  temperature,  known  as  the  Curie 
temperature,  the  domains  in  ferromagnetic  materials  become 
indistinct,  and  the  material  becomes  paramagnetic.  The  Curie 
temperature  and  melting  points  for  the  materials  of  interest  are 
listed  in  Table  III.l.  When  exposed  to  an  external  magnetic 
field,  the  magnetic  moments  in  both  f erromagnet ic  and 
paramagnetic  materials  tend  to  align  with  the  external  field. 

When  cooled  below  the  Curie  temperature,  the  domains  reform  tc 
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Table  III.l.  Curie  Temperatures  and  Melting  Points  (20:723  and 
22:94) 


Material  Curie  Temperature  f°C)  Melting  Point  (°C) 


Co 

1130 

1492 

Fe 

770 

1535 

Ni 

358 

1453 

the  characteristic  of  ferromagnetic  materials  (2:214). 

Magnetostriction .  A  change  in  length  of  ferromagnetic 
materials  occurs  when  the  domains  rotate  and  move  as  they  align 
with  an  external  magnetic  field  (20:632).  Figure  III. 8 
illustrates  a  simplified  model  of  magnetostriction.  The 
rectangles  represent  the  magnetic  domains,  and  the  arrows 
symbolize  the  net  magnetic  moment.  As  shown,  magnetostriction 
refers  to  expansion  and  contraction.  The  amount  of  change  in 
length  in  a  material  depends  upon  the  magnetic  field.  Figure 
III. 9  depicts  the  fractional  change  in  length  (A  =  £.1/1)  for  the 
materials  of  interest  relative  to  an  externally  applied  magnetic 
field . 

The  direction  of  the  external  magnetic  field  relative  to  the 
magnetic  moment  of  a  single  domain  affects  the  fractional  change 
in  length  and  is  giver,  by  (20:  634): 

A 1  =  3/2As(cos2*  -  1/3)  (14) 

where  A1  =  observed  fractional  change  in  length 

A s  =  saturated  factional  change  in  length  value 
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Figure  III. 8.  Magnetostriction  Model  (23:36  and  20:635) 
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I 


length  and  the  external  magnetic  field. 

Equation  14  shows  that  the  maximum  or  saturated  value  of 
magnetostriction  occurs  when  the  change  in  length  is  measured 
along  the  direction  of  the  external  magnetic  field  ( &  -  0°).  The 
above  equation  also  shows  that  when  the  change  in  length  is 


measured  perpendicular  to  the  magnetic  field  ( e  =  90°),  the 
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transverse  change  in  length  is  one  half  the  negative  of  the 
longitudinal  saturation  value. 

For  an  ensemble  of  domains  whose  initial  magnetic  moments 
are  randomly  arranged,  the  single  cosine  dependence  is  replaced 
by  its  average  value.  For  an  initial  random  distribution  of 
domains  and  no  external  magnetic  field,  <cos20>  =  1/3,  which 
yields  no  magnetostriction.  When  a  strong  external  magnetic 
field  is  applied,  <cos20>  =  0  ,  producing  the  saturated  value  of 
magnetostriction . 

If  the  domains  are  not  initially  randomly  arranged,  the 
following  expression  describes  the  fractional  change  in  length 
(20:634)  : 

A1  =  3/2 As(<cos20 >gv  -  <cos2# >0)  (15) 

where  <cos 20>av  refers  to  the  final  distribution  of  domains,  and 
<cos?j>0  refers  to  the  initial  distribution  of  domains.  This 
expression  reveals  that  the  magnetostriction  is  three-halves  that 
of  the  saturated  value,  and  this  result  can  be  obtained  if  the 
domains  are  all  originally  oriented  perpendicular  to  the  external 
magnetic  field. 

Optimized  Magnetostriction . 

As  implied  in  the  previous  discussion,  there  are  methods 
that  can  be  used  to  maximize  the  amount  of  magnetostriction. 

These  techniques  include  material  selection  and  initial 
orientation  the  of  the  domains  by  stressing  the  material  or 
thermally  annealing  the  material  in  the  presence  of  external 
magnetic  fields. 
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Much  of  the  recent  research  associated  with  the  fiber  optic 
magnetic  field  sensor  has  been  focused  on  developing  new 
magnetostrictive  materials.  The  Allied  Chemical  Company's  2605 SC 
metallic  glass  is  an  example  of  these  new  materials.  As  shown  in 
Figure  III. 9,  the  bulk  material  has  a  greater  value  of 
magnetostriction  at  low  external  magnetic  fields  compared  to  the 
elemental  metals.  This  metallic  glass  is  formed  from  a 
combination  of  boron,  carbon,  iron,  and  silicon  which  is  quenched 
from  a  melt  (8:6).  This  fabrication  method  limits  the  use  of 
developing  metallic  glasses  as  a  magnetostrictive  material  for 
the  fiber  optic  magnetic  field  sensor.  Metallic  glasses  cannot 
be  RF  sputtered  in  thin  layers  of  various  geometries  and  preserve 
their  chemical  composition.  The  magnetostrictive  behavior  of 
metallic  glasses  also  changes  with  time  (8:6).  While  the  single 
elements  have  smaller  magnetostriction  values,  they  were 
investigated  because  they  can  be  deposited  in  thin  layers  with 
various  geometries.  Once  the~~  high  purity  layers  of  single 
elements  are  deposited,  they  can  processed  to  initially  align 
their  magnetic  domains. 

The  first  method  which  can  be  used  to  initially  orient  the 
magnetic  domains  of  ferromagnetic  materials  is  to  apply  a  stress 
to  the  material  (20:636).  Figure  III. 10  shows  the  effect  of 
compression  and  tension  on  the  magnetostriction  of  a  nickel 
sample.  The  data  in  Figure  III. 10  shows  that  the  value  of 
magnetostriction  varies  from  zero  to  three-halves  its  saturation 
value  without  tension.  Under  tension,  the  nickel  domains  were 
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Figure  III. 10.  Magnetostriction  of  Nickel  Under  Stress  (20:636) 


oriented  perpendicular  to  the  direction  of  the  magnetic  field  and 
under  compression,  the  magnetic  field  and  the  domains  were 
aligned  parallel.  The  stress  required  to  align  the  domains  is 
greater  than  the  forces  maintaining  alignment  of  the  domains  in 
the  unstressed  material  (20:634).  These  forces  in  the  unstressed 
material  are  attributed  to  crystal  anisotropy.  Once  an  external 
stress  greater  than  the  crystal  anisotropy  force  is  applied,  the 
domains  align  and  the  magnetostriction  is  given  by: 
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X  -  3/2Ascos 2e  (16) 

where  6  =  angle  between  the  external  magnetic  field  and  stress. 
While  tension  has  proven  to  be  effective  for  domain  control,  it 
was  not  investigated  in  this  research  project  since  the  following 
method  of  domain  control  appears  more  practical  for  a  fiber  optic 
magnetic  field  sensor. 

The  second  method  of  domain  control  is  to  heat  the 
ferromagnetic  material  and  then  allow  it  to  cool  in  the  presence 
of  an  external  magnetic  field  (20:637).  When  the  ferromagnetic 
material  is  heated  past  its  Curie  temperature,  the  domain 
boundaries  disappear.  If  the  material  is  then  allowed  to  cool  in 
the  presence  of  an  external  magnetic  field,  the  domains  should 
reform  and  be  aligned  with  the  external  field.  Figure  III. 11 
shows  the  change  in  magnetostriction  for  a  sample  of  68  Permalloy 
(68%  Ni,  32%  Fe) .  Figure  III. 11  illustrates  how  the 
magnetostriction  of  the  material  was  increased  in  a  direction 
perpendicular  to  the  annealing  field,  and  how  it  decreased  in  a 
direction  aligned  with  the  external  field.  The  graph's  abscissa 
(magnetization)  refers  to  the  magnetic  dipole  moment  per  volume 
of  the  sample  (2:217)  which  is  proportional  to  the  magnetic  field 
flux  density  and  has  the  same  units  of  gauss. 

Although  not  addressed  in  the  literature,  it  is  postulated 
that  the  geometrical  configuration  of  the  magnetostr ict ive 
material  may  affect  its  bulk  magnetostrict ive  property.  That  is, 
the  magnetic  domains  may  rotate  more  easily  if  they  are  not  in 
contact  with  other  domains  that  are  also  attempting  to  move  under 
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the  influence  of  an  external  magnetic  field.  Therefore,  by 
investigating  the  influence  of  the  magnetostriction  of  domain¬ 
sized  stripes  of  magnetostrictive  material,  the  magnetostriction 
of  the  collection  may  be  greater  than  that  of  a  uniform  coating. 


04  - - - - 

- - -  - 

M 

20  - 

Transverse  Field 

a 

je 

16  ■ 

n 

t 

o 

12  - 

t 

r 

l 

N  Field 

c 

t 

i 

8  - 

o 

n 

1 

Longitudinal  (  it  Id 

- 

-r. 

IN  10 

) 

()  — — 

-  -  •  ■  — 

0  2 

•1  6  H  10  12  11 

Majinot  r/.at  ion  (hi  i  \  io  i 

Figure  III. 11.  Effect  of  Magnetic  Annealing  on  the 

Magnetostriction  Effect  (20:638) 


III-21 


Fiber  Optic  Cables 

A  description  of  fiber  optic  cables,  in  addition  to  the 
discussion  of  magnetostrict ive  material,  completes  the  analysis 
of  the  components  of  the  sensor  arm.  The  composition  of  the 
fiber  optic  cable  directly  affects  the  detection  mechanism  of  the 
sensor  by  establishing  how  the  light  propagates  in  the 
interferometer . 

The  fiber  optic  cables  used  throughout  the  interferometer 
are  York  Technology,  Model  SM600,  single-mode  fibers.  These 
fibers  are  composed  of  four  concentric  circular  material  layers 
as  shown  in  Figure  III. 12.  The  inner  core  is  composed  of  Ge- 


doped  silicon  dioxide  (Ge:Si02)  which  has  an  index  of  refraction 
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York  Technology,  Model  SM600,  Single-mode  Fiber 
Optic  Cable  Cross  Section  (24) 


1 1 1 -2  2 


(n)  of  1.461  and  a  radius  of  1.75  microns.  The  outer  glass 
cladding  is  composed  of  two  parts:  an  inner  cladding  and  an  outer 
cladding.  These  two  cladding  layers  are  used  to  achieve  a 
continuously-graded  profile  of  the  index  of  refraction.  The 
inner  cladding  has  a  radius  of  25  microns  and  is  a  phosphorous 
fluorine-doped  silicon  dioxide  (PF:Si02)  material,  which  has  an 
index  of  refraction  of  1.457.  The  outer  cladding  is  undoped 
silicon  dioxide  (Si02)  ,  which  has  an  index  of  refraction  of  1.456 
and  a  radius  of  62.5  microns.  The  final  layer  is  a  jacket  of 
mode  absorbing  ultraviolet  cured  acrylate  which  has  a  radius  of 
110  microns.  The  jacketing  material  provides  mechanical  strength 
and  absorbs  undesirable  light  which  may  propagate  in  the  cladding 
layers  (24) . 

Under  the  influence  of  an  external  magnetic  field,  the 
magnetostrictive  material  covering  the  sensor  arm  will  apply  a 
force  on  the  fiber  optic  cable  to  change  the  cable's  length. 
Regardless  of  the  source  of  the  force,  the  length  change  of  the 
optical  fiber  has  been  shown  to  be  given  by  (25:102): 

=  FLn  (17)  , 

EA' 

where  aL  =  change  in  length 
Ln  =  initial  length 

F  =  magnitude  of  the  applied  force 
A  =  cross  sectional  area  of  the  fiber 
E  =  Young's  modulus  of  the  fiber. 
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This  expression  indicates  that  by  reducing  the  cross-sectional 
area  of  the  fiber,  an  increase  in  the  change  of  length  will 
result  when  the  other  variables  remain  constant.  The  cross- 
sectional  area  of  the  fiber  can  be  reduced  by  thinning  the 
cladding  of  the  fiber  optic  cable. 

Cylindrical  Waveguide  Propagation 

The  optical  fiber  previously  described  will  serve  as  a 
waveguide  for  laser  radiation  in  the  sensor.  The  choice  of  an 
optical  source  and  optical  fiber  relative  to  their  influence  on 
the  sensor's  performance  will  be  determined  by  examining 
cylindrical  waveguide  propagation. 

Method  of  Analysis.  This  section  describes  the 
electromagnetic  wave  propagation  phenomena  in  an  optical  fiber. 
The  analysis  begins  by  assuming  the  general  case  cf  a  time- 
varying  (propagating)  electromagnetic  wave.  These  solutions  are 
then  substituted  into  the  cylindrical  coordinate  system 
representation  of  Maxwell's  equations.  All  components  of  the 
electric  and  magnetic  fields  are  then  found  in  terms  of  the 
electric  and  magnetic  field's  z-components .  The  wave  equation  is 
then  derived  and  solved  using  the  separation  of  variables 
technique  to  yield  the  z-components.  Finally,  the  specific 
geometry  and  corresponding  boundary  conditions  of  an  optical 
fiber  are  used  to  determine  the  constants  generated  in  the 
solution  of  the  differential  equations.  These  solutions  provide 
insight  into  the  propagation  of  the  allowed  modes. 

Maxwell's  Equations.  For  nonmagnetic  and  nonconducting 


1 1 1-2  4 


materials  with  no  enclosed  charge,  the  set  of  Maxwell's  equations 
which  relate  the  electric  field  (E) ,  magnetic  field  (H) , 
permittivity  (e  =  £oer),  and  permeability  (p=po)  may  be  represented 
as  (bold  characters  indicate  vector  quantities)  (2:374): 


V 

x  H  =  edE 

(18) 

dt 

7 

x  E  =  ~udH 

(19) 

dt 

V 

•  E  =  0 

(20) 

V 

•  H  =  0 

(21)  . 

(The  italic  d  (d)  denotes  a  partial  derivative  in  this  analysis). 
The  electric  and  magnetic  fields  of  the  laser  radiation  are 
assumed  to  be  propagating  along  the  z-direction  of  the  optical 
fiber  and  are  of  the  form  (26:55): 

E  =  E0(r,*)exp[j  (wt-/?z)  ]  (22) 

H  =  Hc(r,  e )  exp[  j  (u-t-/3z)  ]  ,  (23) 

where  w  =  angular  frequency 

(3  =  axial  phase  constant. 

Substituting  the  assumed  form  of  the  electric  and  magnetic  fields 
(equations  22  and  23)  into  the  set  of  Maxwell's  general  equations 
(equations  18,  19,  20,  and  21)  yields  (2:374): 


t  x  H  -  ji.eE  (2  4  ) 

V  x  E  =  -ju-pH  (25) 

7  •  E  =  0  (26) 

V  •  H  =  0  (27). 
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E  and  H  Field  Components.  Using  the  expressions  for  the 
curl  and  divergence  in  cylindrical  coordinates,  equations  24  and 
25  can  be  expanded  as: 


dUl  -  dH. 
rdf  dz  ' 


+  ©rdHr  -  dHj  +  zlfdrH..  -  dH, 
[dz  dr  j  r  |_  dr'  d f 


and 


r 

r 

~  dE- ! 
dz  e\ 

+  e 

~dEr 

~  dEl 

+  zl 

rdrE. 

~  £LEr" 

r  df 

d  z 

dr  j 

r 

dr 

df 

j  -  c  E 


=  -  j  u-fiH 


(28) 


(29)  . 


Each  component  of  equations  28  and  29  is  a  separate  equation. 
Evaluating  the  partial  derivatives  where  possible,  the  six 


.ions 

become 

(27 : 

49)  : 

dH2  ■ 

ii 

Q. 

1 

dHz 

+  D0K.  -  ju-eEr 

(30) 

rdf 

dz 

rdf 

dHr  - 

^h2  =  - 

■j/5Hr 

-  dH  =  joeE, 

(31) 

dz 

dr 

dr 

drHg 

-  dHr 

=  jw 

£EZ 

(32) 

rdr 

rdf 

_dE: 

-  dE_  =- 

dEz 

+  j£E_  = 

(33) 

rdf 

dz 

rdf 

—r  ~ 

la 

lw 

II 

1 

■j^Er 

-  dE2  =  -ju>/iHs 

(34) 

dz 

dr 

dr 

drE. 

-  dE,  = 

=  -j- 

uh’z 

(35) 

rdr 

rdf 

:r, 

,  Hr ,  H. 

components  car.  be  conveniently  expressed 

ir. 

;  of  the  Ez  a 

nd  H. 

t  components.  That  is,  from  equation 

30,  E 

can  be  expressed  as: 

Er  =  -  ■  j  dH,  +  £Hr 
u'tr  df  u-e 

Correspondingly,  from  equation  34,  H„  can  be  expressed  as: 

h,  =  n.r  -  jm, 

u'H  u;fjdr 


(36)  . 


(37)  . 


Substituting  equation  37  into  equation  36  yields: 
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(38)  . 


Er  -  -  - ^ 

uer  de 


im7  +  _1 

u)  £ 


1  -  -4 


2  1  =  - 


CM 


£Er  -  _OdE2 

lo\i  u^idr 


diM,  -  _Jj3dE 


u>erd0 


v  iiedr 


Er  =  -  _ j  dH^  -  _  ifl  dEr 

«er [1-flV  (c^e/i)  ]dtf  u>2/ie  [  1-/32/  (777m)  ]dr 

The  following  identities  can  then  be  used  to  simplify  the  above 

expression  (27:50): 


(35) 


k2  =  <7e  m  (wave  number  squared) 

/S'  =  k2n2  -  (3  (transverse  wave  number) 

n  =  (ur£r)"  =  (index  of  refraction). 

Therefore,  equation  38  now  can  be  written  as: 

Er  =  _nlkiKfH2  +  £dE2] 
rdS  dr  J 

The  other  components  can  be  derived  in  a  similar  manner,  and  they 
are  given  as  (27:50): 


(40) 


H. 


— 

P'  \ 

[rdfl 

-i  , 

r^uH 

(3 '2 

[  dr 

\iiih 

idE7  -  y^id H 


a  r 


rdt 


(41) 


(42) 


(43) 


/S  ’  l_rdf  Ur  j 

Wave  Equation.  All  the  components  of  the  electric  and 
magnetic  fields  can  be  specified  cr.ce  E,  and  H,  are  known.  The 
z-components  can  be  determined  using  the  wave  equation.  To 
derive  the  wave  equation  for  the  electric  field,  equations 
43  are  substituted  into  equation  32  which  yields: 


j  E 


_dr  It 
rdr [p  ' 


+  ulriEj  j  -  d 

rdf>  ur  j  J 


PdH.  -  e  c.E. 
rdO[p'c{  dr  r at 


4-1=2? 

J  /P  <  ^ 


III-27 


jweE2 

=  -_.j  a/ h2  - 

i^edrdE.  4- 

rp  ,2drd(? 

P '  2rdr2 

=  -l^irf_!~rdE  1 
p  ' £ rdr L  dr 

1  - 

j3'Vdr 

jweE2 

=  -iu’edE.  -  iwerd2E_  - 

p ' 2rdr"  p 

i '  2rdr2*  p 

d2  Ez  + 
d? 

dE  +  d2E  4 

Z  *5  2  p 

rdr  rdr 

■  P,2ez  =  o 

-i e  d2V 


(44)  . 

Equation  44  is  the  wave  equation  for  the  electric  field.  The 
wave  equation  for  the  magnetic  field  is  derived  in  the  same 
manner,  and  it  is  given  by  (27:51): 

0  (45) . 


d2H  +  dH  +  d2H  +  S'2H  = 

- PZ  - Z  - T — Z  ?  ^  Z 

dr  rdr  r  dr 


The  wave  equation  in  equations  44  and  45  can  be  solved  using 
the  method  of  separation  of  variables.  That  is,  a  solution  for 
the  electric  field  can  be  assumed  to  be  of  the  form  (27:51): 


Ez  =  Rz(r)e*<ff)  =  Re  (46)  • 

To  solve  for  each  component,  equation  46  is  substituted  into  the 

wave  equation  (equation  44)  which  yields: 

d2  (R9 )  +  d  (RQ )  +  d-?  (RQ)  +  P  '  2R9  =  0 
dr  rdr  r  2d$2 

9d2R  +  9dR  +  Rd26 .  +  /?  ' 2  =  0 
dr2~  rdr  r2d 5 ' 

d2  R  +  _ dR  +  d2e  +  P  '  2r2  =  0  (47)  . 

Rdr2  Rrdr  9df2 

Equation  47  can  be  separated  into  two  equations  that  are  each  a 
separate  function  of  R  and  8.  That  is,  the  only  way  for  the  sum 
of  the  two  functions  to  equal  a  constant,  is  for  each  individual 
function  to  be  equal  to  a  constant.  Since  equation  47  is  equal 
to  zero,  the  sum  of  these  two  constants  must  be  zero.  Therefore, 
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the  two  constants  are  equal,  but  have  different  signs.  By 

specifying  these  arbitrary  constants  as  N2  and  -N2,  the  two 

equations  are  given  by: 

r 2d2R  +  rdR  +  (3 '  2r2  =  N2 
Rdr2  Rdr 

d20  =  -N2 

edT2 


(48) 

(49) 


Equations  48  and  49  are  rewritten  as  follows  (27:51): 

d2R  +  _dR  +  R  f/3  ' 2  -  N2]  =  0  (50) 

dr 2  rdr  [  PJ 

d©2  +  ©N2  =0  (51) 

7P 


(52) 


The  closed  form  solution  of  equation  51  is  given  by  (27:51): 
f  cos(N0  +  <p) 

ez(n  =  \  . 

[  sin(Nf  +  <p ) 

where  0  is  an  arbitrary  phase  constant  (assumed  to  be  zero  for 
simplicity).  Equation  50  is  a  specific  form  of  Bessel’s 
differentia]  equation,  and  it  has  a  solution  given  by  (27:51): 


Mr)  - 


AJK(P'r)  +  A'YN(/?'r)  when  (3'  is  real 


(53) 


CKN(|/?'jr)  +C '  1N  (  j  f3  '  j  r)  when  (3  '  is  imaginary  (54) 


where : 


N  =  the  azimuthal  mode  number 
A,  A',  C,  C'  =  arbitrary  constants 

JK ,  \'K  =  Nth-order  Bessel  functions  of  the  first  and 
second  kind,  respectively 

Kn,  In  =  Nth-order  modified  Bessel  functions  of  the  first 
and  second  kind,  respectively. 
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The  Bessel  functions  are  illustrated  in  Figures  III. 13  and 
III. 14,  and  the  modified  Bessel  function  are  shown  in  Figures 
III. 15  and  III. 16. 

Boundary  Conditions.  To  determine  the  constants  in  the 
solution  to  the  differential  equations,  the  specific  geometry  of 
an  optical  fiber  is  now  applied  to  the  equations  for  R  and  8.  As 
discussed  earlier,  an  optical  fiber  is  composed  of  two  concentric 
circular  layers  of  dielectric  material.  The  inner  core  has  a 
radius  (a)  and  an  index  of  refraction  (n,).  The  outer  cladding 
has  a  radius  (b)  and  an  index  of  refraction  (n2)  .  The 
corresponding  boundary  conditions  for  the  electric  and  magnetic 
fields  are  (27:56): 


p  core 

L2 

£  cladding 

(55) 

—  core  _ 

r 

£  cladding 

f 

(56) 

E  core 

P  cladding 

-  e2rr 

(57) 

Hzcore  - 

,,  cladding 

HZ 

(58) 

H  -re  = 

„  cladding 

(59) 

Hrc°re  s 

jj,  cladding 

=  Hr 

(60)  . 

The  boundary  conditions  also  dictate  that  the  azimuthal  mode 
number,  N,  and  the  axial  phase  constant,  0,  must  be  consistent 
for  the  fields  in  both  regions  (27:53).  Therefore,  the  solutions 
in  the  cladding  and  core  regions  will  be  examined  separately,  and 
then  they  will  be  reconciled  at  the  core-cladding  interface. 

For  the  cladding  region,  the  R2(r)  solution  must  decay  in  an 
exponential  fashion  in  order  for  light  to  propagate  in  the  core. 
This  requirement  means  that  equation  53  cannot  be  a  solution  for 
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Figure  III.  14.  Nth-Order  Bessel  Functions  of  the  Second  Kind 

(28:391) 
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1  o 


Figure  III. 16.  Nth-Order  Modified  Bessel  Functions  of  the 

Second  Kind  (28:417) 
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the  electric  field  in  the  cladding.  Consequently,  the  constant 
C'  in  equation  54  must  also  be  zero  since  In  diverges  for  large 
values  of  r  (27:53). 

For  the  core  region,  the  solution  cannot  contain  YN  or  K^, 
since  these  functions  diverge  on-axis  when  r  equals  zero.  The 
requirement  for  the  electric  field  to  be  smoothly  connected  at 
the  interface  prevents  the  solution  from  containing  IK.  As  shown 
in  the  graphical  representations,  IN'/1N  (the  '  refers  to 
differentiation  with  respect  to  r)  is  always  positive,  and  KN'/K*, 
(the  cladding  solution  function)  is  always  negative  (27:53). 

Therefore,  the  solution  in  the  core  is  proportional  to  JN, 
and  the  solution  in  the  cladding  is  proportional  to  KN .  The 
boundary  conditions  also  indicate  that  f3 '  must  be  real  in  the 
core  and  imaginary  in  the  cladding.  Therefore,  based  on  the 
definition  of  £> '  (equation  39),  (3  is  bounded  as  follows  (27:53): 

kn2  <  0  <  knr 

The  analysis  for  the  magnetic  field  component  follows  the 
analysis  outlined  for  the  electric  field.  Based  upon  the 
previous  arguments,  the  z-components  of  the  electric  and  magnetic 
fields  are  expressed  as  (27:54): 


E,  =  AJ^ (ur/a) sin (Nf ) 

(61) 

H;  =  BJ^ (ur/a) cos (N? ) 

(62) 

the  core,  and 

Ea  =  CKn (wr/a) sin (Nf ) 

(63) 

Hz  =  DKn  (wr/a)  cos  (Nfl  ) 

(64) 

the  cladding, 
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where 


u  =  p\a  =  (k2n,2  -  p2)'ha  (normalized  transverse  phase 

constant  in  the  core)  (65) 

w  =|/3'2|a  =  (P2  -  k2n22pa  (normalized  transverse  phase 

constant  in  the  cladding)  (66). 

The  choice  of  the  sine  or  cosine  function  in  equations  61  through 
64  was  made  to  facilitate  the  hybrid  mode  solution.  With  E2  and 
Hz  specified,  the  other  components  can  be  derived  by  substituting 
equations  61,  62,  63  and  64  into  equations  40,  41,  42,  and  43. 

As  an  example,  Er  in  the  core  is  derived  as  follows: 


a  r 


Er  = 

P77 


Er  ~  — — -jr 

J 72 


E.  = 


Er  = 


uMirilk  +  £dE2 
r  d$  dr 

u'Ud  (BJ^  (ur/a)  cos  (NS  ) )  +  p_d  A  J  N(ur/a)  s  i  n  ( N  6  )] 
rdS  (  J  dr  [ 

lcuBJ  (ur/a)Nsin(N£)  +  pkuJ  1  (ur/a)  sin  (Nf? ) 


(40) 


-AifiJ,,'  fur/a)  +  BiwuNJt.(ur/a) 
(u/a)  (u/apr 


cos (N£ ) 


(67) 


The  other  components  are  given  as  (27:54-55): 
Es  =  f -Aj /3NJn  ( ur / a )  +  Biu-uJ,/  (ur/a) 


( 68 ; 


;  (u/apr 

(u/a) 

1 

Hr 

= 

'Ai^NJ,  (ur/a)  - 

Bij3J(i'  (ur/a) 1 

cos (N  r ) 

(60) 

I  (u/apr 

(u/a)  J 

U 

i 

1 

r  „  > 

|  i  w  -  -  ’ \  '  *  /  a  . 

'a) 

s  i  n  ( N ::  ) 

(70) 

• 

L  (u/a) 

(u/apr 

the 

core,  and 

Er 

= 

j" _C.jp_kNi_ ( wr /a±  - 

Diu-uNK  (wr/a)  ]  sin  (Nr  ) 

(71) 

[  (w/a) 

(w/a )  r 

i 

j 

E. 

t* 

[  ._Cj/3NKw(wr/a)  - 

DiuuJN'  (wr/a) 

1  cos (N  5 ) 

(72) 

(w/a )  r 


(w/a) 
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«r  = 


Hfi  = 


'-Ci^e-,NKt,  (wr/a)  +  DigK^1  fwr/a)lcosfN^ 
(w/apr  (w/a) 

'Cioe-^ l  (wr/a)  -  Pi (?NKU  1  fwr/a )  ]  sin (Ne  ) 
(w/a)  (w/apr 


(73) 

(74) 


in  the  cladding. 

A  set  of  equations  for  the  constants  A,  B,  C,  and  D  is 
obtained  by  setting  r  equal  to  *'3"  in  the  expressions  for  the 
electric  and  magnetic  field  components  (equations  61  through  74), 
and  substituting  these  equations  into  the  boundary  condition 
equations  (equations  55  through  60) . 

As  an  example,  starting  with  equation  55: 


core  _  £  cladding 


(55)  . 


Substituting  equations  61  and  63  into  equation  55  yields: 

AJn  (ur/a) sin (N8 )  =  CKN  (wr/a ) sin (Ntf ) . 

Setting  r=a  in  the  above  expression  yields: 

AJN(u)  -  CKn(w)  =  0  (75)  . 

The  same  method  produces  the  following  equations  (27:57): 

BJn(u)  -  DK„(w)  =  0  (76) 

Aj<3NJtJ(u)  -  BiwuJ,,'  (u)  +  CjflNKt.fw)  -  Dj.-iiK,/  (w)  =0  (77) 

(u/apa  (u/a)  (w/apa  (w/a) 

1  (u)  -  Bj/5NJh  (u)  Cju  e-,K(j  1  (w)  -  DjfiNK,.  (wj  =  0  (78)  . 

(u/a)  (u/apa  (w/a)  (w/a)  a 

A  normalized  set  of  values  for  the  constants  can  be  derived  by- 

setting  A  equal  to  one  and  using  equations  75  through  78.  From 

equation  75,  C  is  found  to  be: 


C  =  J*lul 
Kn(w) 

From  equation  76,  D  in  terms  of  B  is  given  by: 


(79)  . 
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(80)  . 


D  =J,Jfu)B 
K^(w) 

To  solve  for  B,  equations  79  and  80  are  substituted  into  equation 
77  to  yield: 

jflNJ,,  (U)  -  BioMJ,/  (U)  +  JN(u)Ci/?NK  (w)  -  BC^Ullo/iK,1  (wj  =  0. 

(u/a)V  (u/a)  Kn(w)  (w/a)2a  KN(w)(w/a) 

Solving  for  B  gives: 


BN 


1  +  1 
U2  W2 


B  = 


(81)  . 


:VJul  +  K.Hvl 

[uJN(u)  WKN(W) 


Correspondingly,  substituting  equation  81  into  80  yields  D: 


mNXul 

o)^Kn(w)  [u‘  w 


1  +  1 
T2  .12 


D  = 


(82) 


JN.,.l,u.l.  +  iy_iwi 

UJN(U)  Kn(w)w 


Finally,  all  the  components,  including  the  normalized  constants, 
are  now  specified  for  the  electric  and  magnetic  fields  in  the 
optical  fiber. 

Allowed  Modes.  The  set  of  constant  equations  (equations  75 
through  78)  limit  the  modes  chat  are  will  propagate  in  the  fiber. 
Equations  75  through  78  can  be  written  in  a  matrix  form  as 
(27:57): 


[A] 

[M]  X  I B 

j  C 

ld 


0, 


where 
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JN(u) 

0 

-K,(w) 

0 

0 

J*(u) 

0 

-Kn(w) 

M  = 

(u/aj^a 

(U) 

(u/a) 

i  BNKn  ( w ) 
(w/a)za 

(w/a) 

l^JVlul 

L  (u/a) 

-jfiNJw(u) 

(u/a)  a 

jwi;>V  (w) 
(w/a) 

-1£NKn(w) 
(w/a) 2a 

For  a  solution  to  exist  for  the  A,  B,  C,  and  D  constants,  the 

determinant  of  the  matrix  M  must  be  zero.  Evaluating  the 

determinant  produces  the  following  equation  (27:57): 

J n-'-Lu).  +  K^liXL]  +  V  (w)' 

uJN(u)  wKn(w)  J  Le2uJN(u)  wK^w) 

=  Nzfl  +  1  ]  (e  +  1}  (83)  . 

(u  w  J  [e2u  w  J 

Equation  83  places  restrictions  on  the  allowed  modes.  These 
restrictions  are  examined  by  the  various  mode  classifications. 
That  is,  the  boundary  conditions  at  the  core-cladding  interface 
dictate  that  TM  (transverse  magnetic)  and  TE  (transverse 
electric)  modes  can  only  exist  when  N  equals  zero  in  the 
expressions  for  the  electric  and  magnetic  fields  in  the  core  and 
cladding  (equations  61  through  74)  (27:54).  When  N  is  not  zero, 

the  boundary  conditions  can  be  satisfied  by  a  linear  combination 
of  TM  and  TE  modes  which  are  labeled  HE  and  EH  modes  (hybrid 
modes) .  TM  modes  occur  when  Hz  and  N  are  zero,  which  from 
equations  62  and  64,  requires  that  the  constants  B  and  D  are 
zero . 

For  TM  modes,  equation  83  is  simplified  to: 

f.lSV_Lu 1  +  Ko'.jw).  =  0  (84)  . 

e2uJ0(u)  wK0(w) 
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For  TE  modes,  Ez  and  N  are  zero  which  requires  that  the 

constants  A  and  C  in  equations  61  and  63  are  zero.  Consequently, 

equation  83  for  TE  modes  can  be  simplified  to: 

V  (u)  +  K0 '  (w)  =  0  (85)  . 

uJ0(u)  wK0(w) 

To  simplify  equation  83  for  the  hybrid  modes,  it  is 
necessary  to  invoke  the  weakly-guided  approximation,  which 
assumes  that  e,  is  approximately  equal  to  e2.  With  this 
approximation,  equation  64  for  the  hybrid  modes  can  be  simplified 
to: 

JV-lPl  +  Kh-'JwI  =  ±N2(l  +  1  |  (86)  . 

uJN(u)  wKn(w)  (u7  w2} 

Using  the  weakly-guided  approximation  for  the  TM  modes  simplifies 
equation  84  into  equation  85.  Equation  85  can  be  simplified 
further  using  the  following  identities  (27:59): 

J0' (u)  =  ~J, (u) 

V  (W)  =  -K,(w) 

v2  =  u2  +  w2  =  k2n.2a22t  (normalized  frequency  squared)  (87) 

^  =  n,2  -  n:2  (88)  . 

2n7 

Equation  85  is  now  of  the  form: 

=  _K,(tyg-u2I^l_ . (39). 

uJ0(u)  (v  -u2)  ‘Kc(  (v  -u2)  ) 

Equation  89  can  be  solved  graphically  to  illustrate  the  allowed 
TM  or  TE  modes.  The  v-parameter  is  a  number  which  can  be 
calculated  for  a  particular  fiber  optic  cable  and  excitation 
source  wavelength.  For  example,  if  v  is  set  equal  to  seven, 
Figure  III. 17  yields  a  graphical  solution  when  N  equals  zero.  In 
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Figure  III. 17,  the  right-hand  side  of  equation  89  (T,)  is  an 
exponentially  growing  function,  and  the  tangent-like  function  is 
the  left-hand  side  of  equation  89  ( X1 )  .  The  intersections  of  the 
two  functions  satisfy  equation  87  for  the  corresponding  value  of 
u.  For  a  particular  solution  set  of  u  and  v,  a  propagation 
constant  f3  for  the  allowed  mode  can  be  calculated.  To  generate 
similar  graphical  solutions  for  the  hybrid  modes,  equation  88  is 
simplified  using  the  following  identities  (27:267): 


TMC1,  or 


TM  or 

Oil 


Normaliz'd  1  i  ; »  1 1  s  \  <  1  • <  1  ’  1  \ ;  \ 


igure  III.  17.  Allowed  TK  or  TE  Modes 


Equation  90  and  91  are  substituted  into  equation  86  to  yield  two 
equations  due  to  the  positive  and  negative  signs  in  equation  86; 
that  is: 

=  -JW  (v?-uV)  (92) 

uJN(u)  (v2-u2)\(  (v2-u2)'/l) 

and 

=  K^llyWhl)  (93). 

uJN(u)  (v2-u2)\((v2-u2)H 

Equation  92  describes  the  allowed  EH  inodes,  and  equation  93 
describes  the  allowed  HE  inodes.  Graphical  solutions  to  equations 
92  and  93  for  N  equal  to  one  are  shown  in  Figures  III. 18  and 
III. 19.  The  axis  labels  for  Figures  III. 18  and  III. 19  are  as 
follows:  r2  refers  to  the  left-hand  side  of  equation  92,  X2 

refers  to  the  right-hand  side  of  equation  92,  r3  refers  to  the 
left-hand  side  of  equation  93,  and  X3  refers  to  the  right-hand 
side  of  equation  93. 

Observations .  As  seen  in  Figures  111.17  and  III. 18,  there 
are  values  of  the  normalized  frequency  for  which  there  are  no 
allowed  TE ,  TM,  or  EH  nodes.  However,  as  shown  in  Figure  III. 19, 
there  is  always  an  allowed  HE  node.  Only  the  HE.,  node  will 
propagate  as  long  as  the  normalized  frequency  is  sufficiently 
snail  to  prevent  any  allowed  TMon  or  TEp,  modes  from  propagating. 
From  equation  89,  this  critical  value  of  the  normalized  frequency 
occurs  at  the  first  zero  of  J  ,  which  occurs  at  a  value  of  2.405. 
Correspondingly,  this  value  of  the  normalized  frequency 
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Figure  III. IS.  Allowed  EH  Modes 


Figure  III. 19.  Allowed  HE  Modes 
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constrains  the  values  of  the  fiber's  core  radius  (a),  core  index 
of  refraction  (n,)  ,  cladding  index  of  refraction  (n2)  ,  and 
excitation  wavelength  (t).  More  explicitly,,  equation  89  is 
expanded  to  show  the  relationships  between  these  parameters. 

That  is: 

v  =  kn1a2%AV’  =  2V2mn1aAVi  <  2.405  (94) 

X 

Assuming  that  the  index  of  refraction  difference  between  the  core 
and  inner  cladding  and  the  value  of  the  fiber  core's  index  of 
refraction  are  limited  by  manufacturing  considerat ions ,  the  two 
parameters  of  interest  are  the  core  size  and  excitation 
wavelength.  The  core  size  should  be  maximized  to  facilitate  the 
source-to-f iber  coupling.  Therefore,  the  excitation  wavelength 
should  also  be  maximized.  However,  there  is  an  upper  limit  on 
the  excitation  wavelength  in  the  near  infrared  spectrum  due  to 
absorption  in  the  fiber  itself  (26:87).  In  addition,  there  are 
further  considerations  concerning  the  choice  of  the  optical 
source's  wavelength  that  are  discussed  in  the  next  section. 

Specific  Calculation.  The  electric  and  magnetic  fields 
propagating  in  an  optical  fiber  are  now  calculated  for  the 
particular  optical  source  and  fiber  that  was  used  in  this 
research  effort.  After  the  phase  constants  are  calculated,  the 
A,  B,  C,  and  D  constants  are  determined.  The  final  step  is  to 
calculate  and  plot  the  electric  and  magnetic  fields. 

The  normalized  frequency  (v)  is  calculated  to  be  1.87  (unit¬ 
less)  using  equation  87  and  the  fiber  optic  cable  dimensions  for 
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the  York  Technology,  Model  SM600,  fiber  that  was  described  in  the 
previous  section.  To  determine  the  normalized  transverse  phase 
constant  in  the  core  (u) ,  a  BASIC  computer  program  was  written 
(included  in  Appendix  C)  to  solve  equation  89.  Equation  89  was 
also  solved  graphically  as  shown  in  Figure  III. 20.  As  a  result, 
u  was  found  to  be  1.48  (unit-less).  The  axial  phase  constant  (/3) 
was  next  evaluated  to  be  1.448  x  107  m'1  from  equation  65. 

Additionally,  the  u,  v,  and  /S  constants  were  evaluated  in 
another  BASIC  computer  program  (included  in  Appendix  C)  which 
also  solved  for  the  A,  B,  C,  and  D  constants.  This  computer 
program  used  numerical  series  approximations  for  the  Bessel 
functions  and  modified  Bessel  functions.  The  constants  A,  B,  C, 
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Figure  III. 20.  Single-mode  Operation  of  York  Technology,  Model 

SM600,  Optical  Fiber  With  a  0.633  urn  Laser 
Source 
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and  D  were  found  to  be  1,  -3.832  x  10'3,  1.178,  and  -4.509  x  10'3 
(unit-less),  respectively.  The  same  computer  program  then  solved 
for  the  electric  and  magnetic  fields  in  the  fiber.  The  electric 
field  was  found  to  be  linearly  polarized  in  the  negative  y- 
direction,  and  the  magnetic  field  was  linearly  polarized  in  the 
positive  x-direction.  The  components  of  electric  and  magnetic 
fields  in  the  z-direction  were  found  to  be  negligible  in 
magnitude  compared  to  the  components  in  the  x-y  plane.  Figures 
III. 21  and  III. 22  illustrate  the  electric  and  magnetic  fields 
propagating  in  the  fiber  optic  cable.  The  electric  and  magnetic 
fields  have  decayed  significantly  after  only  a  short  distance 
into  the  cladding,  which  extends  to  over  35  times  the  radius  of 
the  core. 

To  summarize,  solving  for  the  electric  and  magnetic  fields 
that  propagated  in  the  particular  fiber  optic  cable  used  in  this 
research  project  contributed  information  concerning  the  sensor's 
design.  That  is,  chere  is  a  trade-off  concerning  fiber  core 
radius  and  optical  source  wavelength,  with  maximum  wavelength  and 
minimum  core  radius  desireable  for  single-mode  operation,  with  a 
maximum  core  radius  preferred  for  source-to-f iber  coupling.  The 
analysis  also  showed  that  thinning  the  cladding  would  net 
significantly  affect  the  propagation  of  the  electromagnetic  wave 
since  the  wave  has  decayed  after  a  short  distance  into  the 
cladding . 

Interferometer  Operation 

Additional  information  concerning  the  sensor's  design  can  be 
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Figure  III. 21.  Electric  Field  Propagating  in  York  Technology, 

Model  SMGOO,  Optical  Fiber  with  a  0.633  ur.  Laser 
Source 


Figure  111.22.  Magnetic  Field  Propagating  in  York  Technolo 

Model  SM600,  Optical  Fiber  with  a  0.633  ur  La 
Source 
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gained  by  examining  the  behavior  of  interferometers.  A  fiber 
optic  implementation  of  the  Mach-Zehnder  (MZ)  interferometer  can 
be  used  to  measure  small  changes  in  a  fiber's  length,  such  as 
those  induced  by  a  magnetostrict ive  material  coating.  This 
section  explains  how  the  small  change  in  length  in  one  arm  of  a 
MZ  interferometer  can  be  measured  as  a  corresponding  change  in 
the  phase  of  light  propagating  in  that  arm.  This  phase  change  is 
then  converted  into  a  detectable  amplitude  change.  The 
components  of  a  MZ  interferometer  are  illustrated  in  Figure  1.1. 
The  section  concludes  with  a  discussion  of  the  significant 
factors  affecting  magnetostr ictive  interferometer  performance. 

As  shown  in  the  previous  analysis,  the  only  mode  propagating 
in  single-mode  fiber  optic  cables  is  the  HEn  mode,  which  is 
linearly  polarized.  Using  an  equivalent  form  of  equation  27,  the 
electric  field  in  the  reference  and  sensor  arms  can  be  modeled  as 
(29 : 277)  : 

Er(z,t)  =  Erocos(/?rz-_t--;  J  (95) 

Es(z,t)  =  E&0cos(/3sz-wt+4s)  ,  (96) 

where  Er  ( z , t  =  electric  field  in  the  reference  arm 
Es(z,t)  =  electric  field  in  the  sensor  arm 
Erc  =  reference  arm  electric  field  vector  at  time  t, 
and  position  z 

Esc  =  sensor  arm  electric  field  vector  at  time  t, 
and  position  z 
/?  =  axial  phase  constant 

fr  =  initial  phase  constant  in  the  reference  arm 
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£s  =  initial  phase  constant  in  the  sensor  arm. 

For  simplicity,  the  electric  field  is  assumed  to  be  equal  in 
magnitude  and  parallel  in  each  arm,  with  equal  initial  phase 
constants.  While  the  initial  phase  constants  may  not  be  equal 
due  to  the  fact  that  the  sensor  arm  and  reference  arm  are  not 
exactly  the  same  length,  this  assumption  merely  adds  a  constant 
in  the  final  expression.  The  assumption  that  the  electric  fields 
are  equal  in  magnitude  is  valid  since  the  optical  power  in  each 
arm  of  the  sensor  can  be  adjusted  until  they  are  equal.  If  the 
electric  fields  in  each  arm  were  not  parallel,  but  constantly  at 
an  angle  with  respect  to  each  other,  this  assumption  would  again 
only  add  another  constant  to  the  final  result.  However,  if  the 
orientation  of  the  two  electric  fields  varied  with  time,  then  the 
final  result  would  be  modulated  by  a  time-varying  factor.  Since 
the  behavior  of  the  electric  fields  in  a  fused  coupler  has  not 
been  modeled,  it  will  be  assumed  that  the  fields  are  constant 
with  respect  to  each  other.  The  mathematical  representation  of 
these  assumptions  is  stated  as: 

E-o  =  Eso  =  Ec  (9?) 

*r  ^  =  i  (96)  • 

Due  to  the  high  frequencies  of  optical  signals  (-10“  Hr), 
light  detectors  do  not  detect  electric  fields  (29:45).  Instead, 
optical  detectors  sense  irradiance  which  is  given  by  (29:277): 

I  =  e  v<E2>  (99) 

where  I  =  irradiance  (W/cm2) 

e  =  permittivity  of  the  propagation  media 
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v  =  speed  of  the  propagating  optical  energy 

E2  =  <|E2|>  =  < | E - E [ >  =  time-average  value  of  the  magnitude 

of  the  square  of  the  electric 
field. 

The  time  average  of  a  function  f(t)  is  given  by  (29:58): 


<f(t)>  =  1 
T 


t+T 


f (t ' )dt' 


where  T  =  time  interval 

t'  =  dummy  integration  variable. 

The  high  frequency  of  the  electromagnetic  wave  compared  to  the 
frequency  response  of  light  detectors  indicates  that  the  time 
interval  (T)  is  much  greater  than  the  period  of  the 
electromagnetic  wave.  The  square  of  the  electric  field  for  this 
application  is  given  by: 

E2  =  (Er  +  Es)  •  (Er  +  Es) 


E2  =  1  Er  |  2  +  |  Es  1  2  +  2|Ep-Es| 


I  =  I  +  I  +  I 

r  s  rs 


(100) 

(101) . 

Multiplying  by  the  appropriate  constants  and  taking  the  time- 
averaged  value  of  the  above  expression  yields: 

(102)  . 

The  last  term,  I..,.,  is  the  interference  term  which  is  evaluated 
by  substituting  equations  95  and  96  into  equation  99  as  follows: 
Irs  =  2ev<|Er-Ej> 

Irc  =  2  e  v<  |  E  •  E  cos  (fl  z-ut+£)  cos  (/3  z-ut->{|> 

Irs  =  2  e v<  |  Eo •  Ec  [  cos  (0rz+  ~  )  cos.  t  +  sin  (£?rz+i  )  sin_  t  ]  x 
[cos  (/?sz+ *  )  cos-.-.t  +  sin  )  sinwt  J  |  > 
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(102)  . 


Irs  =  2€V<IEo*Eo  {COS2wtCOS(/9rZ  +  OcOS(/?sZ  +  0  + 

sin2wtsin  (0rz+£  )  sin  0?sz+£  )  + 
cos^tsinwt  x  [cos  (Prz+(  )  sin  (0sz+O  + 
sin(/?rz+£  )cos(/?$z+£)  ]  )  |> 

Using  the  fact  that  (29:278): 

<sin2wt>  =  <cos2ut>  =  h 
and  <cos^tsinwt>  =  0, 

the  interference  term  (equation  102)  can  be  rewritten  as: 

1  rs  =  2eV<lEo*Eo  [HcOS(/3rZ+^)cOS(/3sZ  +  C)  + 

^sin(^rz+f )  sin(/?sz+0  ]  I  >  (103). 

Equation  103  can  be  further  simplified  using  the  trigonometry 
identity: 

cos(a)cos(b)  +  sin(a)sin(b)  =  cos(a-b), 
which,  when  used  in  equation  103  yields: 

Irs  =  ev<|Eo-E0  cos(/3rz-/?rz)  |  >  (104). 

Equation  104  is  no  longer  a  function  of  time,  and  the  time- 
averaged  value  is  simply: 

1  rs  =  ^[E02  COS(^rZ-/?rZ)  ] 

Irs^  =  210  cos.5  f  1  0 5  ) 

where  6  =  phase  difference  or  the  difference  in  the  optical  path 
length . 

The  non-interference  terms  in  the  irradiance  expression  are 
correspondingly  evaluated  as  (29:278): 

=  t  =  t  '  Ei- 

Finally,  the  total  irradiance  can  now  be  expressed  as: 

I  =  Ir  +  Ts  +  Irs 
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I  =  I„  +  I„  +  2I„cos<5 
0  0  0 

I  =  2I0(1  +  cos<5 ) 

I  =  4I0cos2(^6)  (106)  . 

The  total  irradiance  function  is  plotted  with  respect  to  the 
phase  difference  (6)  in  Figure  III. 23.  This  Figure  shows  how 
phase  changes  are  converted  into  detectable  amplitude  variations. 

Now  that  the  effect  of  changing  the  phase  difference  has 
been  illustrated,  the  exact  effect  of  strain  on  the  optical 


Figure  III. 23.  Total  Irradiance  Function  Plotted  with 

Respect  to  the  Phase  Difference 
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fiber  is  examined.  In  the  section  on  optical  fibers,  equation  17 
described  optical  fibers  as  a  bulk  piece  of  material  which 
undergoes  a  longitudinal  change  in  length  in  response  to  a 
specifically  applied  force.  To  examine  this  characterization 
more  closely,  the  actual  phase  difference  in  an  optical  fiber 
exposed  to  an  external  force  is  explored. 

The  difference  in  path  lengths  in  the  sensor  and  reference 
arms  is  the  phase  difference  (6),  which  is  given  by: 

<5  =  /3rZ-/?sZ  =  pi  (107)  . 

A  phase  shift  occurs  when  the  phase  difference  changes;  that  is: 

A 6  =  L  (pi)  (108)  . 

When  the  fiber  length  changes,  the  change  in  the  phase  difference 
is  due  to  both  the  length  change  and  an  induced  change  in  the 
index  of  refraction  of  the  glass  fiber.  This  phase  difference 
change  is  given  by  (3:88): 

A6  =  2  m  L  f  A  L  +  An]  (109) 

*  L  nj 

A<5  =  2z_nLi  75  -  hn2  [  (pn  +  P12)  7,  +  P1273]>,  (110) 

where  n  =  index  of  refraction  of  the  core  =  1.461 
L  =  initial  length  of  the  sensor  arm 
>  =  wavelength  of  the  excitation  source 
7,  -  change  in  length  in  the  radial  direction  =  -.>7 
7 3  =  change  in  length  along  the  optical  axis 
pn  =  photoelastic  tensor  element 
p12  =  photoelastic  tensor  element. 

The  constants  pn  and  p,2  are  0.12  and  0.27,  respectively,  for 
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optical  fibers  in  general.  Substituting  the  numerical  values 
into  the  above  expression  yields: 

A<S  =  2n_nL(7j  ~  O.O873)  1  (HI) 

A 

which  shows  that  the  major  contribution  to  phase  difference 
arises  from  a  change  in  the  fiber's  length.  Therefore,  it  is  a 
valid  approximation  to  only  attribute  phase  change  difference  to 
a  change  in  length  along  the  optical  axis  of  the  fiber. 

Other  factors  may  also  affect  a  change  in  the  phase 
difference  for  a  magnetostrict ive  fiber  optic  interferometer. 
Temperature  variations  can  affect  the  fiber  optic  sensor  in  two 
ways:  (i)  linear  expansion  of  the  magnetostrictive  material  and, 
(ii)  index  of  refraction  changes  (30:1478).  Figure  III. 24 
depicts  how  calculated  thermal  variations  may  cancel  the  effects 
of  magnetostriction  in  nickel.  In  this  research  project,  thermal 
variations  were  minimized  by  thermally  insulating  the  optical 
fibers  with  cloth,  and  placing  the  sensor  and  reference  arms  in 
close  proximity  so  that  any  thermal  effects  would  occur  in  both 
arms  to  produce  a  minimal  effect  on  the  sensor's  response. 

The  thickness  of  the  ferromagnetic  material  also  affects  the 
phase  difference  change.  Figure  III. 25  shows  experimental  data 
relating  phase  difference  change  for  different  thicknesses  of  a 
nickel  jacket.  In  this  research  project,  the  magnetostrictive 
coatings  were  designed  to  be  at  least  fifteen  microns  thick  to 
avoid  compromising  the  phase  change  difference. 
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Figure  III.  24.  Thermal  Fluctuations  Producing  a  Phase  Change 

Identical  to  a  Magnetostriction  Phase  Chanae 
(31:3719) 
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Figure  III. 25.  Effect  of  Magnetostrict ive  Coating  Thickness 

on  the  Phase  Change  Differer.ee  (  32:3747) 
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Critical  Performanc e  Parameters 

In  summary,  the  results  of  this  Chapter's  analysis  revealed 
the  critical  parameters  affecting  the  sensitivity  of  a  MZ  fiber 
optic  magnetic  field  sensor.  These  parameters  are  shown  in  Table 
III . 2 . 

The  magnetostrict ive  material  used  in  the  sensor  design 
incorporates  the  most  significant  influence  on  the  sensor's 
performance.  The  material's  magnetostriction  properties  directly 
affect  the  magnitude  of  the  change  in  the  length  of  the  sensor 
arm.  The  amount  of  magnetostriction  depends  on  the  initial 
orientation  of  the  magnetic  domains.  Maximum  magnetostriction 
occurs  when  the  magnetic  domains  are  initially  aligned 
perpendicular  to  the  direction  of  the  magnetic  field  to  be 

Table  III.  2.  Critical  Parameters  Affecting  Sensor  Performance 


-  Magnetostrictive  Material 

-  Initial  Orientation  of  the  Magnetic  Domains 

-  Magnetostrictive  Material  Thickness 

-  Magnetostrictive  Material  Length 

-  Shape  of  the  Magnetostrictive  Material  Coating 

-  Fiber  Optic  Cable's  Cladding  Thickness 

-  Wavelength  of  the  Optical  Excitation  Source 

-  Sensor  Temperature  Stability 
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In  addition  to  processing  the  material,  the  geometry  of  the 
magnetostrictive  coating  affects  the  sensor's  response.  There  is 
experimental  evidence  in  the  literature  that  the  thickness  of  the 
magnetostrictive  material  affects  the  amount  of  magnetostriction. 
The  amount  of  induced  phase  change  increases  with  increasing 
coating  thickness  until  saturating  at  a  coating  thickness 
associated  with  bulk  material.  The  static  physical  length  of  the 
magnetostrictive  material  is  directly  proportional  to  the  phase 
change  difference  as  shown  in  equation  111.  In  a  packaged 
sensor,  the  length  of  the  coating  could  be  maximized  by  winding 
the  sensor  arm  in  a  coil.  For  this  research  project,  the  length 
is  limited  to  10  cm.  While  not  addressed  in  the  literature,  it 

is  postulated  that  the  shape  of  the  magnetostrictive  material  may 
affect  its  magnetostriction . 

The  geometry  of  the  fiber  optic  cable  directly  affects  the 
arc  -r.t  that  the  cable  will  stretch  under  the  force  applied  the 
magnetostrictive  material.  To  maximize  the  change  in  length,  the 
cross-sectional  area  of  the  fiber  optic  cable  should  to 
minimized.  The  cress- sectional  area  can  be  minimized  by  reducing 
t he  fiber  opt i c  cat  1 c ' s  cladding  thickness. 
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wavelength  excitation  source  is  desired.  However,  the  shorter 
the  wavelength,  the  smaller  the  fiber  optic  core  rust  be  f c  r 
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single-mode  operation.  The  smaller  core  leads  to  fiber  optic 
cable  fabrication  problems  and  difficulties  in  coupling  light 
into  the.  fiber. 

Finally,  the  ambient  temperature  environment  surrounding  the 
sensor  will  also  affect  its  performance.  Thermal  fluctuations 
cause  thermal  expansion  and  changes  in  the  index  of  refraction  of 
the  optical  fiber.  These  thermally  induced  changes  in  the  sensor 
arm  may  be  sufficiently  large  so  as  to  mask  any  magnetostrict ive 
induced  changes.  For  this  project,  thermal  fluctuations  were 
minimized  by  covering  the  fibers  with  cloth  and  placing  the 
sensor  and  reference  arms  in  close  proximity. 

Sensor  Arm  Designs 

Based  on  the  preceding  theory,  specific  sensor  arms  were 
designed  to  test  specific  aspects  of  the  fiber  optic  magnetic 
field  sensor.  These  sensor  arms  were  intended  to  evaluate  the 
following  sensor  parameters: 

1)  magnetostrict ive  material 

2)  bulk  material  vs.  thin  films 

3)  initial  orientation  of  magnetic  domains 

4)  magnetostrict ive  material  coating  geometry 

5)  fiber  opt:c  cable  cladding  thickness 

A  series  of  experiments  were  carefully  planned  to  oval uatc  ea  h 
of  these  sensor  parameters. 

E>'Herimcnt _ 1.  The  first  experiment  was  designed  to  evaluate 

the  most  promising  magnetostrict ive  material.  To  aeeompl ish  this 
objective,  three  sensor  arms  were  designed  to  be  10  cr  in  lens*.1: 
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with  a  15  urn  thick  uniform  coating  of  a  particular 
magnetostrictive  material  (cobalt,  iron,  and  nickel). 

Experiment  2.  The  second  experiment  was  designed  to  compare 
the  sensor’s  performance  relative  to  thin  film  coated  fibers  and 
fibers  bonded  to  bulk  material.  To  accomplish  this  objective, 
the  sensor  arms  were  designed  to  be  three  identical  10  cm  long 
bulk  rods  of  the  magnetostrictive  materials  bonded  to  optical 
f ibers . 

Experiment  3.  The  third  experiment  was  designed  to  evaluate 
the  effect  of  initially  orienting  the  magnetic  domains.  To 
accomplish  this  objective,  the  sensor  arms  from  the  first 
experiment  were  planned  to  be  heated  above  their  Curie 
temperature,  and  then  allowed  to  cool  in  the  presence  of  an 
externally  applied  static  magnetic  field  flux  density  of  1000  G. 

Experiment  4.  The  fourth  experiment  investigated  the  effect 
of  coating  geometry  on  the  sensor's  performance.  Using  the  most 
premising  material  from  the  first  experiment,  an  ensemble  of 
sensor  arms  were  designed.  The  length  of  the  magnetostrictive 
material  was  limited  to  10  cm,  and  its  thickness  was  controlled 
to  15  urn.  The  sensor's  response  with  the  uniform, ally  coated 
sensor  arm,  of  the  first  experiment  was  planned  to  be  compared  tc 
the  sensor's  performance  with  fibers  coated  with  the  following 
geometries:  (i)  semi-cy1 4 ndrical  coating,  and  (ii)  spiral 
coatings  with  void  spacings  of  1  mm,  5  mm,  and  10  mm. 

Experiment  5.  This  experiment  was  designed  to  determine  the 
effects  of  fiber  optic  cable  cladding  thickness  on  sensor 
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performance.  The  magnetostrictive  coating  on  each  fiber  in  this 
experiment  was  designed  to  be  carefully  controlled  (10  cm  length 
and  15  um  thickness) .  Prior  to  RF  sputtering  the 
magnetostrictive  material  on  the  fiber,  the  cladding  on  four 
fiber  were  designed  to  be  thinned  to  diameters  of  125,  ICO,  75, 
and  50  um.  These  diameters  correspond  to  normalized  areas  of 
100%,  64%,  36%,  and  16%  with  respect  to  the  125  um  diameter 
fiber. 
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IV.  Experimental  Procedure 


This  Chapter  describes  the  equipment  configuration  and 
experimental  procedures  that  were  used  to  evaluate  the  sensor. 
This  Chapter  describes  the  Helmholtz  coil  calibration  procedure, 
sensor  arm  fabrication,  interferometer  arrangement,  sensor 
sensitivity  measurements,  and  sensor  optimization.  The  results 
of  these  procedures  are  reported  in  the  following  chapter. 
Helmholtz  Coil  Calibration 

In  order  to  verify  that  the  Helmholtz  coil  generated  the 
magnetic  field  flux  density  which  was  predicted  by  its 
theoretical  analysis,  the  Helmholtz  coil  was  calibrated  using  a 
commercial  Hall  effect  magnetometer.  The  magnetometer  available 
at  AFIT  (Walker  Scientific,  Model  MG-7D,  Worcester,  MA)  was  not 
sufficiently  sensitive,  so  a  more  responsive  magnetometer  (F.  W. 
Bell  Inc.,  model  620,  Columbus,  OH)  was  borrowed  from  Tim 
Peterson  of  the  Materials  Laboratory  to  perform  the  calibration. 
After  the  ambient  magnetic  field  noise  was  measured,  the  magnetic 
field  flux  density  produced  by  the  Helmholtz  coil  was 
characterized . 

Since  the  sensor's  sensitivity  goal  was  to  detect  magnetic 
flux  densities  below  the  earth's  ambient  magnetic  field  of 
approximately  0.1  G  (2:360),  a  magnetically  shielded  test  chamber 
was  constructed.  A  set  of  hollow  spheres  fabricated  fr^.n  a 
material  with  high  permeability  would  represent  an  ideal 
magnetically  shielded  chamber  (17:9).  However,  due  to  practical 


IV- 1 


considerations,  a  first  level  magnetically  shielded  chamber  was 
constructed.  The  AFIT  model  shop  built  a  plexiglass  rectangular 
chamber  which  measured  31  cm  x  51  cm  x  31  cm.  This  size 
permitted  the  Helmholtz  coils  to  be  placed  inside  the  chamoer. 

The  enclosure  had  two  2.5  cm  diameter  holes  as  ports  for  the 
optical  fibers  and  a  25  pin  connector  for  electrical  connections 
to  the  Helmholtz  coil.  The  surfaces  of  the  plexiglass  were 
sprayed  with  a  nickel-impregnated  polymer  coating  (Miller- 
Stephenson  Chemical  Co.,  model  MS-485N,  Danbury,  CT) .  However, 
the  conducting  paint  supply  was  exhausted  before  an  even  coat 
could  be  applied  to  the  chamber.  As  a  result,  the  conductivity 
of  the  surfaces  varied;  for  example,  the  resistance  across  the 
chamber's  lid  was  measured  to  be  seven  ohms.  To  further  lower 
the  chamber  surfaces'  resistance,  aluminum  foil  was  placed  on 
both  sides  of  each  surface.  This  foil  was  electrically  connected 
to  the  ground  terminal  on  the  Helmholtz  coil.  The  resistance 
between  any  surface  on  the  enclosure  and  the  Helmholtz  ground 
connection  was  measured  to  be  one  ohm  or  less  on  the  average. 

Before  attempting  to  make  any  measurements,  the  magnetometer 
was  calibrated  for  DC  magnetic  fields.  AC  calibration  was  not 
attempted  since  reference  AC  magnetic  field  sources  were  not 
available.  The  magnetometer  readings  were  taken  when  the  probe 
was  placed  insid^  a  zero  gauss  chamber  and  two  reference  magnets 
of  95  and  978  G,  respectively.  Reference  magnets  with  smaller 
magnetic  flux  densities  would  have  been  more  appropriate,  but 
were  not  available. 
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To  determine  the  ambient  magnetic  field  flux  density,  the 
magnetometer  probe  was  placed  in  various  positions  and 
orientations,  both  inside  the  chamber  and  around  the  laboratory. 
The  magnetometer's  readings  were  then  recorded  with  no  current 
flowing  in  the  Helmholtz  coil. 

After  the  ambient  magnetic  field  flux  density  was 
determined,  power  was  applied  to  the  Helmholtz  coil.  For  the  DC 
measurements,  the  current  was  supplied  with  a  power  supply 
(Hewlett-Packard,  model  HP  6236B,  Palo  Alto,  CA)  and  for  the  AC 
measurements,  the  current  was  provided  with  a  function  generator 
(Wavetek,  model  23,  San  Diego,  CA)  and  a  power  amplifier 
(Hewlett-Packard,  model  HP  467A,  Palo  Alto,  CA) .  The  power 
amplifier  was  used  when  currents  between  120  and  500  mA  were 
reguired.  (Larry  Reitz  of  the  Electronic  Technology  Laboratory 
provided  the  power  amplifier) .  The  freguency  of  the  AC  signal 
was  chosen  to  be  25  Hz  since  this  freguency  is  in  the  0  to  30  Hz 
range  of  interest,  and  was  not  an  AC  power  harmonic.  During 
measurements,  the  current  was  monitored  with  a  multimeter  (Fluke, 
model  8010A ,  Everett,  WA) . 

Once  power  was  applied  to  the  Helmholtz  coil,  DC  and  AC 
measurements  were  made.  The  first  measurement  was  accomplished  to 
determine  the  magnetic  field  flux  density  at  the  center  point 
between  the  two  coils  for  different  amounts  of  applied  current. 
For  the  second  set  of  measurements,  a  constant  current  of  250  mA 
was  applied  to  the  Helmholtz  coil.  The  magnetic  field  flux 
density  at  the  center  point  between  the  two  coils  was  then 
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determined  for  different  numbers  of  turns  in  the  coils.  The 
final  measurement  was  implemented  to  evaluate  the  magnetic  field 
flux  density  at  various  spatial  locations  between  the  two  coils 
for  a  constant  applied  current  and  fixed  number  of  turns.  The 
spatial  variations  of  the  magnetic  field  flux  density  was 
determined  by  placing  the  2  mm  x  5  mm  x  0 . 5  mm  probe  tip  in 
various  positions  on  and  off  the  axis  between  the  coils. 

Sensor  Arm  Fabrication 

After  the  model  biomagnetic  field  source  was  generated  and 
characterized,  the  sensor  arms  were  fabricated.  The  two  classes 
of  sensor  arms  used  in  this  research  project  included:  (i) 

fibers  that  were  coated  with  thin  films  of  magnetostrictive 
material,  and  (ii)  fibers  that  were  attached  to  bulk 
magnetostrictive  material.  The  preparation  of  the  fibers, 
application  of  the  magnetostrictive  material  to  the  fiber,  and 
sensor  arm  characterization  are  described  in  this  section.  The 
sensor  arms  fabricated  were  based  on  the  designs  presented  in  the 
previous  chapter.  However,  there  are  some  differences  between 
the  sensor  arms  designed  and  those  fabricated.  The  difference- 
resulted  from  fabrication  difficulties  and  actual  sensor 
performance  results.  A  table  summarizing  the  characteristics  of 
the  fabricated  sensor  arms  is  included  in  the  next  chapter. 

Preparation .  Both  classes  of  sensor  arms  required  that  the 
optical  fiber  be  prepared.  One  meter  was  chosen  as  the  length  ol' 
the  sensor  arm  since  this  distance  was  adequate  to  make 
connections  in  the  interferometer  without  having  excess  fiber 
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length.  The  one  meter  length  was  also  sufficient  for  the  mode¬ 
stripping,  jacketing  material  to  absorb  extraneous  modes 
propagating  in  the  cladding.  The  extraneous  modes  are  introduced 
into  the  fiber  by  non-ideal  coupling  of  light  into  the  fiber. 
Surplus  fiber  length  in  addition  to  the  one  meter  standard  would 
just  add  environmental  noise  to  the  sensor. 

In  order  for  the  magnetostrictive  material  to  have  a  maximum 
effect  on  the  fiber,  the  magnetostrictive  material  must  be  as 
close  to  the  fiber's  core  as  possible.  The  first  step  in 
achieving  this  proximity  was  to  remove  the  jacketing  material 
around  a  10  cm  interior  segment  of  the  fiber.  This  length  was 
chosen  to  match  the  extent  of  the  Helmholtz  coil's  uniform 
magnetic  field.  The  entire  length  of  fiber  was  not  stripped 
since  the  jacketing  material  provides  mechanical  strength  and 
absorbs  extraneous  modes.  Without  the  jacket,  the  fiber  is 
extremely  fragile.  Since  an  interior  section  of  fiber  was  to  be 
stripped,  a  microstripper  was  not  appropriate.  Instead,  the 
section  of  fiber  to  be  stripped  was  placed  in  a  chemical  solvent 
(Star  Bronze,  Model  Zipstrip,  Alliance,  OH)  for  a  few  minutes. 

The  chemical  solvent  softened  the  jacketing  material  enough  such 
that  this  layer  peeled  away  after  being  gently  rubbed.  Although 
this  method  of  stripping  an  interior  length  produced  the  best 
results,  the  method  was  not  ideal  since  approximately  one  quarter 
of  all  the  fibers  broke  while  being  stripped. 

Once  the  fiber  was  stripped  of  its  jacketing  material,  the 
bare  silica  was  cleaned.  The  cleaning  process  involved  a  multi- 
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step  procedure  consisting  of  sequentially  dipping  the  fiber  in 
acetone,  methanol,  buffered  hydrofluoric  acid,  and  deionized 
water.  The  fiber  was  allowed  to  remain  in  each  solution  for  two 
minutes.  The  acetone  and  methanol  served  to  remove  the  chemical 
solvent  and  an^  remnants  of  the  jacketing  material.  The  4:1 
NH4F:HF  buffered  hydrofluoric  acid  solution  slightly  etched  the 
bare  silica,  which  was  done  to  enhance  adhesion  of  the 
magnetostrictive  material  to  the  fiber.  Following  cleaning,  the 
fiber  was  placed  in  an  oven  at  70°C  for  twenty  minutes  to  remove 
absorbed  moisture.  The  caustic  nature  of  the  cleaning  process 
required  gloves  to  be  worn  while  cleaning  the  fibers.  The 
handling  required  to  implement  the  cleaning  resulted  in 
approximately  one  third  of  the  fibers  breaking. 

Thin  film  sensor  arms.  In  the  first  class  of  sensor  arms, 
the  thin  film  of  magnetostrictive  material  was  applied  using  a  RF 
sputtering  system  (Denton  Vacuum,  Inc.,  Model  DV  602,  Cherry 
Hill,  NJ) .  This  following  discussion  describes  the  procedures 
used  to  sputter  the  metal,  thermally  anneal  the  fibers,  and 
characterize  them. 

To  fit  the  fiber  in  the  sputtering  chamber,  the  jacketed  end 
lengths  of  the  fiber  were  coiled  into  approximately  four 
centimeter  circles,  which  were  loosely  held  together  with  masking 
tape.  The  coiled  ends  of  the  fiber  were  then  placed  on  two 
stacks  of  four  microscope  slides  such  that  the  bare  section  of 
fiber  was  suspended  approximately  0.5  cm  above  the  deposition 
system's  vacuum  chamber  floor.  The  coiled  ends  were  then  masked 
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from  the  deposited  metal  by  placing  microscope  slides  over  them. 

Masking  was  also  used  on  the  bare  interior  section  of  one 
fiber  to  produce  a  geometrical  design  with  the  deposited  metal. 

In  order  to  produce  a  spiral  coating  configuration,  a  0.05  mm 
nickel  wire  was  wrapped  around  the  fiber  with  a  pitch  of  2  mm. 
After  the  metal  was  deposited,  the  nickel  wire  was  carefully 
removed.  In  order  to  produce  an  approximate  semi-cylindrical 
pattern,  the  fiber  was  positioned  on  a  piece  of  masking  tape. 

Once  the  fibers  were  placed  in  the  chamber  directly  beneath 
the  sputtering  target,  the  enclosure  was  evacuated,  and  the  RF 
power  was  applied.  The  three  metal  sputtering  targets  (cobalt, 
iron,  and  nickel)  and  fresh  fiber  samples  were  interchanged 
between  sputtering  runs.  After  the  chamber  was  pumped  down  to  2 
x  10"5  Torr,  increasing  amounts  of  argon  gas  were  introduced  into 
the  chamber.  With  the  RF  forward  power  set  at  500  W,  the  argon 
pressure  was  allowed  to  increase  until  the  plasma  formed.  Plasma 
formation  was  marked  by  a  blue  glow  in  ^he  chamber  and  a  drastic 
reduction  in  the  RF  reflected  power.  The  argon  pressure  required 
ranged  from  7  to  15  Torr.  Once  the  plasma  was  formed,  the  argon 
pressure  was  decreased  to  just  sustain  the  plasma.  For  nickel, 
this  pressure  was  0.5  Torr,  but  for  cobalt  and  iron,  this 
pressure  was  6  to  7  Torr.  The  argon  >  ressure  was  decreased  to 
prevent  the  cryo  pump  from  becoming  saturated  before  the 
sputtering  run  was  finished.  After  the  sputtering  process  had 
continued  for  the  time  required  for  the  desired  film  thickness, 
the  RF  power  was  shut-off,  and  the  fibers  were  allowed  to  cool 
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for  one  hour.  Once  the  fibers  cooled,  they  were  removed  and 
immediately  attached  to  plexiglass  strips  for  mechanical  support. 

In  order  to  test  the  effects  of  thermal  annealing  on  the 
sensor's  performance,  three  of  the  sputtered  fibers  were  later 
placed  in  an  oven  (Sheldon  Manufacturing,  Model  1610D,  Cornelius, 
OR) .  The  fibers  were  arranged  on  a  glass  plate  and  heated  to  the 
maximum  oven  temperature  of  350°  C  for  two  hours.  Although  this 
temperature  is  below  the  Curie  temperature  of  the 
magnetostrictive  materials  investigated,  it  was  the  maximum 
temperature  available  in  an  oven.  In  order  to  test  the  effects 
of  thermal  annealing  in  the  presence  of  an  external  magnetic 
field,  the  fibers  were  again  placed  on  the  glass  plate  which  was 
resting  on  a  laboratory  horseshoe  magnet.  The  magnetic  field 
flux  density  of  the  horseshoe  magnet  was  measured  to  be  1000  G  at 
the  center  line  between  its  two  poles.  The  fibers,  glass  plate, 
and  magnet  were  all  placed  in  the  oven  and  heated  to  350°  C  for 
two  hours.  After  this  thermal  annealing  c'7  :le,  the  power  to  the 
oven  was  shut-off,  and  the  oven  was  allowed  to  cool  overnight 
before  the  fibers  were  removed. 

The  metalizea  fibers  were  characterized  for  film  thickness 
and  quality.  The  film  thickness  was  determined  using  a  glass 
slide  that  was  in  the  sputtering  chamber  with  the  fibers.  This 
glass  slide  was  partially  masked  with  a  grating.  A  profilometer 
(Sloan  Dektak,  Model  900051,  Santa  Barbara,  CA)  was  used  to 
determine  the  metal  step  height  between  the  masked  and  unmasked 
portions  on  the  slide.  To  determine  the  film  quality,  an  optical 
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microscope  (Heerbrugg,  Model  Wild,  Switzerland)  was  used.  The 
microscope  was  used  to  examine  the  surface  of  the  fiber  for  film 
roughness.  In  order  to  determine  if  the  metal  was  uniformly 
deposited  on  the  top  and  bottom  sides  of  the  fiber,  a  metalized 
fiber  was  cleaved  and  then  inspected  with  the  optic  axis  of  the 
fiber  facing  the  microscope. 

Bulk  sensor  arms.  The  second  class  of  sensor  arms  was 
composed  of  fibers  bonded  to  a  bulk  piece  of  magnetostrict ive 
material.  A  30  cm  long,  2.5  cm  diameter  nickel  rod  was  obtained 
from  Art  Beecraft  of  the  Electronic  Technology  Laboratory.  This 
rod  was  labeled  as  2  0  1%  nickel.  Other  samples  of  high  purity 
magnetostrict ive  materials  could  not  be  obtained  before  the 
research  project  was  completed.  In  order  to  perform  an  accurate 
comparison  of  the  different  magnetostrictive  materials,  it  would 
be  necessary  to  evaluate  identically  sized  samples  of  the 
different  materials. 

Once  the  fiber  was  cleaned  and  prepared,  it  was  bonded  along 
the  length  of  the  nickel  rod  with  cyanoacrylate  adhesive 
(Loctite,  Model  MIL-A-46050C,  Puerto  Rico).  The  procedure  used 
to  adhere  the  fiber  to  the  rod  was  accomplished  by  holding  the 
fiber  flat  against  the  rod  and  gluing  one  end  of  the  fiber  to  rod 
and  allowing  it  to  dry.  The  fiber  was  then  stretche  1  taut  and 
the  other  end  was  bonded  and  allowed  to  dry.  After  both  ends 
were  secure,  the  entire  length  of  the  fiber  was  glued  to  the  rod. 
Variables  in  this  process  include  the  amount  of  glue  used  and  the 
degree  of  tension  on  the  fiber  as  it  was  being  attached  to  the 
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rod.  After  the  rod  and  fiber  combination  was  evaluated,  the 
fiber  and  glue  was  scrapped  off,  and  the  rod  was  rubbed  with 
acetone . 

With  only  one  magnetostrict ive  bulk  material  sample,  the 
various  sensor  arms  were  fabricated  seguentially .  The  first 
three  sensor  arms  were  used  to  test  the  effect  of  thermal 
annealing  on  sensor  performance.  The  first  bulk  sensor  arm  was 
simply  a  single  length  of  fiber  bonded  to  the  nickel  rod.  For 
the  second  sensor  arm,  the  nickel  rod  was  placed  in  an 
oxidation/diffusion  furnace  (Thermco,  Model  201,  0 range,  CA)  at 
758°  C  for  three  hours.  This  temperature  was  chosen  because  it 
was  above  the  Curie  temperature  of  nickel.  After  allowing  the 
rod  to  cool  overnight,  a  fiber  was  attached  to  the  rod.  The 
third  bulk  sensor  arm  was  processed  identically  to  the  second 
except  that  the  nickel  rod  cooled  while  resting  on  the  glass 
plate  and  under  the  influence  of  the  1000  G  magnet. 

The  final  two  bull:  sensor  arm  configurations  involved  no 
further  processing  of  the  nickel  rod,  but  did  require 
modifications  of  the  fiber  attached  to  it.  The  fourth  bulk 
sensor  arm  incorporated  a  thinned  fiber.  The  thinned  fiber  was 
prepared  by  etching  it  in  a  1:1  HF:NH4F  solution  for  forty 
minutes.  The  etched  fiber  was  compared  to  a  full-sized  fiber 
with  an  optical  microscope.  The  thinned  fiber  was  then  attache! 
to  the  rod.  Thinned  fibers  were  riot  RF  sputtered  because  they 
would  not  have  survived  the  handling.  The  last  bulk  sensor  arm 
was  fabricated  to  test  the  effect  of  the  active  region  of  the 
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fiber's  length  on  the  sensor's  performance.  To  accomplish  this 
objective,  a  single  fiber  with  three  10  cm  bare  silica  lengths 
was  attached  to  the  nickel  rod. 

Interferometer  Arrangement 

After  the  sensor  arms  were  fabricated,  the  renainder  of  the 
sensor  was  assembled.  This  section,  which  describes  the  sensor's 
configuration,  is  subdivided  into  two  parts.  The  first  part, 
fiber  optic  cable  fundamentals,  describes  the  eguipment  and 
procedures  used  to  work  with  the  fiber  optic  cables.  The  second 
part  details  the  remaining  components  and  procedures  used  to 
fabricate  three  versions  of  a  fiber  optic  Mach-Zehnder 
interferometer . 

Fiber  optic  cable  fundamentals.  Working  with  fiber  optic 
cables  reguired  developing  technigues  to  cleave  the  enas  of 
fibers,  couple  light  into  the  end  of  the  fiber,  and  make  fiber 
splices . 

In  order  to  minimize  loss  at  each  fiber  connection,  the  ends 
of  the  fiber  must  be  smooth,  flat,  and  perpendicular  to  the 
optical  axis  of  the  fiber.  Properly  prepared  fiber  ends  were 
achieved  by  cleaving  the  fiber.  The  first  step  in  cleaving  a 
fiber  was  to  use  a  microstripper  (Utica,  Model  MS,  Orangeburg, 

SC)  to  remove  the  fiber  optic  cable's  jacket  material  along  a 
length  that  is  within  4  cm  of  the  fiber's  end.  Microstrippers 
operate  just  like  wire  strippers.  The  fiber  was  then  cleaved 
with  a  fiber  cleaver  (Fujikura,  Model  CT-02,  New  York,  NY),  which 
is  a  small  device  that  scores  the  fiber  with  a  carbide  edge,  and 
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then  bends  it  to  cause  a  clean  break. 


The  fiber  end  was  then  inspected  using  the  fusion  splicer 
(Orionics,  Model  FW-303,  Bozeman,  MT) .  The  fusion  splicer 
possesses  the  necessary  optics,  lamp,  and  viewer  to  accomplish 
this  task.  If  the  fiber's  end  surface  was  not  as  smooth  when 
compared  to  the  example  pictures  included  in  the  fusion  splicer 
manual,  a  small  amount  of  heat  was  applied  to  the  end  of  the 
fiber  to  flatten  its  surface.  Both  the  duration  and  amount  of 
heat  applied  are  contronabie  variables  associated  with  the 
fusion  splicer.  Another  method  for  examining  the  quality  of  a 
cleaved  fiber  is  to  examine  the  quality  of  laser  light  exiting  a 
fiber.  A  high  quality  cleave  produces  a  uniform  and  sharply- 
defined  circular  spot,  whereas  a  low  quality  cleave  produces  a 
distorted  and  irregularly-shaped  region  of  light. 

The  last  method  o?  testing  a  fiber's  cleave  requires 
launching  light  into  the  fiber.  Coupling  laser  light  into  the 
end  of  a  fiber  involves  aligning  a  laser  beam  whose  diameter  is 
approximately  1  mm  and  focusing  it  down  to  the  fiber's  core 
diameter  of  3.5  urn.  This  task  was  achieved  using  a  fiber  light 
coupler  (Newport,  Model  F-916,  Fountain  Valley,  CA) ,  which  was 
borrowed  from  Dale  Stevens  of  the  Electronic  Technology 
Laboratory.  The  fiber  light  coupler  consists  of  a  20X  microscope 
objective  and  a  brass  ferrule  fiber  holder.  The  fiber  holder  has 
micrometer  adjustments  in  three  perpendicular  directions  and  tilt 
adjustments  in  two  directions. 

To  couple  light  into  the  fiber,  it  was  important  that  the 
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beam  be  perpendicular  tc  the  plane  of  the  microscope  objective. 

A  slight  departure  from  the  perpendicular  condition  prevents 
efficient  coupling  of  the  light  into  the  fiber.  Before  the  fiber 
holder  was  positioned  in  the  light  coupler,  the  beam  was  examined 
to  ensure  that  it  was  passing  through  the  center  of  the  aperture 
where  the  fiber  holder  is  secured.  Since  rhe  light  is  rapidly 
diverging  after  being  focused,  a  clean  circular  spot  indicated 
proper  alignment.  One  end  of  the  fiber  was  then  positioned  in 
the  fiber  holder  with  approximately  2  mm  of  the  fiber  protruding 
from  the  holder,  while  the  other  end  was  attached  to  the  optical 
power  meter  (described  later) .  The  2  mm  length  of  the  extended 
fiber  was  essential  for  positioning  the  fiber's  end  the  correct 
distance  from  the  microscope  objective. 

The  technique  used  to  couple  light  into  the  fiber's  core  was 
similar  to  the  method  used  to  couple  light  ehrough  a  spatial 
filter.  That  is,  the  fiber  holder  was  initially  positioned  a 
great  distance  from  the  microscope  objective  such  that  the 
diverging  beam's  size  was  very  large,  which  also  ensures  that 
some  of  the  light  enters  the  fiber  core.  The  other  two 
positioning  micrometers  were  adjusted  to  maximize  the  light  on 
the  output  power  detector  which  served  to  center  the  core  on  the 
beam  since  the  TEM0  0  beam  had  more  power  towards  the  center  of 
the  beam.  The  two  tilt  micrometers  were  also  adjusted  to 
maximize  the  output  power.  After  the  two  tilt  micrometers  were 
adjusted  once,  it  was  not  necessary  to  change  them.  The  fiber 
holder  was  then  moved  several  full  rotations,  and  the  two 
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positioning  micrometers  were  used  to  maximize  the  output  power. 
This  procedure  of  rotating  the  fiber  holder  and  adjusting  for 
maximum  output  power  was  repeated  approximately  six  times  until 
the  output  power  was  maximized. 

Most  applications  of  fiber  optic  cables  require  that  two 
fibers  be  joined  together  in  a  splice.  Two  pieces  of  equipment 
were  investigated  to  make  fiber  splices:  the  fusion  splicer  and 
glass  ferrule  tubes  (Norland,  Model  20125,  New  Brunswick,  NJ) . 

The  best  method  for  formina  a  splice,  regardless  of  which 
splicing  equipment  was  used,  was  to  begin  with  laser  light  being 
launched  into  the  first  of  the  two  fibers  that  are  being  spliced. 
The  light  from  the  first  fiber  was  measured  with  the  optical 
power  meter  before  being  joined  to  ‘he  second  fiber.  The  light 
from  the  second  fiber  was  then  monitored  with  the  optical  power 
meter  as  the  two  fibers  were  being  adjusted  and  spliced. 

The  fusion  splicer  joins  two  fibers  together  by  simply 
melting  their  ends  together.  The  glass  ferrule  consists  of  a 
precision  made  tube  whose  internal  diameter  ports  are  only 
slightly  greater  than  the  diameter  of  the  fiber.  To  form  a 
splice  with  the  glass  ferrules,  each  fiber  end  was  dipped  in 
index  matching  gel  and  inserted  into  the  glass  ferrule  until  the 
two  fibers  met  in  th-  n ter.  The  glass  ferrule  was  then  mounted 

on  a  splice  holder  (Norland,  Model  Splicemate,  New  Brunswick,  NJ ) 
to  prevent  unwanted  movement  of  the  glass  ferrule  or  fibers. 

Fiber  optic  Mach-Zehnder  interferometers.  Once  the 
techniques  required  to  work  with  fibers  were  established,  three 
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different  fiber  optic  Mach-Zehnder  interferometers  were 
fabricated  to  investigate  particular  aspects  of  the  sensor.  The 
three  interferometers  are  illustrated  in  Figures  IV. 1,  IV. 2  and 
IV. 3.  The  GaAs  fused  coupler  interferometer  was  only  constructed 
to  verify  the  performance  of  the  fused  couplers.  This  part 
describes  the  remaining  unique  equipment  used,  criteria  for 
choosing  the  equipment,  and  the  arrangement  of  the 
interferometers.  Standard  optical  equipment  such  as  lenses, 
mirrors,  and  beam  splitters  will  not  be  described  in  detail. 

The  first  requirement  for  assembling  an  interferometer  was 
to  have  an  optical  bench  that  was  isolated  from  floor  vibrations. 
The  3.5  m  x  1.2  m  air-cushioned  optical  bench  in  Room  127, 
Building  640,  provided  adequate  vibration  isolation.  When  the 
air  table  was  floating  correctly,  the  interferometer  produced 
useful  results,  such  as  a  visible  fringe  pattern.  When  the  air 
table  rested  on  its  support  platform  due  to  loss  of  the 
compressed  air  utility,  no  fringe  pattern  was  visible.  Even 
partial  floating  of  the  table  caused  problems.  The  fringe 
pattern  was  degraded  significantly  when  one  of  the  compressed  air 
valves  malfunctioned  and  caused  one  leg  to  loose  its  support 
capacity . 

The  optical  components  were  arranged  on  five  one  meter  long 
rectangular  rails  (Oriel,  Stratford,  CT)  that  were  bolted  onto 
the  air  table.  Each  component  was  mounted  on  the  rail  using 
translation  stages  (Oriel,  Models  11621,  14781,  and  16021, 
Stratford,  CT) .  These  stages  provided  translation  in  three 
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Figure  IV. 1.  HeNe  Fused  Coupler  Interferometer 
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Figure  IV. 2.  GaAs  Fused  Coupler  Interferometer 
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perpendicular  directions.  The  laser  was  mounted  on  a  special 
stage  (Oriel,  Model  16162,  Stratford,  CT)  that  was  specifically 
designed  to  hold  a  laser  tube  and  provided  translational  and  tilr 
adjustments.  The  rails  and  stages  were  necessary  to  achieve  the 
precise  optical  alignment  required  to  couple  light  into  the  fiber 
and  form  interference  patterns. 

The  optical  source  used  for  the  magnetic  field  fiber  optic 
sensor  was  a  2  mW  HeNe  laser  (Oriel,  Model  79360,  Stratford,  CT) . 
The  HeNe  laser  was  selected  because  of  its  availability  and 
performance  characteristics.  The  wavelength  of  the  laser's 
output  was  0.633  urn,  which  is  compatible  with  single-mode 
operation  of  the  York  Technology,  Model  SM600,  fiber  optic 
cables.  As  explained  in  the  interferometer's  operation  section 
of  Chapter  III,  the  laser's  output  power  must  also  be  stable, 
polarized,  and  coherent.  The  laser's  output  power  and 
polarization  were  measured  with  a  radiometer  and  polarizer.  By 
examining  the  fringe  pattern  of  the  interferometer  with  various 
dif f ci er.cc ~  in  the  ler.jth  of  the  two  arms,  -t-he  coherence  length 
of  the  laser  was  examined. 

Two  different  light  detectors  were  used  with  the  different 
interferometer  arrangements.  The  first  detector-  ™?>s  the  optical 
power  meter,  a  device  made  specifically  to  be  used  with  optical 
fibers.  The  components  of  the  optical  power  meter  include:  an 
optical  adapter  (Hewlett-Packard ,  Model  HP  81510C,  Palo  Alto,  C A) 
which  serves  to  hold  the  optical  fiber,  an  optical  head  (Hewlett- 
Packard,  Model  HP  8 1 5 1 1 A ,  Palo  Alto,  CA)  which  houses  the  optics 
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and  detector,  and  the  optical  pulse  power  meter  (Hewlett-Packard 
Model  HP  8 1 5 1 A ,  Palo  Alto,  CA)  which  incorporates  the  detector': 
signal  display  electronics.  With  this  combination  of 
accessories,  the  optical  power  meter  had  the  necessary 
characteristics  to  detect  light  in  the  fiber  optic 
interferometer.  The  detector's  wavelength  sensitivity  spanned 
550  to  950  nm,  which  included  the  laser's  wavelength  (633  nr) . 
The  optical  power  meter's  bandwidth  of  250  MHz  permitted 
detection  of  the  time  varying  light  from  the  interferometer. 
Light  in  the  pW  range  could  be  detected  with  the  optical  power 
meter.  Finally,  the  optical  power  meter  hac  an  optical -to- 
electrical  transducer  that  provided  the  electrical  signal  used  t 
characterize  the  interferometer's  output  light. 

A  radiometer  ( EG&C- ,  Model  450,  Salem,  MA)  was  the  second 
type  of  light  detector  used.  The  radiometer  was  used  when  bulk 
optics  were  used  to  recombine  the  light  beams  in  the 
interferometer.  The  optical  adapter  prevented  using  the  optical 
power  meter  for  this  application.  Like  the  optical  power  meter, 
the  radiometer's  characteristics  were  appropriate  for  the 
interferometer.  The  radiometer  was  sensitive  to  the  0.633  urn 
laser  radiation,  and  its  silicon  detector's  bandwidth  of  6.25  MH; 
was  sufficient  to  follow  changes  in  the  interferometer's  light 
output.  Although  not  as  sensitive  as  the  optical  power  meter, 
the  nW  sensitivity  of  the  radiometer  was  adeguate  to  make 
measurements.  The  radiometer  also  had  an  opt ica 1 -to-electr  ica ! 
transducer . 
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The  final  unique  component  used  in  the  fiber  optic 
interferon*"^  „r  was  the  fiber  optic  fused  coupler  (York 


Technol  dy,  Model  SC-633-50-0.5,  Hampshire,  England  and  Allied, 
Model  945-120-1001,  Spring,  TX)  .  The  fused  coupler,  as  shown  ir. 
Figure  IV.  4,  is  a  four  port  device  that  was  used  to  recombine  tbr 
light  in  two  of  the  interferometers.  It  is  composed  of  dual 
single-mode  fibers  whose  cladding  layers  has  been  thinned  to  loss 
than  a  micron  so  that  the  two  cores  are  in  close  proximity  for  a 
few  centimeters  (33:477).  Through  evanescent  field  coupling, 
light  from  one  core  can  leak  through  the  thinned  cladding  into 
the  next  core.  The  insertion  loss,  w:hich  assumes  light  enters 
only  from  port  1,  is  defined  as  10  log10  [  (P2+P3) /P1  ]  was  specified 
as  -0.34  dB.  The  coupling  ratio,  10  log10  [  P3/ ( P2+P3)  ]  ,  was 
reported  as  -2.91  dB  (34). 
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Once  all  the  major  components  were  chosen,  the 
interferometers  were  arranged  as  shown  in  Figures  IV. 1,  IV. 2,  an 
IV. 3.  Techniques  required  for  all  the  interferometers  are 
described,  and  then  the  special  methods  required  for  each 
particular  interferometer  configuration  are  reported. 

For  all  the  interferometers,  it  was  critical  that  the  laser 
beam  propagate  at  the  same  height  along  or  perpendicular 
to  the  rails.  The  beam's  alignment  was  checked  using  two 
identical  apertures  after  it  passed  through  each  optical  element- 
The  two  apertures  initially  were  positioned  close  to  the  optical 
element.  The  alignment  of  the  optical  element  was  altered  to 
keep  the  beam  passing  cleanly  through  both  apertures  when  the 
second  aperture  was  moved  a  distance  of  one  meter  from  the 
optical  element. 

It  was  also  critical  that  the  fiber  optic  cable,  in  all  the 
interferometer  conf igurat ions ,  to  be  securely  fastened  to  a  flat 
surface.  Any  fiber  suspended  in  the  air  would  vibrate  and  cause 
additional  noise  in  the  sensor.  All  possible  lengths  of  fiber 
were  taped  down  to  a  flat  surface.  The  sections  that  could  nor 
be  in  contact  with  a  surface,  such  as  the  lengths  protruding  t ro 
the  fiber  light  coupler,  were  helc  taut  to  minimize  movement, 
order  to  minimize  the  ambient  air  flow  in  the  laboratory,  the 
ventilation  vents  in  the  ceiling  were  obstructed.  The  fiber 
optic  cable  portion  of  the  sensor  was  also  covered  with  cloth  to 
minimize  any  remaining  ambient  air  turbulence. 

The  HeNe  fused  coupler  interferometer  did  not  require  my 
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arrangement  replicates 


further  techniques  to  be  assembler, 
most  fiber  optic  interferometer  sensors  described  in  the 
literature  with  one  exception.  The  mirror  and  beam  splitter  were 
used  instead  of  another  fused  coupler  with  two  additional 
splices.  The  rational  for  this  modification  was  the  need  for 
equal  power  in  each  arm  of  the  sensor.  While  the  fused  coupler 
would  evenly  split  the  light,  there  would  be  uncertainty  in  the 
loss  of  each  splice.  The  loss  of  the  first  splice  could  be 
characterized,  but  the  instant  the  second  splice  was  connected, 
the  measurement  was  affected  by  the  operation  of  the 
interferometer.  By  using  the  beam  splitter  and  mirror,  the  light 
in  each  arm  was  easily  temporarily  blocked  to  measure  the  optical 
power  in  the  other  arm. 

The  GaAs  fused  coupler  interferometer  required  developing 
techniques  for  working  with  infrared  radiation  from  a  multi-mode 
fiber.  The  optical  source  for  this  interferometer  was  a  GaAs 
laser  diode  which  was  attached  to  a  multi-mode  fiber  pigtail 
(Spectra  Diode  Labs,  Model  SDL-2410-H1,  San  Jose,  GA) .  The  bean 
was  viewed  with  night  vision  goggles  (ITT,  Model  4934 -A,  Roanoke, 
VA) .  The  rapidly  diverging  light  from  the  fiber  pigtail  could 
not  be  directly  coupled  into  the  single-mode  fiber  with  the  fiber 
light  couplers.  To  collimate  the  radiation,  a  graded  index 
(GRID)  lens  (Nippon  Sheet  Glass,  Model  SRL,  Somerset,  NJ)  was 
used.  The  profile  of  the  index  of  refraction  of  a  GRIN  len^  is 
designed  to  collimate  an  incoming  diverging  beam.  Once  the 
infrared  radiation  was  coll  imat.ed ,  no  further  special  techniques 
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were  used  to  fabricate  the  GaAs  fused  coupler  interferometer. 

Unlike  the  HeNe  fused  coupler  interferometer,  the  HeNe 
fringe  pattern  interferometer  required  more  techniques  to  be 
developed  since  it  required  the  light  in  the  sensor  and  reference- 
arms  to  be  collimated  and  recombined.  With  respect  to  the  laser 
and  the  sensor  and  reference  arms,  the  two  interferometers  are 
identical.  For  the  fringe  pattern  interferometer,  the  fibers 
were  terminated  in  fiber  holders  which  were  similar  to  the  fiber 
light  couplers,  except  that  they  did  not  have  the  microscope 
objective.  Since  the  fiber's  core  approximated  a  point  source  of 
light,  the  beams  were  collimated  by  placing  a  25  mm  focal  length 
lens  one  focal  length  from  the  end  of  the  fiber.  The  exact 
position  of  the  lens  was  adjusted  until  the  spot  size  of  the 
laser  beam  was  the  same  just  after  the  lens  and  several  meters 
from  the  lens. 

A  fringe  pattern  was  formed  on  the  radiometer  by  recombining 
the  beams  using  a  mirror  and  beam  splitter.  It  was  necessary  to 
use  a  mirror  mount  with  precision  micrometer  tilt  adjustments 
(Newport,  Model  600A--2,  Fountain  Valley,  CA)  to  form  a  fringe 
pattern.  The  beam  splitter  could  be  mounted  in  a  standard 
laboratory  mount  (Newport,  GM-2,  Fountain  Valley,  CA)  since  once 
it  was  aligned,  it  was  not  altered.  The  fine  adjustments  were 
achieved  using  the  mirror  mount.  No  fringe  pattern  was  visible 
while  the  mirror  was  being  adjusted.  A  small  adjustment  was  made 
and  then  the  interferometer  was  allowed  to  stabilize.  The 
spacing  and  curvature  of  the  fringe  lines  were  determined  by  the 
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alignment  and  collimation  of  the  beams.  Reflections  from  the 
front  and  back  side  of  the  beam  splitter  caused  a  secondary 
fringe  pattern  to  be  formed.  This  secondary  fringe  pattern 
resulted  in  smaller  light  and  dark  bands  within  the  larger  light 
and  dark  bands. 

Once  the  beams  were  recombined,  a  polarizer  and  10X 
microscope  objective  were  used  to  view  the  fringe  pattern.  The 
polarizer  was  necessary  to  ensure  that  the  polarization  of  the 
two  recombined  beams  was  identical.  Although  the  light  from  the 
single  mode  fiber  was  measured  to  be  more  than  ninety  percent 
polarized  in  one  direction,  the  polarization  direction  was 
dependent  on  the  orientation  of  the  fiber.  For  the  highest 
fringe  visibility,  the  output  power  from  each  arm  after  the 
polarizer  was  made  eguivalent.  The  optical  power  was  decreased 
in  one  arm  until  it  matched  the  optical  power  in  the  other  by 
defocusing  one  of  the  fiber  light  couplers  which  caused  less 
light  to  be  launched  into  the  fiber.  The  10X  objective  was  used 
to  magnify  a  portion  of  the  fringe  pattern.  A  0.5  mm  x  0.5  mm 
slit  was  placed  before  the  radiometer  to  sample  a  small  portion 
of  the  fringe  pattern. 

Sensitivity  Measurements 

Once  the  sensor  arms  were  fabricated  ana  the  interferometers 
were  assembled,  the  sensor's  sensitivity  was  measured.  Signal 
processing  equipment  and  techniques  were  required  to  quantify  the 
sensor's  response.  This  section  describes  the  equipment  and 
techniques  used  to  make  measurements  of  the  sensor's  response  to 
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DC  and  AC  magnetic  fields. 


DC  measurements.  Both  the  HeNe  fused  coupler  interferometer 
and  the  HeNe  fringe  pattern  interferometers  were  used  to  make  DC 
magnetic  field  measurements.  The  method  for  quantifying  the 
sensor's  response  was  to  make  measurements  of  output  power  for 
different  current  levels  in  the  Helmholtz  coil. 

To  verify  that  the  sensor  sensitivity  measurements  were 
valid,  a  sensor  arm  with  a  known  response  was  used.  This  new 
sensor  arm  was  a  fiber  wrapped  around  a  5  cm  wide  piezoelectric 
( PZT )  ring  (Veritron,  Model  24-32200-5A,  New  Bedford,  OH).  As 
explained  before,  placing  a  voltage  across  the  PZT  ring  causes 
its  outer  dimension  to  change,  which  would  simulate  a  fiber  being 
stretched  by  a  magnetost r ict ive  material.  The  PZT  ring  was 
tightly  wrapped  with  ten  turns  of  fiber  with  the  two  end  points 
being  glued  to  the  ring.  Different  DC  voltage  levels  were 
applied  to  the  PZT  ring  to  simulate  the  effect  of  DC  magnetic 
fields  on  the  magnetostr ict ive  material  sensor  arms. 

AC  measurements.  The  HeNe  fused  coupler  was  the  only 
interferometer  used  to  make  AC  measurements.  The  equipment  used 
to  make  the  measurements  is  first  reported,  which  is  followed  by 
a  description  of  the  method  used  to  quantify  the  sensor's 
performance . 

The  signals  from  the  optical  power  meter  and  function 
generator  served  as  inputs  to  the  dual-channel  signal  analyzer 
(Bruel  &  Kjaer,  Model  2032,  Marlborough,  MA) ,  which  served  as  a 
digital  storage  oscilloscope  and  spectrum  analyzer.  The 
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information  from  the  signal  analyzer  was  output  via  a  GPIB 
interface  card  (MetraByte,  Model  IE-488,  Taunton,  MA)  to  a 
personal  computer  (Zenith,  Model  Z-248,  St.  Joseph,  MI).  Two 
BASIC  computer  programs,  which  are  included  as  Appendix  D,  were 
written  to  control  the  signal  analyzer  and  store  the  measured 
data.  The  first  program  permitted  the  operator  to  switch  between 
viewing  single  and  dual  displays  of  the  signals'  time-  and 
f reguency-domain  characteristics.  This  program  also  stored  the 
digitized  displays.  The  second  program  was  used  to  store  the 
spectral  information  spanning  23  to  27  Hz. 

The  method  used  to  determine  the  sensor's  sensitivity  was  to 
examine  the  freguency  content  of  the  interferometer's  output. 
Before  the  Helmholtz  coil  was  energized,  the  sensor's  output 
spectrum  was  measured  thirty  times,  averaged,  and  stored  as  a 
baseline  noise  file  in  a  BASIC  data  file.  The  magnetic  field  was 
then  generated  and  the  sensor's  output  spectrum  was  again 
measured  thirty  times,  averaged,  and  stored  as  a  signal  file  in 
another  BASIC  data  file.  The  two  files  were  then  imported  into 
the  LOTUS  123  spreadsheet  program  so  that  the  noise  file  could  be 
subtracted  from  the  signal  file  to  form  the  sensor's  response 
file.  As  was  done  with  the  DC  measurements,  the  PZT  ring  sensor 
arm  w as  used  to  verify  that  the  instrumentation  configuration  and 
data  processing  scheme  used  to  guantify  the  sensor's  response  was 
valid. 

With  the  magnetostrict ive  material  sensor  arms  as  part  of 
the  interferometer,  the  sensor's  response  was  investigated  to 
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ensure  that  the  response  was  due  to  the  magnetostrict ive  effect 
in  the  sensor  arm,  and  not  a  noise  source.  Three  experiments 
were  conducted  to  prove  that  the  sensor  was  operating  correctly. 
The  initial  step  of  the  first  experiment  involved  energizing  the 
Helmholtz  coil  and  examining  the  frequency  content  of  the  output 
signal.  The  output  signal  was  then  examined  when  the  light  was 
only  allowed  to  propagate  in  either  the  sensor  or  reference  arm. 
This  experiment  proved  that  the  sensor's  response  was  due  to  an 
interference  effect,  and  not  some  extraneous  noise  source,  such 
as  cross  talk  between  the  instrumentation's  electrical  cables. 

The  second  experiment  required  replacing  the  sensor  arm  with  an 
uncoated  fiber  and  then  examining  the  sensor's  response.  The 
final  experiment  investigated  the  reproducibility  of  measurements 
taken  with  the  bulk  sensor  arms.  To  ensure  that  any  changes  in 
the  sensor's  performance  were  not  due  to  variations  in  the 
bonding  of  the  fiber  to  the  rod,  the  sensor's  response  was 
recorded  after  the  same  type  of  bulk  sensor  arm  was  constructed 
three  times.  After  each  measurement,  the  fiber  was  removed  from 
the  rod  and  re-bonded  to  the  magnetostrictive  material. 

After  the  method  of  quantifying  the  sensor's  response  was 
validated,  and  the  response  of  the  sensor  was  verified  to  be  the 
result  of  the  external  magnetic  field,  the  sensor  was  evaluated 
with  the  different  sensor  arms.  The  output  spectrum  of  the 
sensor  with  each  different  sensor  arm  was  recorded  for  a  series 
of  different  external  magnetic  field  values. 
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Sensor  Optimization 

The  AC  measurements  made  thus  far  involved  comparing  the 
response  of  the  sensor  with  different  sensor  arms  relative  to  the 
same  series  of  magnetic  fields.  Two  final  experiments  were 
conducted  to  determine  if  the  sensor's  response  could  be  improved 
without  changing  the  sensor  arm.  In  practice,  the  sensor  should 
be  more  sensitive  to  a  small  AC  magnetic  field  superimposed  on  a 
DC  bias  magnetic  field  compared  to  only  a  small  AC  magnetic 
field.  This  optimization  concept  corresponds  to  operating  the 
sensor  in  the  region  of  greatest  slope  depicted  on  the  factional 
change  in  length  vs  the  magnetic  field  plot  shown  in  Figure 
III. 9. 

To  examine  this  hypothesis,  the  sensor's  response  was  tested 
for  a  0.2  G  AC  magnetic  field  with  various  DC  bias  magnetic 
fields.  Both  magnetic  fields  were  generated  at  the  same  time 
with  the  Helmholtz  coil  by  adding  a  DC  current  level  to  the  AC 
current  from  the  function  generator. 

As  the  final  experiment,  the  ultimate  sensitivity  of  the 
sensor  was  evaluated.  The  sensor  was  composed  of  the  sensor  arm 
which  proved  to  be  the  most  responsive  to  the  external  AC 
magnetic  fields.  The  sensor's  response  was  recorded  for 
decreasing  AC  magnetic  fields  that  were  superimposed  on  a 
constant  DC  bias  magnetic  field. 
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In  this  Chapter,  the  experimental  results  and  an  analysis  of 


those  results  are  discussed.  The  outcome  of  the  Helmholtz  coil 
calibration,  censor  arm  fabrication,  interferometer  arrangement, 
sensor  sensitivity  measurements,  and  sensor  optimization  are 
reported . 

Helmholtz  Coil  Performance 

The  first  measurements  recorded  were  the  characteri zat ion  of 
the  Helmholtz  coil.  Once  the  Hall  effect  magnetometer  was 
calibrated,  and  the  ambient  magnetic  field  was  determined,  the 
performance  of  the  Helmholtz  coil  was  evaluated. 

To  calibrate  the  Hall  effect  magnetometer,  the  probe  of  the 
device  was  placed  inside  the  zero  gauss  chamber,  followed  by  the 
95  G  and  978  G  reference  magnets,  respectively.  As  a  result,  the 
magnetometer  measured  0  G,  90  G,  and  980  G,  respectively.  The 
uncertainty  of  the  magnetometer  analog  display  was  approximately 
one  percent.  Some  of  the  calibration  error  may  be  due  to  the 
fact  that  the  reference  magnets  have  not  been  calibrated  for 
several  years.  For  the  purposes  of  establishing  the  ambient 
magnetic  field  flux  density  and  calibrating  the  Helmholtz  coil, 
the  magnetometer  was  deemed  sufficiently  accurate. 

After  the  Hall  effect  magnetometer  was  calibrated,  the 
ambient  magnetic  field  flux  density  inside  the  magnetically 
shielded  chamber  and  around  the  laboratory  was  established. 

Inside  the  test  chamber,  the  DC  magnetic  field  flux  density 
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measured  a  constant  0.215  G.  Around  the  laboratory,  the  DC 
magnetic  field  flux  density  measured  between  0.2  and  1  G.  The  DC 
magnetic  field  flux  density  did  not  vary  with  time  at  any 
location  in  the  laboratory.  The  AC  magnetic  field  flux  density 
was  correspondingly  measured  to  be  a  constant  0  G  both  inside  the 
test  chamber  and  in  various  locations  in  the  laboratory.  The 
magnetically  shielded  chamber  was  not  used  for  four  reasons:  (i) 

the  ambient  DC  magnetic  field  flux  density  did  not  vary  with 
time,  (ii)  the  chamber  did  not  significantly  reduce  the  magnitude 
of  ambient  dc  magnetic  field  flux  density,  (iii)  there  was  no 
detectable  ambient  AC  magnetic  field  flux  density,  and  (iv)  using 
the  chamber  would  add  additional  environmental  noise  due  to  the 
longer  lengths  of  fibers  required. 

Following  the  investigation  of  the  ambient  magnetic  field 
conditions,  the  Helmholtz  coil  was  characterized.  The  first 
experiment  conducted  with  the  Helmholtz  coil  was  to  record  the 
magnetic  field  flux  density  at  the  on-axis  center  point  between 
the  two  coils  while  varying  the  amount  of  AC  and  DC  current 
flowing  through  the  coils.  The  Helmholtz  coil  was  wired  such 
that  the  current  flowed  through  250  turns. 

Figure  V.l  illustrates  the  Helmholtz  coil's  experimental 
performance  compared  to  its  theoretical  performance,  which  was 
predicted  by  equation  11  in  chapter  III.  As  shown,  the  DC 
magnetic  field  flux  density  is  slightly  greater  than  the 
theoretically  predicted  value.  For  example,  at  550  mA,  the  dc 
magnetic  field  flux  density  is  102%  of  its  theoretically 


V-2 


below  the  theoretically  predicted  value.  At  550  roA,  the  measured 
value  was  87%  of  the  predicted  value.  This  discrepancy  can  be 
attributed  to  the  distortion  of  the  25  Hz  sine  wave  introduced  by 
the  power  amplifier.  This  distortion  altered  the  frequency 
content  of  the  signal.  In  particular,  the  distortion  introduced 
frequencies  above  the  400  Hz  detection  limit  of  the  magnetometer. 
At  lower  currents,  the  distortion  was  not  as  severe,  which 
explains  why  the  experimental  and  predicted  values  converge  at 
lower  currents.  Figures  V.2,  V.3,  V.4,  and  V.5  depict  the 
distortion  of  the  sine  wave  caused  by  the  power  amplifier  when 
500  mA  was  drawn  from  the  power  amplifier. 

The  second  experiment  conducted  with  the  Helmholtz  coil  was 
accomplished  to  record  the  magnetic  field  flux  density  at  the  on- 
axis  center  point  between  the  coils  for  different  numbers  of 
turns  of  wire  with  a  constant  current  of  250  mA.  Figure  V.6 
shows  the  result  of  this  experiment.  The  reasons  for  the 
discrepancy  between  the  predicted  and  experimental  values  in  this 
experiment  are  the  same  as  for  the  first  experiment. 

The  final  experiment  conducted  with  the  Helmholtz  coil  was 
accomplished  to  examine  the  spatial  variations  of  the  magnetic 
field  flux  density  generated.  No  variations  in  magnetic  field 
flux  density  could  be  detected  along  the  10.5  cm  centerline 
between  the  coils  for  AC  or  DC  magnetic  fields.  No  variations 
could  be  detected  for  4  cm  in  a  radial  direction  measured  outward 
along  the  center  line  relative  to  the  center  point  between  the 
two  coils. 


V— 4 


Figure  V.2.  Tine-Dor.ain  Display  of  the  Function  Generator  '  r 
Output 
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Figure  V.6.  Magnetic  Field  Flux  Density  General 
Helmhcltz  Coi]  for  Varying  Numbers  : 

Tc  summarize  the  Helmholtz  coil's  performance 
or. eluded  that  the  difference-  between  the  predict*: 
f  an  ideal  Helmholtz  coil  and  the  performance  o: 
n  a  laboratory  was  sufficiently  small  to  be  acccj 


Helmholtz  coil  was  capable  of  simulating  biomagnetic  fields,  and 
consequently,  it  was  used  to  challenge  the  fiber  optic  magnetic 
field  sensor. 

Fabricated  Sensor  Arms 

Following  the  characterization  of  the  Helmholtz  coil,  the 
sensor  arms  were  fabricated.  This  section  reports  the  results  of 
fabricating  the  two  classes  of  sensor  arms:  thin  film  sensor  arms 
and  bulk  material  sensor  arms.  Table  V.l  summarizes  the  sensor 


Table  V.l.  Fabricated  Sensor  Arms 


Sensor  Arm 

Thin/Bulk 

Material 

Film 

Thickness 

Comment 

1 

Bulk 

PZT 

N/A 

Validation  tool 

2 

Thin 

Ni 

9.5  urn 

- 

3 

Thin 

Ni 

4.8  urn 

- 

4 

Thin 

Ni 

9.5  um 

Spiral  design 

5 

Thin 

Ni 

8  um 

9  lengths 

6 

Thin 

Co 

0.3  um 

- 

7 

Thin 

Fe 

<  0.1  um 

Possible  Cu  ccriton 

8 

Bulk 

Ni 

- 

As  received 

5 

Bulk 

Ni 

- 

Therm.a  1  ly  a rne a  1  e  d 

10 

Bulk 

Ni 

- 

Thermally  annealed 

with  B  field 

11 

Bulk 

Ni 

- 

Thinned  fiber 

12 

Bulk 

Ni 

- 

3  lengths 
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arms  fabricated  for  this  research  project. 


Thin  film  sensor  arms.  Thin  films  of  cobalt,  iron  and 
nickel  were  deposited  on  the  fiber  optic  cables.  It  was  an 
objective  to  fabricate  the  fibers  with  equal  metal  thicknesses. 
This  part  reports  the  difficulties  encountered  which  prevented 
the  designed  sensor  arms  from  being  fabricated,  and  discusses  the 
results  of  characterizing  the  metal  coated  fibers. 

Cobalt  and  iron  were  difficult  metals  to  deposit.  Three  and 
one  half  hours  of  sputtering  cobalt  resulted  in  a  film  thickness 
of  only  0.3  um.  Sputtering  iron  for  one  and  one  half  hours 
produced  a  film  thickness  that  was  so  thin  that  it  could  not  be 
measured;  it  was  sufficiently  thin  to  be  partially  transparent. 
Based  on  the  sensitivity  limit  of  the  profilometer  when  used  to 
determine  a  film  thickness  on  a  glass  microscope  slide,  an 
estimate  of  the  iron  film  thickness  was  determined  to  be  less 
than  0.1  um.  The  sputtering  system  may  have  not  been  operating 
correctly  when  the  iron  layer  was  deposited  since  some  copper 
from  the  sputtering  target  base  may  have  been  deposited  with  the 
iron.  Copper  content  in  the  deposited  film  was  suspected  due  to 
the  bronze  tint  of  the  deposited  metal.  Additional  cobalt  and 
iron  coated  fibers  were  not  produced  due  to  the  disappointing 
sensor  results  of  the  nickel  sputtered  fibers. 

As  a  softer  metal  than  cobalt  and  iron,  nickel  was  more 
easily  sputtered.  A  9.5  um  thick  layer  of  nickel  was  deposited 
in  eight  hours.  Based  on  the  relatively  quick  deposition  rate 
and  Figure  III. 9,  which  indicated  that  nickel  experiences  the 
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largest  magnetostriction  effects,  many  different  nickel  fibers 
were  fabricated  as  reflected  in  Table  V.l. 

To  test  the  effect  of  sensor  arm  metal  thickness  on  sensor 
performance,  two  fibers  of  different  nickel  film  thickness  were 
created.  The  first  fiber  was  coated  with  a  9.5  urn  thick  layer  of 
nickel.  A  second  fiber  was  sputtered  for  four  hours,  which 
resulted  in  a  film  thickness  of  4.8  urn.  The  fact  that  the  filr. 
thicknesses  were  not  exactly  proportional  to  the  duration  of  the 
metal  deposition  is  not  surprising  since  sputtering  is  not  a 
linear  deposition  process. 

To  create  fibers  with  different  geometrical  patterns,  the 
fibers  were  partially  masked  during  metal  deposition.  The  seri- 
cylindrical  patterns  could  not  be  fabricated  because  the  tape 
used  to  mask  the  fiber  could  not  be  removed  without  breaking  the 
fiber.  Only  one  fiber  was  created  with  a  spiral  pattern  due  to 
the  handling  difficulties  of  wrapping  a  0.05  mm  wire  around  the 
bare  fiber. 

Three  attempts  were  made  to  thermally  anneal  a  sputtered 
nickel  fiber,  but  each  time  the  fiber  broke  inside  the  even.  As 
the  fiber  cooled,  the  thermal  expansion  coefficient  mismatch  must 
have  been  sufficient  to  break  the  fiber.  For  example,  when  the 
temperature  of  a  nickel  sample  is  raised  from  293  to  400c  K,  the 
change  in  length  relative  to  a  nominal  length  is  0.15%  (35:225). 

For  a  sample  of  silica  experiencing  the  same  change  in 
temperature,  the  change  in  length  relative  to  a  nominal  length  is 
0.12%  (35:290). 
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The  broken  thermally  annealed  fibers  do  indicate  that  the 
metal  adhered  well  to  the  fiber.  The  results  of  the  optical 
microscope  inspection  indicate  that  metal  film's  surface  was 
smooth  and  uniform.  Figure  V.7  illustrates  an  example  of  the 
nickel  surface  which  also  shows  the  notch  from  the  spiral  design 
on  the  fiber.  The  surface  quality  of  the  iron  and  cobalt  coated 
fibers  was  similar  to  the  nickel  fiber.  Various  metalized  fibers 
were  cleaved  to  determine  the  uniformity  of  the  deposited  layers. 
The  metal  coating  was  not  visible  from  the  end  view  of  the  fiber, 
which  was  probably  due  to  the  thin  coating  thickness  and  uneven 
cleaving  of  the  metal  coated  fiber. 

Bulk  material  sensor  arms.  Bulk  material  sensor  arms  were 
fabricated  as  summarized  in  Table  V.l.  The  fiber  in  one  sensor 


Figure  V.7.  Nickel  Coated  Fiber  Surface 


Showing  the  Spiral  Notch 
Pattern 
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arm  was  thinned  as  shown  when  compared  to  a  full-sized  fiber  in 
Figure  V.8.  Based  on  this  photograph,  the  fiber's  diameter  was 
reduced  approximately  48%,  which  corresponds  to  a  23%  reduction 
in  cross-sectional  area. 

Overall,  fabricating  the  sensor  arms  from  bulk  materials  was 
much  easier  than  RF  sputtering  due  to  the  fiber  handling 
considerations  and  metal  bonding  techniques.  The  RF  sputtering 
system  was  not  an  ideal  system  to  apply  metal  to  fibers.  The 
most  serious  problem  with  the  sputtering  system  was  the  amount  of 
time  required  to  deposit  layers  greater  than  10  um  thick.  The 
sputtering  system  also  suffered  numerous  mechanical  malfunctions 
that  not  only  delayed  processing,  but  occasionally  destroyed  the 
samples . 


Figure  V.8.  Comparison  of  the  Thinned  Fiber 
with  a  Full  Size  Fiber 
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Interferometer  Parameters 


Once  the  sensor  arms  were  fabricated,  all  the  components 
were  available  to  assemble  the  interferometers.  This  section 
reports  the  results  of  fabricating  the  interferometers  in  two 
parts.  The  first  part  describes  the  results  of  the  procedures 
developed  to  work  with  fiber  optic  cables.  The  results 
associated  with  fabricating  the  interferometers  are  included  in 
the  second  part. 

Fiber  optic  cable  results.  This  part  documents  the 
consequences  of  the  techniques  developed  to  cleave  the  ends  of 
fibers,  couple  light  into  the  end  of  the  fiber,  and  make  fiber 
splices.  Except  when  otherwise  noted,  the  maximum  uncertainty  in 
the  optical  power  measurements  was  approximately  five  percent. 

The  initial  fiber  cleaver  obtained  did  not  produce  adequate 
quality  cleaves.  When  viewed  with  the  fusion  splicer  optics,  the 
ends  of  the  fibers  were  observed  to  be  not  smooth  and 
perpendicular  to  the  optic  axis  of  the  fibers.  The  poor  quality 
cleaved  edges  were  attributed  to  the  finite  lifetime  of  the 
carbide  edge  and  bending  mechanism  of  the  fiber  cleaver.  A  newer 
fiber  cleaver  was  borrowed  from  Dale  Stevens  of  the  Electronic 
Technology  Laboratory.  This  cleaver  produced  high  quality 
cleaved  ends  practically  every  time. 

High  quality  cleaved  fiber  ends  were  necessary  to  couple 
light  into  the  end  of  the  fiber.  Using  the  method  outlined  in 
the  previous  chapter,  400  uW  from  a  1.914  mW  laser  was  coupled 
into  the  fiber.  This  loss  of  7  dB  was  due  to  two  following 


V-13 


reasons:  i)  mismatch  in  the  area  of  the  focused  laser  spot  size 

and  the  fiber's  core  and  ii)  mismatch  in  the  numerical  aperture 
(NA)  between  the  microscope  objective  and  the  optical  fiber.  As 
a  first  level  approximation,  the  laser  spot  size  is  considered  a 
uniform  circular  beam  of  light  whose  area  is  4.91  x  10‘8  m2. 

After  passing  through  the  20X  microscope  objective,  the  focused 
spot  size  would  then  have  an  area  of  1.23  x  10‘10  m2.  With  a 
fiber's  core  area  of  only  9.62  x  10’12  m2,  the  expected  loss  would 
be  11  dB.  Ignoring  the  gaussian  nature  of  the  laser  beam  and  its 
actual  interaction  with  optics  accounts  for  the  difference 
between  the  achieved  light  coupling  and  the  expected  light 
coupling.  The  numerical  aperture  (NA)  affects  the  light  coupling 
since  it  describes  the  solid  angle  of  light  focused  from  the 
microscope  and  the  solid  angle  of  light  that  will  be  accepted 
into  the  optical  fiber  (29:51).  The  NA  of  the  microscope 
objective  was  0.4,  and  the  optical  fiber's  NA  was  0.1.  This 
mismatch  in  NA  indicates  that  a  portion  of  the  light  focused  from 
the  microscope  objective  will  enter  the  fiber  at  an  angle  too 
large  as  measured  from  the  optic  axis  to  be  propagated  along  the 
fiber.  The  only  undesirable  feature  of  the  fiber  light  coupler 
was  the  plastic  screw  micrometers  which  required  a  fine  touch 
during  the  final  adjustment  process. 

Low  loss  splices  also  required  high  quality  cleaved  fiber 
ends.  The  fusion  splicer  proved  to  be  inadequate  to  form 
splices.  The  available  fusion  splicer  had  only  micrometer 
adjustments  along  the  optical  axis  of  the  two  fibers.  The  lack 
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of  adjustments  in  the  outward  radial  directions  prevented 
reproducible  low  loss  splices,  with  most  fusion  splices  resulting 
in  over  5  dB  of  loss.  The  fusion  splicer  would  be  adequate  for 
multi-mode  fiber  splices  since  the  fiber  core  size  is  much  larger 
than  in  single-mode  fibers.  The  glass  ferrule  splice  tube 
established  itself  as  the  superior  technology  needed  for  making 
single-mode  optical  fiber  splices.  Using  the  glass  tubes, 
splices  with  less  than  0.5  dB  of  loss  were  formed.  A  further 
advantage  of  using  the  glass  ferrules  is  that  the  splices  were 
temporary.  However,  as  the  glass  tubes  were  reused,  index 
matching  gel  bubbles  accumulated  in  the  tube  causing  higher 
losses  of  1  to  2  dB. 

In  summary,  the  initial  techniques  required  to  use  optical 
fibers  were  very  time  consuming.  However,  once  the  techniques 
were  refined  and  practiced,  optical  fibers  could  be  treated  as 
any  other  optical  component  in  +-he  interferometer. 

Interferometer  Performance.  The  three  interferometer 
configurations  were  assembled  once  the  techniques  to  work  with 
the  fibers  were  developed.  This  part  discusses  the  results  of 
the  characterization  of  the  remaining  optical  equipment  and 
describes  the  performance  of  the  three  interferometer 
configurations . 

The  optical  power  of  the  HeNe  laser  measured  with  the 
radiometer  was  1.914  mW.  The  laser's  optical  power  varied  less 
than  one  percent  of  this  value  after  it  was  allowed  to  stabilize 
for  thirty  minutes.  HeNe  laser  output  stability  problems  were 
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experienced  when  reflections  from  the  end  of  the  optical  fiber 
were  directed  back  into  the  laser.  This  problem  was  corrected  by 
providing  a  sufficient  distance  between  the  HeNe  laser  and  the 
optical  fiber.  Since  the  laser  could  not  be  perfectly  aligned 
with  the  optical  fiber,  the  reflection  deviated  from  the  axis  of 
the  laser's  output  beam.  When  the  distance  between  the  laser  and 
optical  fiber  was  greater  than  0.3  m,  the  reflected  laser  power 
did  not  couple  back  into  the  laser. 

In  addition  to  optical  power,  the  laser's  output  was 
investigated  for  polarization  and  coherence.  The  polarization  of 
the  HeNe  laser's  output  was  measured  to  be  linear.  With  a 
polarizer  placed  with  its  transmission  axis  parallel  to  the 
laser's  polarization,  the  output  power  was  measured  to  be  1.913 
mW.  When  the  laser's  polarization  was  perpendicular  to  the 
transmission  axis  of  the  polarizer,  the  output  power  was  measured 
to  be  0.001  mW.  The  exact  coherence  length  was  not  determined 
because  the  longest  difference  in  path  lengths  was  less  than  1  r. 
Observed  interference  effects  indicated  that  the  coherence  length 
of  the  laser's  output  must  be  greater  than  1  m. 

The  final  optical  component  to  be  characterized  was  the 
fused  coupler.  The  coupling  ratio  agreed  with  the  manufacturers 
reported  figure  of  2.91  dB  within  ten  percent  (34).  An 
experiment  was  conducted  to  measure  the  coupling  ratio  by 
inserting  light  into  port  1  of  the  fused  coupler,  and  then 
measuring  the  output  power  in  ports  2  and  3.  The  output  power 
was  measured  to  155  uW  and  180  uW  for  ports  2  and  3, 
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respectively.  These  values  indicate  a  coupling  ratio  of  -2.70 
dB.  The  insertion  ratio  was  more  difficult  to  measure  since  it 
involved  other  variables  such  as  splicing  efficiency  or  input 
light  coupling.  However,  the  insertion  loss  of  the  device  was 
probably  close  to  the  stated  value  based  on  the  average  output 
power  of  the  interferometer.  More  information  on  the  fused 
coupler  is  included  in  the  following  discussion  concerning  the 
interferometer  configurations. 

After  all  the  optical  components  were  evaluated  separately, 
the  HeNe  fused  coupler  interferometer  was  assembled  as  shown  in 
Figure  V.9.  Table  V.2  summarizes  the  measured  power  budget  for 
the  HeNe  fused  coupler  interferometer.  The  value  associated  with 
the  optical  power  measured  by  the  detector  is  an  average  value 
(the  minimum  and  maximum  values  were  1  uW  and  67  uW) . 


Figure  V.9.  HeNe  Fused  Coupler 
Interferometer 
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Table  V.2.  Example  of  the  Power  Budget  for  the  HeNe  Fused 
Coupler  Interferometer  Configuration 


Component 

Attenuation 

Power 

Laser 

0  dB 

1.91  mw 

Beam  Splitter 

-0.05  dB 

1.89  mw 

Fiber  Light  Coupler 

-11  dB  (2) 

75  +  75  =  150 

UW 

1  m  of  Optical  Fiber 

-0.01  dB  (2) 

74.8  +  74.8  = 

149.6  uV 

Splice 

-3  dB  (2) 

37.5  +  37.5  = 

7  5  UW 

Fused  Coupler 

-0.3  dB 

7  0  uW 

Detector 

-0.2  dB 

3  4  uW 

Unfortunately,  when  the  interferometer  was  completely 
assembled,  the  output  power  measured  by  the  detector  was  not 
constant.  However,  with  light  only  propagating  in  one  arm  of  the 
interferometer,  the  output  power  was  constant  within  ten  percent. 
As  soon  as  light  propagated  in  both  arms  to  complete  the 
interferometer  configuration,  the  output  power  was  no  longer 
constant.  For  example,  with  30  uW  of  light  inserted  into  each 
arm  of  the  interferometer,  the  output  power  measured  by  the 
optical  power  meter  would  vary  from  1  to  60  uK  with  a  randomly 
changing  frequency  of  a  few  hertz.  It  was  observed  that  when  the 
power  from  one  output  leg  of  the  fused  coupler  was  a  maximum,  the 
power  from  the  other  output  leg  of  the  fused  coupler  was  a 
minimum.  Combining  the  power  from  both  legs  yielded  a  more 
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constant  total  output  power  value,  which  now  varied  to  within 
twenty  percent  of  its  nominal  value. 

The  problem  of  excessive  noise  in  the  fused  coupler  was 
initially  attributed  to  be  a  defect  in  the  fused  coupler.  To 
test  this  hypothesis,  eight  different  samples  of  the  same  model 
of  fused  coupler  were  inserted  into  the  interferometer;  all  the 
couplers  produced  identical  results.  Consequently,  the  GaAs 
fused  coupler  interferometer  was  then  constructed  to  test  whether 
the  excessive  noise  was  due  to  any  particular  component  in  the 
HeNe  fused  coupler  interferometer.  The  different  optical  source 
wavelength  requirement  dictated  that  the  GaAs  fused  coupler 
interferometer  be  constructed  with  a  different  optical  source, 
optical  fiber,  and  fused  coupler  compared  to  the  HeNe  fused 
coupler  interferometer.  However,  the  same  excessive  noise 
occurred  in  the  GaAs  fused  coupler  interferometer.  As  a  final 
check  of  the  equipment,  Dr.  Dale  Sarrafzaeh  of  the  American 
Research  Corporation  of  Virginia,  who  is  also  researching  the 
fiber  optic  magnetic  field  sensor  concept,  was  contacted  (36). 

Dr.  Sarrafzaeh  indicated  that  the  same  effect  occurred  in  the 
interferometer  that  he  constructed.  This  result  indicates  that 
the  phenomena  is  not  equipment  dependent,  but  rather  due  tc  the 
inherent  physics  of  the  fused  coupler  Mach-Zehnder  interferometer 
configuration . 

Two  possible  explanations  for  the  excessive  noise  were 
offered  by  two  engineers'  in  industry  who  are  familiar  with  the 
fused  coupler .  Dr.  Sarraf zaeh  attributed  the  excessive  noise  to 
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the  fact  that  the  interferometer  was  very  sensitive  to  any 
environmental  perturbations  such  as  infinitesimal  temperature 
variations  and  vibrations.  In  addition,  Brad  Sousa  of  York 
Technology,  the  firm  that  manufactures  the  HeNe  fused  couplers, 
suggested  a  different  explanation  (37).  He  stated  that  the 
excessive  noise  was  due  to  interference  between  the  different 
polarizations  of  light  in  each  input  of  the  coupler.  Mr.  Sousa 
stated  that  the  excessive  noise  would  not  occur  if  the 
interferometer  incorporated  polarization-preserving,  single-mode 
optical  fiber  and  fused  couplers.  Due  to  the  unavailability  of 
the  polarization-preserving  components,  this  hypothesis  was  not 
tested . 

Since  the  excessive  noise  in  the  fused  couplers  could  not  be 
eliminated,  the  HeNe  fringe  pattern  interferometer  was 
constructed.  In  this  interferometer  configuration,  the  output 
beams  are  recombined  using  bulk  optics  instead  of  a  fused 
coupler.  The  HeNe  fringe  pattern  interferometer,  as  -hewn  in 
Figure  IV. 3,  was  produced.  Reflections  from  the  front  and  back 
side  of  the  beam  splitter  caused  a  secondary  fringe  pattern  tc  to 
formed.  This  secondary  fringe  pattern  produced  smaller  light  and 
dark  bands  within  the  primary  fringe  pattern. 

The  secondary  fringe  pattern  could  have  been  removed  by  two 
methods.  The  first  technique  is  to  use  a  beam  splitter  that  was 
slightly  wedged-shaped  such  that  the  secondary  fringe  pattern 
would  form  at  an  offset  relative  to  the  primary  fringe  pattern. 
This  technique  was  not  used  since  a  wedged-shaped  beam  spl itter 
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was  not  available.  A  second  option  for  correcting  this 
aberration  is  to  utilize  a  pellicle,  or  thin  membrane  bear, 
splitter.  By  reducing  the  distance  between  the  front  and  back 
surfaces  of  the  beam  splitter,  the  secondary  fringe  pattern  is 
minimized.  The  pellicle  method  was  rejected  since  tne  pellicle 
is  very  sensitive  to  air  movement  and  vibration  which  would  have 
added  another  noise  source  to  the  sensor. 

The  existing  fringe  pattern  presented  two  sources  of 
uncertainty  in  the  sensor.  First,  the  fringe  pattern  drifted 
with  time.  As  the  fringes  swept  by  the  slit,  the  output  power 
measured  by  the  radiometer  varied  from  a  maximum  to  minimum  va lur 
with  a  randomly  changing  frequency  of  a  few  hertz.  Second,  each 
time  a  sensor  arm  was  changed,  a  new  fringe  pattern  had  to  be 
formed.  The  fringe  pattern  could  not  be  exactly  reproduced  duo 
to  the  slight  curve  and  width  of  each  line  in  the  fringe  patter:.. 
As  a  consequence,  the  response  of  the  sensor  would  be  slightly 
different  for  each  of  the  different  fringe  patterns  because  the 
slit  would  be  sampling  a  different  amount  of  each  fringe  pattern. 

In  summary,  the  three  interferometer  configurations  were 
assembled.  The  GaAs  fused  coupler  interferometer  was  only  u  •.  . : 
to  examine  whether  the  excessive  noise  contributed  i  rom  the  :  -  ■. 
coupler  was  dependent  on  the  model  of  the  fused  cout.  1c:  .  e 
HeNe  fringe  pattern  was  assembled,  but  due  to  the  uncertainty 
involved  with  the  fringe  pattern,  was  only  used  to  in  an  attest t 
examine  the  sensor's  response  to  DC  magnetic  fields.  Although 
the  HeNe  fused  coupler  interferometer  suffered  from  excessive 
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noise  associated  with  the  fused  coupler,  this  interferometer  was 
used  to  measure  the  sensor's  response  to  DC  and  AC  magnetic 
fields . 

Sensitivity  Measurements 

After  the  interferometers  were  assembled,  the  sensor's 
response  for  DC  and  AC  magnetic  fields  was  quantified.  After 
describing  the  DC  measurements,  the  verification  of  the 
measurement  method  for  the  AC  sensor  response  is  detailed.  This 
section  concludes  with  the  results  associated  with  the  sensor's 
response  to  several  magnetic  field  strengths  when  the  sensor 
incorporated  of  the  different  sensor  arms. 

DC  measurements.  DC  measurements  were  initially  attempted 
with  the  PZT  ring  as  the  sensor  arm  in  both  the  HeNe  fused 
coupler  interferometer  and  the  fringe  pattern  interferometer. 
Both  steady  voltage  levels  and  ramping  voltages  were  applied  to 
PZT  ring.  Due  to  the  excessive  noise  from  the  fused  coupler  and 
the  fringe  pattern  drift,  the  external  stimulus  could  not  be 
detected.  With  the  magnetostrict ive  material  sensor  arms  placed 
in  both  type  of  interferometers,  the  same  negative  result  was 
observed.  This  result  indicates  that  the  fiber  optic 
interferometers  constructed  for  this  research  project  are  not 
appropriate  for  detecting  DC  magnetic  fields. 

AC  measurements.  Prior  to  measuring  the  response  of  the 
sensor  to  AC  magnetic  fields,  the  measurement  method  was 
validated  using  the  PZT  ring  sensor  arm.  After  the  measurement 
method  was  confirmed,  and  the  sensor's  response  with 
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magnetostrictive  coating  sensor  arms  was  verified  to  be  from  the 
magnetic  field,  the  sensor's  performance  was  evaluated  for 
different  sensor  arms. 

Due  to  the  uncertainty  involved  with  the  fringe  pattern 
interferometer,  only  the  HeNe  fused  coupler  interferometer  was 
used  to  accomplish  the  AC  measurements.  The  sensor  was  first 
characterized  with  the  PZT  ring  as  the  sensor  arm.  Figure  V.10 
illustrates  the  averaged  spectrum  of  the  interferometer's  output 
with  no  voltage  applied  to  the  PZT  ring.  Figure  V.ll  depicts  the 
averaged  spectrum  of  the  sensor's  output  with  10  Vp-p  signal 
applied  to  the  PZT  ring,  and  Figure  V.12  shows  the  difference 
spectrum  which  is  the  result  of  subtracting  the  noise  spectrum 
from  the  response  spectrum.  As  shown  in  Figure  V.12,  the 
excitation  signal  at  25  Hz  and  its  harmonics  were  easily 
detected.  The  harmonics  are  due  to  the  fact  that  the  PZT  ring 
induced  a  phase  change  difference  greater  than  one  wavelength 
with  respect  to  the  optical  source.  Figure  V.13  illustrates  the 
sensor's  difference  spectrum  for  different  voltage  levels  applied 
to  the  PZT  ring.  Figure  V.14  depicts  a  plot  of  the  amplitude  of 
the  interferometer's  difference  spectral  response  at  25  Hz  for 
decreasing  values  of  the  externally  applied  voltage.  Figure 
V.14  verifies  that  the  amplitude  of  the  difference  spectrum  at  25 
Hz  is  a  direct  indication  of  the  strength  associated  with  an 
external  stimulus.  The  change  in  the  radius  of  the  PZT  ring  as  a 
function  of  the  externally  applied  voltage  signal  level  is  given 
by  (38 : 19)  : 
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Figure  V.10.  Output  Spectrum  of  the  HeNe  Fused  Coupler 


Interferometer  with  the  PZT  Ring  Sensor  Arm  but 
No  Externally  Applied  Voltage  Signal 
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a  25  Hz,  10  Vp-p  Externally  Applied  Voltage  Signal 
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Figure  V. 12 . 


HeNe  Fused  Coupler  Interferometer's  Output  Noise 
Spectrum  (Figure  V.10)  Subtracted  from  the  Signal 
Spectrum  (Figure  V.ll) 
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Figure  V.13.  HeNe  Fused  Coupler  Interferometer  Difference 
Response  for  Different  Applied  Voltages  to 
the  PZT  Ring  Sensor  Arm 
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Figure  V . 14 . 


HeNe  Fused  Coupler  Interferometer  Difference 
Spectral  Component  Magnitude  at  25  Hz  for 
Different  Applied  Voltage  Signals  to  the  PZT 
Ring  Sensor  Arm 


(112) 


Ar  =  d3JV 

where  Ar  =  change  in  radius  of  the  ring 
d33  =  piezoelectric  constant 
V  =  applied  voltage. 

Therefore,  the  change  in  length  of  the  PZT  sensor  arm  is  given 
by: 

Al  =  2tt  Ar  =  27rd33V  (113). 

Equation  113  states  that  the  sensor's  response  in  Figure  V.14  is 
representative  of  a  linear  input  response. 

In  addition  to  verifying  the  measurement  scheme  with  the  PZT 
ring  sensor  arm,  the  magnetostrictive  material  sensor's  response 
was  investigated  to  ensure  that  the  response  was  due  to  the 
magnetostrictive  effect  in  the  sensor  arm  and  not  some  other 
source.  Two  experiments  were  conducted  to  prove  that  the  sensor 
was  operating  correctly. 

The  first  experiment  involved  energizing  the  Helmholtz  coil, 
which  produced  a  sensor  response  similar  to  Figure  V.ll.  When 
the  laser  light  was  only  allowed  to  propagate  in  one  arm  at  a 
time,  the  output  signal  at  25  Hz  became  indistinguishable  from 
the  noise,  which  indicated  that  the  response  was  due  to  an 
interference  effect  and  not  some  extraneous  source,  such  as 
cross-talk  between  electrical  cables. 

The  second  experiment  was  implemented  by  replacing  the 
sensor  arm  with  an  uncoated  fiber  and  measuring  the  sensor's 
response.  Initially,  the  sensor  responded  as  depicted  in  Figure 
V . 15 .  This  response  was  attributed  to  the  small  vibrations 
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Figure  V. 15 . 


HeNe  Fused  Coupler  Interferometer  Difference 
Response  Spectrum  Due  to  Low-level,  25  Hz 
Vibrations 
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generated  by  the  Helmholtz  coil  that  are  similar  to  transformer 
hum.  Once  the  Helmholtz  coil  was  suspended  above  the  air  table, 
the  vibrational  response  disappeared.  This  experiment 
illustrated  the  extreme  sensitivity  of  the  sensor  to 
environmental  noise.  Based  on  these  two  experiments,  it  was 
concluded  that  the  sensor's  response  was  due  to  the 
magnetostrictive  effect  of  the  sensor  arm. 

After  the  sensor's  response  was  demonstrated  to  be  directly 
related  to  the  external  magnetic  field,  the  sensor's  response  was 
evaluated  for  both  thin  film  coatings  and  buik  material 
configured  sensor  arms.  Unfortunately,  none  of  the  sputtered 
thin  film  sensor  arms  produced  a  measurable  response,  even  with 
the  maximum  AC  magnetic  i^wiu  flux  density  of  11.8  G. 

The  failure  of  the  sputtered  fibers  to  respond  to  the  AC 
magnetic  field  was  attributed  to  the  thinness  and  adhesion  of  the 
metal  coating  on  the  fiber  not  generating  enough  force  to 
^compress  or  stretch  the  fiber.  The  thickest  sputtered  film  (S.5 
urn)  coated  a  125  urn  diameter  silica  fiber.  Using  equation  17 
from  the  theory  discussion  of  fiber  optic  cables  in  Chapter  III, 
the  force  required  to  change  the  fiber's  length  is  given  by: 

F  =  hLEA  (17) 

where  F  =  magnitude  of  the  applied  force 
hL  =  change  ir.  length 
L0  =  initial  length  (0.1  m) 

E  =  Young's  modulus  of  the  silica  (6.89  x  10’c  N/nk) 

A  =  cross  sectional  area  of  the  fiber  (1.23  x  10'8  m2)  . 
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To  induce  a  quarter  wavelength  equivalent  change  in  the  fiber's 
length  (158  nm)  requires  the  magnetostrictive  material  to  exert  a 
force  of  1.34  mN.  Although  1.34  mN  is  not  a  large  force,  it  was 
more  than  the  thin  metal  film  coatings  could  generate. 

Fortunately,  the  nickel  rod  generated  a  sufficient  force 
level  for  the  sensor  to  detect,  which  allowed  a  series  of 
experiments  to  be  conducted.  Figure  V.16  depicts  an  example  of 
the  difference  response  spectrum  of  the  sensor  with  the  as 
received  nickel  rod  for  several  magnetic  field  levels.  The  first 
and  second  derivatives  of  the  response  curves  provided  no  further 
information  due  to  the  noise  in  the  signal. 

Before  attributing  any  changes  in  the  sensor’s  response  with 
different  bulk  sensor  arms,  a  reproducibility  study  was 
conduced.  The  sensor's  response  was  recorded  after  the  same 
sensor  a rm  was  bonded  three  times.  Figure  V.17  shows  that  the 
sensor's  response  for  each  trial  was  similar  for  small  AC 

1 

magnetic  flux  densities,  which  indicates  that  changes  a  the 
sensor's  response  are  not  due  to  the  bonding  techniques  used  in 
the  fabrication  of  the  sensor  arm. 

After  the  effect  of  the  fiber-nickel  bonding  technique  was 
investigated,  the  influence  associated  with  thermal  and  magnetic 
field  processing  of  the  nickel  rod  was  determined.  The  nickel 
rod  was  first  characterized  as  it  was  received,  after  thermal 
annealing,  and  after  thermae  annealing  in  the  presence  of  an 
external  magnetic  field.  Figure  V.18  illustrates  the  results  of 
that  experiment.  Thermal  annealing  without  an  external  magnetic 
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Figure  V.17.  Spectral  Response  of  the  Sensor  for  Three 
Similarly  Fabricated  Nickel  Sensor  Arms 
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Figure  V.18.  Effect  of  Thermal  Annealing  with  and  without  an 
External  Magnetic  Field  on  the  Sensor's 
Performance 
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field  in  addition  to  the  earth's  ambient  magnetic  field  did  not 
have  an  effect  on  the  performance  of  the  sensor.  This  result  was 
nor  unexpected  since  the  rod  was  heated  in  its  formation. 

However,  when  the  rod  was  cooled  in  the  presence  of  a  1000  G 
external  magnetic  field  flux  density,  the  magnetic  domains 
partially  aligned  with  this  external  field,  which  resulted  in  an 
improved  magnetostriction  effect. 

The  next  experiment  was  conducted  to  determine 
experimentally  the  effect  of  the  optical  fiber's  geometry  on  the 
sensor's  response.  The  two  sensor  arms  evaluated  in  this 
experiment  were  the  thinned  fiber  sensor  arm  and  the  sensor  arm 
incorporating  three  lengths  of  fiber.  Figure  V.19  illustrates 
the  results  of  this  experiment.  Thinning  the  fiber  did  not 
affect  the  sensor's  performance.  When  the  bulk  nickel  rod 
experienced  magnetostriction,  it  must  have  generated  enough  force 
to  compress  the  fiber  regardless  of  the  fiber's  cross-sectional 
area.  The  sensor's  performance  improved  with  longer  optical 
fiber  lengths  as  the  theory  predicted. 

Figures  V.18  and  V.19  proviae  information  concerning  the 
sensor's  response  and  measurement  methods.  For  small  magnetic 
field  flux  densities  (<  4  G) ,  the  sensor's  response  with  the 
ragnctostrictive  sensor  arms  was  similar  to  the  sensor's 
performance  with  the  PZT  ring  sensor  arm.  This  fact  indicates 
that  the  magnetostri ct i ve  effect  was  approximately  linear  icr 
small  magnetic  flux  densities.  This  situation  corresponds  to 
operating  on  a  small  linear  region  on  the  fractional  change  in 
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Figure  V.19.  Effect  of  the  Optical  Fiber's  Geometry  on  the 
Sensor's  Response 


V  -  3  7 


length  vs  magnetic  flux  density  curve  (Figure  III. 9).  Figures 
V.18  and  V.19  also  indicate  that  the  sensor's  response  saturates 
above  4  G.  From  Figure  III.  9,  this  saturation  behavior  is 
attributed  to  the  measurement  scheme,  and  not  the 
magnetostrict ive  effect.  Figure  V.17  also  illustrates  a  flaw  in 
the  measurement  scheme.  Each  set  of  data  corresponding  to  a 
curve  in  Figures  V.17,  V.18,  and  V.19  should  have  been  recorded 
several  times  and  averaged  to  yield  a  more  representative  curve 
for  the  test  conditions.  This  observation  was  not  made  until  the 
nickel  rod  had  been  irrevocably  processed,  which  prevented 
additional  data  from  being  gathered. 

Sensor  Optimization 

The  final  experiments  conducted  involved  improving  the 
sensor's  performance  without  altering  the  fiber  or  the 
magnetostrictive  material.  Instead,  the  effect  of  a  DC  bias 
magnetic  field  flux  density  on  the  sensor's  response  to  an  AC 
magnetic  field  was  investigated.  The  sensor  arm  used  for  this 
experiment  contained  three  lengths  of  optical  fiber.  Figure  V.20 
shows  the  effect  of  different  DC  bias  magnetic  flux  densities  on 
the  sensor's  response  for  an  AC  magnetic  field  flux  density  of 
0.2  G.  Although  the  curve  is  not  smooth,  the  trend  indicates 
that  the  sensor  is  more  sensitive  under  the  influence  of  a  larger 
DC  bias  magnetic  field  flux  density.  From  the  data,  it  cannot,  be 
determined  if  tnere  is  an  optimum  value  of  DC  bias  magnetic  flux 
density.  If  a  larger  DC  magnetic  flux  density  could  have  been 


generated,  the  peak  sensitivity  in  the  vicinity  of  10  G  might 


Figure  V.2 0.  Effect  of  a  DC  Bias  Magnetic  Flux  Density  on  the 
Sensor's  Response  for  a  Given  AC  Magnetic  Flux 
Density 
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have  been  identified  as  a  local  or  absolute  maximum. 

The  final  sensor  characterization  experiment  involved 
determining  the  smallest  AC  magnetic  field  flux  density  that 
could  be  detected.  The  optimum  configuration  of  the  sensor 
utilized  the  nickel  rod  sensor  arm  that  had  been  thermally 
annealed  in  the  presence  of  an  external  magnetic  field  (1000  G) 
with  three  lengths  of  optical  fiber  attached  tc  it.  A  DC  bias 
magnetic  field  flux  density  of  0.428  G  was  then  applied  to  the 
sensor  arm.  With  only  10  turns  of  the  Helmholtz  coil  energized 
(instead  of  250),  this  situation  represented  the  largest  DC 
magnetic  field  flux  density  that  could  be  generated.  Only  10 
turns  of  the  H<_lm'noltz  cull  were  energised  to  produce  a  small  AC 
magnetic  flux  density.  Figure  V.21  depicts  the  sensor's  response 
for  decreasing  values  of  the  AC  magnetic  field  flux  density.  The 
smallest  AC  magnetic  field  flux  density  that  could  be  generated 
and  observed  was  10  mG,  which  was  5  dB  above  the  average  noise 
level . 

In  summary,  this  research  effort  was  successful  in 
characterizing  certain  aspects  of  the  sensor's  performance. 
Although  the  sensor  did  not  detect  DC  or  AC  magnetic  fields  with 
RF  sputtered  sensor  arms,  the  sensor  did  detect  AC  magnetic 
fields  using  a  bulk  nickel  sensor  arm  configuration.  Using  this 
one  bulk  sensor  arm,  it  was  determined  that  the  sensor  was  more 
sensitive  when  the  magnetostrict ive  material  was  thermally 
annealed  in  the  presence  of  a  magnetic  field.  The  sensor  was 
also  more  sensitive  with  longer  lengths  of  optical  fibers  bonded 
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to  the  sensor  arm  and  when  the  AC  magnetic  field  to  be  detected 
was  superimposed  on  a  DC  magnetic  field  flux  density.  Thinning 
the  fiber  and  thermally  annealing  the  nickel  rod  did  not  affect 
sensor  performance.  A  complete  summary  of  the  research  project 
is  included  in  the  next  Chapter. 
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VI.  Conclusion  and  Recommendations 


This  Chapter  summarizes  the  results  of  this  research  project 
and  provides  recommendations  for  further  study  of  the  sensor. 
Conclusion 

A  magnetostrictive  fiber  optic  interferometric  sensor  was 
fabricated  and  evaluated.  To  evaluate  the  sensor,  a  Helmholtz 
coil  was  designed,  fabricated,  and  calibrated.  A  measurement 
scheme  was  then  developed  to  determine  the  sensor's  sensitivity. 
Once  the  sensor's  response  could  be  quantified,  the  sensor's 
performance  was  evaluated  for  different  sensor  arms.  Finally, 
the  sensor's  effectiveness  was  optimized  using  a  DC  bias  magnetic 
field  flux  density. 

The  goal  of  this  research  project  was  to  investigate  a 
magnetic  field  sensor  that  has  the  potential  to  detect  a  pilot's 
biomagnetic  fields.  A  Helmholtz  coil  was  used  to  simulate  the 
biomagnetic  fields.  This  Helmholtz  coil  was  designed  to  generate 
a  uniform  magnetic  field  flux  density  whose  magnitude  could  be 
varied.  To  verify  the  Helmholtz  coil's  performance,  it  was 
calibrated  with  a  Hall  effect  magnetometer,  and  was  found  to 
generate  magnetic  field  flux  densities  to  within  fifteen  percent 
of  the  theoretically  predicted  values. 

Once  the  source  of  the  magnetic  field  was  characterized, 
three  different  interferometer  configurations  were  assembled,  and 
a  measurement  method  was  developed  and  validated.  After  each 
interferometer  configuration  was  fabricated,  its  output  optical 
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power  level  and  stability  were  examined.  The  fused  coupler, 
which  was  used  to  recombine  the  light  propagating  in  each  arm  of 
the  interferometer,  contributed  the  largest  source  of  noise.  A 
piezoelectric  ring  sensor  arm  was  used  to  verify  that,  as  the 
signal  strength  at  25  Hz  decreased,  the  amplitude  of  the 
frequency  component  at  25  Hz  of  the  sensor's  response  also 
decreased . 

Sensor  sensitivity  measurements  were  then  made  with  various 
sensor  arms.  No  DC  magnetic  field  flux  density  measurements 
could  be  made  due  to  the  noise  contributed  by  the  fused  coupler. 
However,  the  sensor  did  detect  a  25  Hz  AC  magnetic  field  using 
sensor  arms  composed  of  a  bulk  nickel  rod.  Sensors  consisting  of 
RF  sputtered  thin  films  of  cobalt,  iron,  and  nickel  on  the 
optical  fiber  did  not  detect  an  AC  magnetic  field  flux  density 
whose  magnitude  approached  12  G.  The  failure  of  the  thin  films 
was  attributed  to  fact  that  the  thin  film  could  not  generate 
sufficient  force  to  change  the  fiber's  length  a  detectable 
amount . 

Experiments  were  also  conducted  to  determine  the  effect  on 
the  sensor's  response  as  a  result  of  processing  the  nickel  rod 
and  altering  the  optical  fiber.  Thermally  annealing  the  nickel 
rod  at  758°  C  did  not  change  the  sensor's  performance.  However, 
as  theoretically  predicted,  thermally  annealing  the  nickel  rod  in 
the  presence  of  an  external  magnetic  field  (1000  G)  did  produce 
improved  sensor  sensitivity.  Thinning  the  optical  fiber  that  was 
bonded  to  the  nickel  rod  did  not  affect  the  sensor's  response 
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because  the  nickel  rod  produced  sufficient  force  to  sLr«=tch  the 
fiber  whether  the  fiber  was  full  size  or  reduced.  However,  a 
three  fold  increase  in  the  length  of  the  fiber  bonded  to  the 
nickel  rod  improved  the  sensor's  sensitivity. 

To  optimize  the  sensor's  performance,  a  DC  bias  magnetic 
field  flux  density  was  generated  in  addition  to  the  AC  magnetic 
field  flux  density  that  was  to  be  detected.  Increasing  the  DC 
bias  magnetic  field  flux  density  improved  the  sensor's  response 
for  a  particular  AC  magnetic  field  flux  density  level.  To 
determine  the  ultimate  sensitivity  of  the  sensor,  a  DC  bias 
magnetic  field  flux  density  of  0.4  G  was  applied  to  the  nickel 
rod  that  had  been  thermally  annealed  in  the  presence  of  an 
external  magnetic  field,  and  it  was  bonded  to  a  length  of  fiber 
that  was  three  times  the  nominal  length  of  fiber  (10  cm).  With 
this  optimal  conf iquration ,  the  sensor  detected  an  AC  magnetic 
field  of  10  mG  at  25  Hz. 

Recommendations 

This  research  effort  demonstrated  that  a  magnetostrict ive 
fiber  optic  interferometer  can  detect  low-level  magnetic  fields, 
and  it  also  provided  some  insight  into  the  behavior  of  the 
sensor.  However,  much  of  the  sensor's  characteristics  are  still 
unknown.  Additional  work  needs  to  be  performed  to  optimize  the 
sensor  so  that  it  may  one  day  achieve  its  theoretical 
sensitivity.  Specifically,  the  following  research  topics  should 
be  investigated: 

1.  The  magnetostrict ive  material  needs  to  be  more 
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thoroughly  understood.  A  model  needs  to  be  developed  that  can 
predict  the  amount  of  force  generated  by  a  magnetostrictive 
material.  This  model  should  include  the  following  material 
parameters:  type  (including  alloys),  structure  (crystalline, 

polycrystalline ,  and  amorphous),  thickness,  and  geometry. 

2.  The  fused  coupler,  which  provided  the  largest  noise 
source,  should  be  investigated.  A  model  should  be  developed  to 
explain  the  interaction  of  the  fields  in  the  device.  This  model 
could  then  be  used  to  vary  the  design  of  the  fused  coupler  to 
minimize  the  noise. 

3.  Another  approach  to  reduce  the  noise  from  the  fused 
coupler  involves  evaluating  polarization-preserving  optical 
fibers  and  fused  couplers. 

4.  Alternate  methods  of  depositing  metal  (such  as 
electroless  plating)  on  the  optical  fibers  should  be  evaluated. 
The  combination  of  two  or  more  metal  deposition  methods  may 
produce  adequate  quality  films  of  sufficient  thickness. 

5.  Different  cooling  procedures  for  the  fabrication  of 
thermally  annealed  thin  film  coated  optical  fibers  should  be 
explored  to  prevent  the  fibers  from  breaking.  A  microprocessor 
controlled  oven  should  be  used  to  carefully  control  the 
temperature  processing  profile  of  the  fiber. 

6.  Bulk  magnetostrictive  materials  should  also  be  furtner 
investigated.  Different  materials  (such  as  met-glass  alloys) 
should  be  evaluated.  The  method  of  bonding  the  fiber  to  the  bulk 
material  should  be  refined. 
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7.  The  method  of  determining  sensor  sensitivity  should  be 
reexamined.  The  use  of  different  signal  analysis  equipment,  such 
as  a  lock-in  amplifier,  may  provide  an  improved  signal-to-noise 
ratio . 

8.  To  improve  the  sensor's  sensitivity,  the  length  of  the 
fiber  in  contact  with  rhe  magnetostrictive  material  should  be 
increased. 

9.  To  examine  the  response  of  the  sensor  with  thin  film 
sensor  arms,  a  larger  magnetic  field  may  be  required.  To 
generate  a  larger  magnetic  field,  a  Helmholtz  coil  should  be 
fabricated  with  more  turns  of  wire. 

10.  The  noise  in  the  sensor  could  be  reduced  by  minimizing 
the  separation  between  the  sensor  arm  and  reference  arm.  One 
method  to  accomplish  this  that  should  be  investigated  is  to  use 
two  fused  couplers  and  bond  the  sensor  and  reference  arms 
together . 

11.  To  further  reduce  the  sensor's  noise,  an  improved  test 
chamber  should  be  constructed.  This  test  chamber  should  shield 
the  optical  fiber  components  from  the  ambient  magnetic  field, 
thermal  variations,  and  vibrations. 

12.  The  possibility  of  a  fabricating  a  single-mode  metallic 
glass  waveguide  should  be  investigated.  This  device  would  remove 
the  f iber-magnetostrict ive  material  bonding  problems  by  allowing 
the  waveguide  to  interact  directly  with  the  magnetic  field  that 
is  to  be  detected. 
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This  appendix  reviews  the  units  of  magnetism  that  are 
relevant  to  the  discussion  of  fiber  optic  magnetic  field  sensors 
(10:57).  For  this  thesis,  an  important  quantity  is  the  magnetic- 
field  flux  density,  (B)  which  is  measured  in  units  of  Gauss. 

This  figure  of  merit  is  directly  proportional  to  magnetic  field 
(H)  since  B  =  mH,  where  n  is  the  permeability  of  the  material. 
For  this  application,  the  permeability  of  the  medium  inside  the 
Helmholtz  coil  is  approximated  to  be  that  of  free  space 
(u=juojur=juo ( 1)  =47r  x  lO^Hm’1).  Table  A  describes  the  relationship 
between  these  two  quantities  in  the  various  system  of  units. 


Table  A.  Magnetic  Unit  Conversion 


System  of  Units 


Quantity 

cqs 

SI 

other 

Magnetic  Flux 
Density,  B 

1  Gauss(G) 

1  O'4  Tesla  (T) 

10  4  Weber/m2 

10"5  Gann 

Magnetic 
Field,  H 


1  Oersted(Oe) 


79.6  A/m 
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MATHCAD  Template  Used  for  Helmholtz  coil  off-axis  magnetic  flux 

density  calculation 

This  is  a  MATHCAD  template  which  was  used  to  calculate 

the  magnetic  flux  density  for  a  Helmholtz  coil  arrangement 

using  a  numerical  summation  technique  applied  to  the 
differential  form  of  the  Biot-Savart  law. 

An  explanation  of  the  Biot-Savart  law  and  a  figure  (Figure  III.  3) 

illustrating  the  geometry  of  the  problem  is  included  in  the  text 

of  this  thesis  (Chapter  III). 

In  the  text,  an  expression  for  the  magnetic  flux  density 

along  the  axis  between  the  two  coils  was  derived  in  closed 

form.  In  order  to  calculate  the  off-axis  magnetic  flux 

density,  a  numerical  summation  technique  was  used. 
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FILENAME:  FIELD 

THIS  PROGRAM  WAS  USED  TO  CALCULATE  THE  ELECTRIC  AND 
MAGNETIC  FIELDS  PROPAGATING  IN  THE  FIBER  OPTIC  CABLE  USED 
USED  IN  THIS  RESEARCH  PROJECT.  THE  FIELDS  ARE  CALCULATED 
FOR  A  FIXED  ANGLE  OF  0  WITH  VARIABLE  DISTANCE  FROM  THE 
OPTICAL  AXIS.  NUMERICAL  SERIES  APPROXIMATIONS  FOR  THE 
BESSEL  FUNCTIONS  AND  MODIFIED  BESSEL  FUNCTIONS  ARE  USED. 
THESE  NUMERICAL  APPROXIMATIONS  ARE  ONLY  VALID  WHEN  THEIR 
ARGUMENT  IS  LESS  THAN  OR  EQUAL  TO  3 . 

'k'k'k-k’k'k'k'klc'k-k'k-k-kic'k-k-klc-k-k'k'k'k'k'kic-k’k'k'k'k'k'k'k-k-kjc'k-k-k-k'k'k-k-k-k-k-k'k  k  -k  -k  -k  -k  -k 


170  ' 

180  '  CONSTAN 

190  ' 

1000  V  =  1.88 
1010  U=1 . 4  8 
1012  ' 

1020  W= ( V * 2-U *  2 )  " . 5 
1022  ' 

1030  A=1.75*10'-6 
1040  L= . 6328*10~-6 
1050  Nl=l .461 
1060  N2=1.457 
1062  EO=8 . 85E-12 
1064  El— N1 *  2  *EO 
1066  E2=N2  *  2 *EO 
1070  C=3  *  10  *  8 
1080  UO=1260*10*-9 
1090  PI=3. 14159 
1100  K=2*PI/L 
SOURCE 

1110  B=( (K*N1) *2-(U/A) “2) 

1120  F=2  *PI *C/L 

(RADS) 

1122  ' 

1124  '******************* 

1126  1 

1128  '  A , B , C , AND 

1130  ' 

1190  X=U : GOSUB  3000 
] 195  X=W : GOSUB  3036 


CONSTANTS  DEFINED 


'NORMALIZED  FREQUENCY 
'NORMALIZED  TRANSVERSE  PHASE 
CONSTANT  IN  THE  CORE 
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'CLADDING  INDEX  OF  REFRACTION 
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'PI 
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'PHASE  CONSTANT 
'FREQUENCY  OF  OPTICAL  SOURCF! 


A , B , C , AND  D  CONSTANTS  CALCULATION 

=U : GOSUB  3000  'JO, Jl, J2, Jl'  CALCULATED 
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1254  LPRINT  "  R 

HY  HZ" 

FOR  R  -  0  TO  1.1  STEP  .1 


EX 


EY 


CORE 


EZ 


HX 

1255 

AXIS 

1257 

1260 

1262 

1265 
AXIS 

1266 
1270 
1280 
1282 
1284 
1290 
1300 
1304 
1310 
1320 

1322 

1323 

1324 

1325 

1326 

1327 

1328 


'SCALED  DISTANCE  FROM  OPTICAL 

'OFFSET  FOR  SINGULARITY 
'ARGUMENT  FOR  J  AND  K 
' JO, Jl, J2 , J1 '  CALCULATED 
'ACTUAL  DISTANCE  FROM  OPTICAI 

' ANGLE 


R  =  R  +.0000001 
X  =  U*R 
GOSUB  3000 
RS  =  R*A 

T=T*PI/180 

ER  =  (-B*JlP/(U/A)+BB*F*UO*Jl/( (U/A) ~2*RS) ) *SIN(T) 
ET  =  (-3*Jl/( (U/A) ‘2*RS)+BB*F*U0*J1P/ (U/A) ) *COS (T) 
HR  =  (F*E1*J1/ ( (U/A) * 2*RS) -BB*B*J1P/ (U/A) ) *COS (T) 
HT  =  ( -F*E1*J1P/ (U/A) +BB*B*J1/ ( (U/A) *  2*RS) ) *SIN (T) 
EX  =  ER*COS (T) -ET*SIN (T) 

EY  =  ER*SIN (T) +ET*COS (T) 

EZ=J1*SIN(T) 

HX  =  HR*COS (T) -HT*SIN (T) 

HY  =  HR*SIN (T) +HT*COS (T) 

HZ=BB* Jl*COS (T) 


********************************************************* 
CALCULATION  FOR  FIELDS  IN  CLADDING 


1330 

1332 

1334 


" ; R ; EX , EY , EZ , HX , HY , HZ 


LPRINT  TAB (10)  ; 

LPRINT  USING  "  +  ##.###'## 

NEXT  R 

LPRINT  "  ************************ 

***************************" 

1336  LPRINT  " 

HY  HZ" 

=  1  TO  3.1  STEP  . 1 

+  .0000001 


CLADDING 


R 


EX 


EY 


HX 

1338 

1339 

1350 

1351 

1352 
1360 
1370 


FOR  R 
R  =  R 
X=W*R 
GOSUB 
GOSUB 
ERC  - 
ETC  = 


3000 

3036 

(CC*B*K1P/ (W/A) -DD*F*U0*K1/ ( (W/A) ‘2*RS) ) *SIN(T) 
( CC*B*K1/ ( (W/A)  "  2  *RS ) -DD*F*U0*K1P/ (W/A)  ) *C0S (T) 


C  — 2 
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2994  '  SUBROUTINES  THAT  CALCULATE  JO , J1 , J2 , J1 ' , KO , K1 , K2 , K1 1 

2995  ' 

3  0  0  0 

J  0  =  1  —  2 . 2499997#* (X/3) ‘2  +  1. 2656208#* (X/3)  *4-. 3163866* (X/3)  "6+ . 0444 
479*(X/3)'8-.0039444*(X/3)‘10+.00021*(X/3)*12 

3  0  10 

J1=X* ( . 5-. 56249985#* (X/3)  *  2+ . 2 109  3 57 3 # * (X/3)  *  4-3 . 9 54 2 89E-02 * ( X/ 3 ) 
*6  +  4 . 43319E-03* (X/3)  *  8-3 . 17 61E-04 * (X/3 )  *  10+1 . 109E-05* ( X/3 )  *12) 
3020  J2=2*J1/X-J0 
3030  J1P=. 5* (J0-J2) 

3032  'RETURN 

3036  IF  X>2  THEN  4000 

3  0  3  7 

10  =  1+3.51562  29#* (X/3 .75)  *2  +  3 . 08994  24#* (X/3 .75)  *  4  +  1 . 2 067492 # * ( X/ 3 . 
75)  *6+. 2659732* (X/3. 75)  *8+ . 03  607  68* (X/3 . 75)  *  10+ . 004 58 13  * ( X/ 3 . 7 5 )  * 
12 

3  0  3  8 

1 1=. 5+. 87  89  0594#* (X/3 . 75)  *  2+ . 514  988  69#* (X/3 . 7 5 )  *  4+ . 15  08  4 9 3 4 # * ( X/ 3 
.75)  *6+2 . 658733E-02* (X/3 . 75)  *  8  +  3 . 01532E-03* (X/3 . 75)  *  10  +  3 . 24 1 1 1E-0 

4  * ( X/ 3.75)  *12 
3039  I1=I1*X 

3  0  4  0 

KO=-LOG(X/2) *10-. 577  2 1566#+. 4  2 2  784  2  * (X/2)  " 2+ . 2 3 0697  56 # * ( X/2 )  ‘4+. 0 
3  4  8859* (X/2)  *  6+2 . 62 69 8E-0 3  * ( X/2 )  *  8+ . 000107 5* ( X/2 )  *  10+ . 0000074  * ( X/ 
2)  *  12 

3  0  4  5 

Kl=X*LOG (X/2 ) *I1+1+. 15443144#* (X/2) *2-. 67278579#* (X/2) *4-. 181 5680 
7  #  * ( X/2 )  *6-1.919402E-02* (X/2)  *  8 - 1 .  104 04E-03 * (X/2 )  ‘10 

3050  K1=X * (-1) *K1 
3060  K2=2*K1/X+K0 
3070  K1P=- . 5* ( K0+K2 ) 

3080  IF  X<=2  THEN  4070 

4  0  0  0 

K0=1 . 25331414#-. 07832* (2/X) +2 . 189568E-02* (2/X) *3-1 . 062446E-02* (X 
2) *3+5 . 87872E-03* (2/X) *4 

4 C 10  K0  =  X*  (-.5) *EXP(-X) *K0 

4  0  2  0 


C-3 


Kl=l. 253 3 14 14 #+.2 3498619#* (2/X) -.0365562* (2/X)  2+1. 5U4268E-02* (2/ 

X)  ~  3-7  . 8  03  53E-03* (2/X)  '4 

4030  K1=X~ (-.5) *EXP(-X) *K1 

4040  K2=2*K1/X+K0 

4050  K1P=- . 5* ( K0+K2 ) 

4070  RETURN 
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********************************************************* 


FILENAME:  UCALC 

THIS  PROGRAM  WAS  USED  TO  CALCULATE  THE  NORMALIZED 
TRANSVERSE  PHASE  CONSTANT  IN  THE  CORE  (U)  IN  THE 
OPTICAL  FIBER  USED  IN  THIS  RESEARCH  PROJECT. 

THE  PROGRAM  USES  NUMERICAL  SERIES  APPROXIMATIONS 
FOR  THE  BESSEL  FUNCTIONS  AND  MODIFIED 

BESSEL  FUNCTIONS.  THESE  NUMERICAL  APPROXIMATIONS  ARE  ONLY 
VALID  WHEN  THEIR  ARGUMENT  IS  LESS  THAN  OR  EQUAL  TO  3. 

TO  SOLVE  FOR  U  REQUIRES  SOLVING  A  TRANSCENDENTAL  EQUATION. 
THE  EQUATION  WAS  SOLVED  BY  SOLVING  EACH  SIDE  FOR  VARIOUS 
VALUES  OF  U  UNTIL  THE  DIFFERENCE  WAS  MINIMIZED. 

********************************************************* 


1000  A  =  0  'START  OF  U  VALUE  RANGE 

1010  B  =  1.8  'END  OF  U  VALUE  RANGE 

1020  S  =  .1  'STEP  SIZE  IN  U  VALUE  RANGE 

1022  LPRINT  "  " 

1023  LPRINT  "  " 

1025  OLDDIFF  =  1E+30  'INITIAL  DIFFERENCE 


1027  V  =  1.87 

1028  LPRINT  TAB ( 10 )  ; 

1029  LPRINT  "  U  DIFFERENCE" 

1030  FOR  U  =  A  TO  B  STEP  S 
1035 

1040  X  =  U  +  .0000001 
1050  GOSUB  3000 
1060  LEFT  =  J0/(X*J1) 

1070  X  =  (V*  2-U~  2 )  * . 5 

1080  GOSUB  3000 

1090  RIGHT  =  K0/(X*K1) 

1100  NEWDIFF  =  ABS (LEFT-RIGHT) 

1102  LPRINT  TAB (10) ; 

1105  LPRINT  USING  "#.###  ";U; 

1107  LPRINT  USING  "##.## - " ; NEWDIFF 

1110  IF  NEWDIFF>OLDDIFF  THEN  1200 
1112 

1115  OLDDIFF=NEWDIFF 


'NORMALIZED  FREQUENCY 


'LOOP  USED  TO  CALCULATE 
'EACH  SIDE  FOR  A  DIFFERENCE 
'OFFSET  FOR  SINGULARITY 
'CALCULATE  J0,J1,K0,K1 
'LEFT  SIDE  OF  EQUATION 
'SCALED  VARIABLE 
'CALCULATE  J0,J1,K0,K1 
'RIGHT  SIDE  OF  EQUATION 
'LEFT  -  RIGHT 


'CONTINUE  LOOP  UNTIL 
'LEFT  -  RIGHT  MINIMIZED 
'RESET  DIFFERENCE 


1120  NEXT  U 

1200  IF  S  <=  .001  THEN  1300  'CHECK  STEP  SIZE 


C-4 


1210  A  =  U-2*S  'NEW  INITIAL  U  RANGE 

1220  B  =  U  'NEW  FINAL  U  RANGE 

1230  S  =  S/10  'NEW  STEP  SIZE 

1235  OLDDIFF  =  OLDDIFF  +  1 

1240  GOTO  1030 

1300  END 

2990  ' 

2992  '********************************************************* 

2993  ' 

2994  ' 

2995  ' 

2996  '  SUBROUTINES  USED  TO  CALCULATE  J0,J1,K0,K1 
2998  ' 

3  0  0  0 

J0=l-2 .2499997#* (X/3 ) * 2+1 . 2656208#* (X/3 ) * 4- . 3 163  866* ( X/3 )  *  6+ . 04 4 4 
479* (X/3) '8-. 0039  44  4* (X/3)  *  10+ . 0002 1* (X/3 ) "12 

3  0  10 

J1=X* ( .5-. 56249985#* (X/3) *  2+ . 2 109 3 57 3 #* (X/3) *4-3 . 954289E-02* (X/3 ) 
*  6+4 . 43319E-03* (X/3 )  *  8-3 . 17 6 IE-04 * (X/ 3 )  *  10+1 . 109E-05* (X/3)  *12) 
3020  J2  =  2 * Jl/X-J 0 
3030  J1P=. 5* (J0-J2) 

3032  'RETURN 

3036  IF  X>2  THEN  4000 

3  0  3  7 

10  =  1+3 . 5156229#* (X/3 .75)  *2  +  3 . 0899424#* (X/3 . 75)  *  4  +  1 . 2 0674 92 # * ( X/ 3 . 
75)  *6+. 26597 32* (X/3. 75)  *8+ . 03  60768* (X/3 . 75)  *  10+ . 004 58 13 * ( X/ 3 . 7 5 )  ' 
12 

3  0  3  8 

1 1=. 5+. 87890594#* (X/3 .75)  *  2+ . 51498  8  69 #* (X/3 . 7  5 )  " 4+ . 15084934 #*(  X/ 3 
.75)  '6+2.658733E-02* (X/3 .75)  *8  +  3 . 01532E-03* (X/3 . 75)  *  10+3 . 24 1 1 1 E-0 
4* (X/3 .75) *12 
3039  I1=I1*X 

3  0  4  0 

KO=-LOG ( X/2 ) *10-. 57721566#+. 4227842* (X/2) *2+. 23069756#* (X/ 2) '4+ . 0 
348859* (X/2)  *  6  +  2 . 62 6O8E-03 * ( X/2 )  * 8+ . 0001075* (X/2 )  *  10+ . 0000074 *( X/ 
2)  “12 

3  0  4  5 

Kl=X*LOG (X/2 ) *11+1+ .  1544314  4#* (X/2)  *  2- . 67278579 #* (X/2 )  *4-. 181568  9 
7  #  * ( X/2 )  *6-1.919402E-02*(X/2)  *  8-1 . 104 04E-03 * (X/2 )  *10 

3050  K1=X* (-1) *K1 
3060  K2=2*K1/X+K0 
3070  K1P=- . 5* (K0+K2 ) 

3080  IF  X<=2  THEN  4070 

4  0  0  0 

K0=1. 25331414#-. 07832* (2/X) +2 . 189568E-02* (2/X) ‘3-1 . 062446E-02* ( X/ 
2) *3+5.87872E-03*(2/X) *4 

4010  K0=X* (-. 5) *EXP(-X) *K0 

4  0  2  0 

Kl=l . 25331414#+. 23498619#* (2/X) -. 0365562* (2/X) " 2+1 . 5042 68E-02 * ( 2/ 
X) *3-7.80353E-03*(2/X) ‘4 
4030  K1=X * (-.5) *EXP( -X) *K1 
4040  K2=2*K1/X+K0 


C-5 


4050  K1P=- . 5* (K0+K2 ) 
4070  RETURN 


C-6 


*  *  * 

+  0  . 
+  0  . 
+  0  . 
+  0  . 
+  0. 
+  0  . 
+  0. 
+  0  . 
+  0  . 
+  0  . 
+  1  • 
★  -k  * 

+1. 
+1. 
+1. 
+  1 . 
+ 1  . 
+ 1  . 
+ 1  . 
+  1 . 
+ 1  . 
+ 1  . 
+  2  . 
+  2  . 
+  2  . 
+  2  . 
+  2  . 
+  2  . 
+  2  . 
+2  . 
+  2  . 
+  2  . 
+  3  . 


****************** 

CORE  *********************** 

R 

EX 

EY 

EZ 

HX 

HY 

0000 

+0 . UwOO 

-17 . 1938 

+0. 0000 

+  0 . 0666 

+0 . 0000 

-0 

1000 

+0.0000 

-17.0996 

+0.0000 

+0.0662 

+0.0000 

-0 

2000 

+0.0000 

-16.8185 

+0.0000 

+0.0651 

+0. 0000 

-0 

3000 

+0.0000 

-16.3551 

+0.0000 

+0.0633 

+0 . 0000 

-0 

4000 

+0.0000 

-15.7170 

+  0 . 0000 

+0.0609 

+0.0000 

-0 

5000 

+0.0000 

-14 .9147 

+0.0000 

+0.0578 

+  0 . 0000 

-0 

6000 

+0.0000 

-13 .9613 

+0.0000 

+0.0541 

+  0 . 0000 

-0 

7000 

+0.0000 

-12 .8721 

+  0 . 0000 

+0.0499 

+0.0000 

-0 

8000 

+0.0000 

-11.6649 

+0.0000 

+0.0452 

+0.0000 

-0 

9000 

+0.0000 

-10.3592 

+0 . 0000 

+0 . 0402 

+  0 . 0000 

-0 

0000 

+0.0000 

-8 . 9757 

+0.0000 

+0.0348 

+0. 0000 

-0 

************* 

***** 

CLADDING 

******************** 

R 

EX 

EY 

EZ 

HX 

HY 

0000 

+0.0000 

-8 . 9758 

+0.0000 

+0. 0346 

+0. 0000 

-0 

1000 

+0.0000 

-6.8640 

+0.0000 

+0 . 0264 

+  0 . 0000 

-0 

2000 

+0.0000 

-5.3338 

+  0 . 0000 

+0.0205 

+  0 . 0000 

-0 

3000 

+0.0000 

-4 . 2002 

+0.0000 

+0.0162 

+0. oooo 

-0 

4000 

+0.0000 

-3 . 3450 

+0.0000 

+0. 0129 

+0 . 0000 

-0 

5000 

+  0 . 0000 

-2 . 6895 

+0.0000 

+0.0104 

+0. 0000 

-0 

6000 

+  0 . 0000 

-2 . 1804 

+0 . 0000 

+0.0084 

+  0 . 0000 

-0 

7000 

+0. 0000 

-1.7803 

+0. 0000 

+0.0069 

+  0 . 0000 

-0 

8000 

+0.0000 

-1.4605 

+0.0000 

+0.0056 

+  0 . 0000 

-0 

9000 

+0.0000 

-1.1972 

+0.0000 

+0.0046 

+  0 . 0000 

-0 

0000 

+0.0000 

-0.9861 

+0.0000 

+0.0038 

+0 . 0000 

-0 

1000 

+0.0000 

-0.8156 

+0.0000 

+0.0031 

+  0 . 0000 

-0 

2000 

+  0.0000 

-0.6769 

+0.0000 

+  0 . 0026 

+  0 . 0000 

-0 

3000 

+0.0000 

-0.5635 

+0.0000 

+0.0022 

+  0 . 0000 

-0 

4000 

+0.0000 

-0.4702 

+0.0000 

+0.0018 

+0 . 0000 

-0 

5000 

+0. 0000 

-0.3932 

+0.0000 

+0.0015 

+0 . 0000 

-0 

6000 

+0.0000 

-0.3293 

+0.0000 

+0.0013 

+0.0000 

-0 

7000 

+0.0000 

-0.2762 

+0.0000 

+0. 0011 

+0. 0000 

-0 

8000 

+0.0000 

-0 .2318 

+0.0000 

+0 . 0009 

+0 . 0000 

-0 

9000 

+0.0000 

-0 . 1947 

+0.0000 

+0. 0007 

+0 . 0000 

-0 

0000 

+0 . 0000 

-0. 1635 

+0. 0000 

+0.0006 

+  0 . 0000 

-0 

*  *  * 

HZ 

0000 
0003 
0006 
0008 
0011 
0013 
0016 
0018 
0019 
0021 
0022 
*  *  * 

HZ 

0022 

0018 

0015 

0013 

0011 

0009 

0008 

0007 

0006 

0005 

0004 

0004 

0003 

0003 

0002 

0002 

0002 

0002 

0001 

0001 

0001 
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u 

DIFFERENCE 

0.000 

2  .  OOE-f-14 

0.100 

1.99E+02 

0 .200 

4.93E+01 

0.300 

2. 15E+01 

0.400 

1 . 18E+01 

0 . 500 

7 . 30E+00 

0.600 

4 . 85E+00 

0 .700 

3 . 37E+00 

0.800 

2 . 40E+00 

0.900 

1.73E+00 

1.000 

1.25E+00 

1.100 

8.79E-01 

1.200 

5.90E-01 

1.300 

3 . 54E-01 

1.400 

1. 51E-01 

1.500 

3 . 32E-02 

1.600 

2 . 17E-01 

1.400 

1. 51E-01 

1.410 

1 . 32E-01 

1 .420 

1 . 13E-01 

1.430 

9 . 48E-02 

1 .440 

7 . 62E-02 

1 .450 

5.78E-02 

1.460 

3 . 94E-02 

1 .470 

2 . 12E-02 

1.480 

3 . 03E-03 

1.490 

1 . 51E-02 

1 .470 

2 . 12E-02 

1.471 

1.94E-02 

1.472 

1.76E-02 

1.473 

1 . 57E-02 

1 .474 

1 . 39E-02 

1.475 

1 . 2 IE-02 

1 .476 

1. 03E-02 

1 .477 

8.48E-03 

1 . 478 

6.66E-03 

1.479 

4 . 85E-03 

1.480 

3 . 03E-03 

1.481 

1.22E-03 

1 . 482 

5 . 96E-04 

1.483 

2 . 41E-03 
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1000 

1010 

1020 

1100 

1110 

1120 

1130 

1170 

1180 

1190 

1210 

1240 

1250 

1260 

1270 

1280 

1290 

1490 

1492 

1493 

1494 
1496 
1498 
1500 
1510 
1520 
1530 
1540 


FILENAME:  BKDISP 

THIS  PROGRAM  WAS  USED  TO  TAKE  DATA  FROM  THE  BREUL  &  KJAER 
TYPE  2032  SIGNAL  ANALYZER  VIA  A  METROBYTE  IEEE  488 
CARD  AND  A  ZENITH  248  COMPUTER.  THE  PROGRAM  USES 
SUBROUTINES  TO  SEND  COMMANDS  TO  AND  RECEIVE  DATA  FROM  THE 
B&K  2032.  THE  PROGRAM  INITIALIZES  BOTH  THE  IEEE  CARD  AMD 
THE  B&K  2032.  AFTER  INITIALIZATION,  THE  PROGRAM  PRESENTS 
A  MENU  WITH  CHOICES  OF  EXAMINING  OR  RECORDING  THE  TIME  OR 
SPECTRAL  INFORMATION  IN  THE  EXCITATION  OR  RESPONSE. 


INITIALIZATION 


OPTION  BASE  1 
GOSUB  9000 
GOSUB  9100 
0$  =  "ED  DS , 0 " 
GOSUB  9300 
0$  =  "EM  MM, 1" 
GOSUB  9300 
0$  =  "EM  AF , 0 " 
GOSUB  9300 
0$  =  "EM  AU , 0 " 
GOSUB  9300 
0$  =  "EM  BF, 0" 
GOSUB  9300 
0$  =  "EM  BU , 0" 
GOSUB  9300 
0$  =  "EM  MC , 2 " 
GOSUB  9300 


'ARRAYS  BEGIN  WITH  ELEMENT  1 
'INITIALIZE  THE  INTERFACE 
'  REMOTE 

'INPUT  DATA  SOURCE 
'SPECTRUM  AVERGING  ZERO  PAD 
'CH  A  FILTER  OFF 
'VOLT  UNITS  FOR  CH  A 
'CH  B  FILTER  OFF 
'VOLT  UNITS  FOR  CH  B 
'DUAL  MEASUREMENT  CHANNEL 


THE  MENU  IS  DISPLAYED 


CLS 

PRINT  "ENTER  CHOICE: 
PRINT  "  " 

PRINT  " 

PRINT  "  " 


1)  DISPLAY  DUAL  TIME" 

2)  DISPLAY  RESPONSE  TIME  ONLY" 


1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1710 

1720 

1722 

1723 

1724 

1725 

1726 
1730 
1740 

1990 

1991 

1992 

1993 

1994 
2000 
2010 
2020 
2030 
2040 
2050 
2052 
2054 
2080 
2090 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2  2  50 
2300 
2310 
2312 
2314 
2480 
2490 


PRINT  " 
PRINT  " 

II 

3) 

DISPLAY  EXCITATION  TIME  ONLY" 

PRINT  " 
PRINT  " 

II 

4) 

DISPLAY  DUAL  SPECTRUM  " 

PRINT  " 
PRINT  " 

II 

5) 

DISPLAY  RESPONSE  SPECTRUM  ONLY" 

PRINT  " 
PRINT  " 

It 

6) 

DISPLAY  EXCITATION  SPECTRUM  ONLY 

PRINT  " 
PRINT  " 

II 

7) 

STORE  TIME  FILE" 

PRINT  " 
PRINT  " 

If 

8) 

STORE  SPECTRUM  FILE" 

INPUT  " 

9) 

END" ;R 

GOSUB  9100 

' REMOTE 

i 


'•a********************************************************* 

i 

'  BASED  ON  MENU  CHOICE,  PROGRAM  BRANCHES  TO  SUBROUTINES 

I 

ON  R  GOSUB  2000, 2500, 3000, 3500, 4000, 4500, 5000, 5500, 7000 
GOTO  1500 

i 

•  ********  -;************************************************* 

I 

'  DUAL  TIME  DISPLAY  COMMANDS 


0$  =  "DF  UM" 

'UPPER  DISPLAY 

GOSUB  9300 

0$  =  "ED  FU , 1 " 

'TIME  CHANNEL  B 

GOSUB  9300 

0$  =  "ED  FC , 0" 

'REAL  FUNCTION  COORD 

GOSUB  9300 

0$  =  "ED , CT, 0" 

'MAIN  CURSOR 

GOSUB  9300 

0$  --  "ED  CU,0" 

'CURSOR  POSITION  =  0 

GOSUB  9300 

0$  =  "EM  TM, 1" 

'TRIGGER  ON  CHANNEL  A 

GOSUB  9300 

0$  =  "EM  FS, 1600" 

'FREQUENCY  SPAN  =  1 . 6 KHZ 

GOSUB  9300 

0$  =  "DF  LM" 

'LOWER  DISPLAY 

GOSUB  9300 

0$  =  "ED  FU , 0" 

'TIME  CHANNEL  A 

GOSUB  9300 

0$  =  "ED, CT, 0" 

'MAIN  CURSOR 

GOSUB  9300 

0$  =  "ED  CU,0" 

'CURSOR  POSITION  =  0 

GOSUB  9300 

0$  =  "DF  UL" 

'DISPLAY  UPPER/LOWER  GRAPHS 

GOSUB  9300 

GOSUB  9200 

' LOCAL 

RETURN 

D-2 


2491 

2492 

2493 

2494  '  RESPONSE  TIME  DISPLAY  COMMANDS 

2495  ' 

2500  0$  =  "DF  UM" 

2510  GOSUB  9300 
2520  0$  =  "ED  FU, 1" 

2530  GOSUB  9300 
2532  0$  =  "ED , CT , 0 " 

2534  GOSUB  9300 
2540  0$  =  "ED  CU,0" 

2550  GOSUB  9300 
2560  0$  =  "EM  TM, 1" 

2570  GOSUB  9300 
2580  0$  =  "EM  FS, 1600" 

2590  GOSUB  9300 
2592  0$  =  "DF  UF" 

2594  GOSUB  9300 
2600  GOSUB  9200 
2610  RETURN 

2901  ' 

2902  '*****•***************************************************** 

2903  ' 

2904  '  EXCITATION  TIME  DISPLAY  COMMANDS 

2905  ' 

3000  0$  =  "DF  UM" 

3010  GOSUB  9300 
3020  0$  =  "ED  FU , 0" 

3030  GOSUB  9300 
3032  0$  -  "ED, CT, 0" 

3034  GOSUB  9300 
3040  0$  =  "EM  TM , 1 " 

3050  GOSUB  9300 
3060  0$  =  "EM  FS , 1600" 

3070  GOSUB  9300 
3072  0$  =  "DF  UF" 

3074  GOSUB  9300 
3080  GOSUB  9200 
3090  RETURN 

3490  ' 

3491  *  ★  •*•**★•*★★★★★★•**•*•*★★★★**•**■**★*★■***•*****■*  a  h***************** 

3492  * 

3493  ’  DUAL  SPECTRUM  DISPLAY  COMMANDS 

3494  ' 

3500  0$  =  "DF  UM"  'UPPER  DISPLAY 

3510  GOSUB  9300 

3520  0$  =  "ED  FU, 11"  'SPECTRUM  CHANNEL  B 

3521  GOSUB  9300 

3522  0$  =  " ED , CT , 1 "  'HARMONIC  CURSOR 

3530  GOSUB  9300 

3540  0$  =  "EM, TM, 0"  'TRIGGER  SET  TO  FREE  RUN 


'UPPER  DISPLAY 

'TIME  CHANNEL  A 

'MAIN  CURSOR 

'TRIGGER  ON  CHANNEL  A 

'FREQUENCY  SPAN  =  1 . 6KHZ 

'UPPER  FULL  DISPLAY 

'  LOCAL 


'UPPER  DISPLAY 
'TIME  CHANNEL  B 
'MAIN  CURSOR 
'CURSOR  POSITION  =  0 
'TRIGGER  ON  CHANNEL  A 
'FREQUENCY  SPAN  =  1 . 6KHZ 
'UPPER  FULL  DISPLAY 
'  LOCAL 
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3550 

GOSUB 

9300 

3560 

0$  -  " 

EM  FS, 100"  ' 

3570 

GOSUB 

9300 

3572 

0$  =  " 

ED  CU,25"  ’ 

3574 

GOSUB 

9300 

3580 

0$  =  " 

DF  LM"  ’ 

3590 

GOSUB 

9300 

3600 

0$  =  " 

ED  FU, 10"  ' 

3602 

GOSUB 

9300 

3604 

0$  -  " 

ED  CT , 1 "  1 

3606 

GOSUB 

9300 

3612 

0$  =  " 

ED  CU,25"  ’ 

3614 

GOSUB 

9300 

3620 

0$  =  " 

EM  FS , 100"  ’ 

3630 

GOSUB 

9300 

3640 

0$  =  " 

DF  UL"  1 

3650 

GOSUB 

9300 

3660 

GOSUB 

9200  ’ 

3670 

RETURN 

3990 

I 

3991 

i ************************** 

3992 

1 

3993 

! 

RESPONSE  SPECTRUM 

3994 

f 

4000 

0$  =  " 

DF  UM"  • 

4010 

GOSUB 

9300 

4020 

0$  =  " 

ED  FU ,11"  ' 

4030 

GOSUB 

9300 

4032 

0$  =  " 

ED  CU ,25" 

4034 

GOSUB 

9300 

4040 

0$  =  " 

EM  FS, 100" 

4050 

GOSUB 

9  3  OC 

4052 

0$  =  " 

EM , TM , 0 "  ' 

4054 

GOSUB 

9300 

4060 

0$  =  " 

DF  UF" 

4070 

GOSUB 

9300 

4  08  0 

GOSUB 

9200 

4090 

RETURN 

4  4  9  0 

f 

4491 

1  ******************** c****** 

4  4  92 

1 

4  4  9  3 

1 

EXCITATION  SPECTRUM 

4  4  9  4 

f 

4  5  0  0 

0$  -  " 

DF  UM"  ' 

4  "->  1  0 

GOSUB 

9  3  0  0 

4  5  1  2 

0$  »  " 

ED  CU , 75"  ' 

4  5  1  4 

GOSUB 

9  3  9  0 

4  5  2  0 

0$  -  " 

E T  FU ,10"  ' 

4  5  3  0 

GOSUB 

9  3  00 

4  5  3  2 

0$  =-  " 

EM , TM , 0 " 

*  r  i  a 

t  >  >  'i 

GOSUB 

9  3  00 

4  5  4  0 

0$  -  " 

EM  FS , 3 00"  ' 

FREQUENCY  SPAN  =  100  HZ 
CURSOR  SET  TO  25  HZ 
LOWER  DISPLAY 
SPECTRUM  CHANNEL  A 
HARMONIC  CURSOR 
CURSOR  SET  TO  25  HZ 
FREQUENCY  SPAN  =  100  HZ 
UPPER/ LOWER  GRAPHS 
LOCAL 

******************************** 

DISPLAY  COMMANDS 
UPPER  DISPLAY 
SPECTRUM  CHANNEL  B 
CURSOR  SET  TO  25  HZ 
FREQUENCY  SPAN  =  100  HZ 
TRIGGER  SET  TO  FREE  RUN 
UPPER  FULL  GRAPH 
LOCAL 

■k-k'k'k-k'k'k’k'k’k'k’k'kic'k’k'k'k'k'k'k’k'k'k'k'k'k’k’k'k'k'k 

DISPLAY  COMMANDS 
UPPER  DISPLAY 
CURSOR  SET  TO  25  HZ 
SPECTRUM  CHANNEL  A 
TRIGGER  SET  TO  FREE  RUN 
FREQUENCY  SPAN  -  100  HZ 


D-  4 


' LOCAL 


TIME  DISPLAY  STORAGE  COMMANDS 


4550  GOSUB  9300 

4560  0$  =  "DF  UF"  'UPPER  FULL  GRAPH 

4570  GOSUB  9300 

4580  GOSUB  9200  'LOCAL 

4590  RETURN 

4990  ' 

4991  1  ★*★***★★★*★★★★★***★★*★★★**★***★★**★■*★★*★■*:★★***★★********** 

4992  ' 

4993  '  TIME  DISPLAY  STORAGE  COMMANDS 

4 99  ■.  ' 

5000  INPUT  "ENTER  DATA  FILE  NAME" ; DF$ 

5010  OPEN  "0" , # 1 , DF$ 

5020  FOR  I  =  1  TO  500  '  500  POINTS  OF  DATA 

5030  0$  =  "AF  " +STR$ ( I ) 

5040  GOSUB  9300 
5050  GOSUB  9400 

5060  DAT  =  VAL(MID$ (E$ , 2 , 4) ) *10'VAL(MID$ (E$, 8 , 3) )  'DAT/ 

CONVERSION 

5080  WRITE  #1,1, DAT 

5090  NEXT  I 

5400  CLOSE  #1 

54x0  GOSUB  9200 

5440  RETURN 

5490  ' 

5491  *********************************************************** 

5492  ' 

5493  '  SPECTRUM  DISPLAY  STORAGE  COMMANDS 

5494  ' 

5500  INPUT  "ENTER  DATA  FILE  NAME" ; DF$ 

5510  OPEN  "O" , # 1 , DF$ 

5530  FOR  I  =  1  TO  800  '  800  POINTS  OF  DATA 

5540  0$  =  "AF  " +STR$ ( I ) 

5550  GOSUB  9300 
5560  GOSUB  9400 

5570  DAT  =  VAL(MID$ (E$ , 2 , 6) )  'DATA  CONVERSION 

5580  WRITE  #1,DAT 

5590  NEXT  I 

5600  CLOSE  #1 

5605  GOSUB  9200 

5610  RETURN 

G990  ' 

6991  '  ****************************************************  *****  * 

6992  ' 

6993  '  END  OF  PROGRAM 

6994  ' 

7000  GOSUB  9200 
7010  END 
8983  ' 

3939  '******************************************************^*** 

3  9  9  0  ' 

8992  '  INITIALIZING  THE  IEEE  INTERFACE 

3  9  9  4  ' 


SPECTRUM  DISPLAY  STORAGE  COMMANDS 


DATA  CONVERSION 


INITIALIZING  THE  IEEE  INTERFACE 
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9000 

DEF  SEG  = 

9010 

CINITS  =  " 

9020 

CREMOTE$  = 

9030 

CLOCALS  = 

9040 

COUTPUTS  = 

9050 

CENTERS  = 

9060 

IE488  =  0: 

90/0 

CALL  IE488 

9080 

IF  FLG%<>0 

9090 

RETURN 

9095 

i 

9 

'  **** 

‘k-k-k'k'k'k-k'k*-k 

9097 

1 

9098 

'  SET  THE 

9099 

1 

9100 

CALL  IE488 

9110 

IF  FLG%<>0 

9120 

RETURN 

9195 

1 

9196 

*  -k-k'k-k-k-k-k-k'k 

9197 

l 

9198 

'  SET  THE 

9199 

i 

9200 

CALL  IE488 

9210 

IF  FLG%<>0 

9220 

RETURN 

9295 

i 

9296 

*  ★★*-****★★ 

9297 

I 

9298 

'  OUTPUT  A 

9299 

1 

&HDOOO 


0 :  FLG% 


0:  BRD%  =  0 


9395 

9396 

9397 

9398 

9399 
q 


9300  OUTPUTS  =  0$+CHR$ ( 10 ) 

9310  CALL  IE483  (COUTPUT$,  OUTPUTS,  FLG% ,  BRD% ) 

9320  IF  FLG%<>0  THEN  PRINT  "OUTPUT  ERROR:  " ; HEX$ ( FLG% ) : END 
9330  RETURN 


INPUT  DATA  FROM  THE  B&K 


E$  -  G PACES (15) 

9410  CALL  IE488  (CENTERS,  E$,  FLG% ,  BRD%) 

9420  IF  FLG%<>0  THEN  PRINT  "ENTRY  ERROR:  " ; HEX$ ( FLG% ) : END 
9430  RETURN 

•k-k'k-k^-k-k-k'k-k'k-k'k-k'k-k^'k-k-k-k'k-k'k-k-k-k'k-k-k'k'kii-k'k'k'k'k'k'k-k-k'k'k'k'k-^-k'k^-k-k-k-k'k'k-k-k-k^-k-kicr: 
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10  *********************************************************** 

20  ' 

30  ' 

40  '  FILENAME:  BK 

50  ' 

60  '  THIS  PROGRAM  WAS  USED  TO  TAKE  DATA  FROM  THE  BREUL  &  KJAER 

70  1  TYPE  2032  SIGNAL  ANALYZER  VIA  A  METROBYTE  IEEE  438 

INTERFACE 

80  '  CARD  AND  A  ZENITH  248  COMPUTER.  THE  PROGRAM  INITIALIZES 

90  '  THE  INTERFACE  CARD,  BUT  DOES  INITIALIZE  THE  B&K  2032.  THF 

100  '  PROGRAM  ASSUMES  THAT  THE  DATA  TO  BE  STORED  IS  ALREADY 

110  '  ON  THE  SCREEN  OF  THE  B&K  2032. 

120  ' 

130  '  THE  PROGRAM  WILL  STORED  THE  SPECTRUM  FROM  23  TO  27  HZ  It 

THE 

140  '  B&K  2032  IS  DISPLAYING  THE  0  TO  50  HZ  SPECTRUM. 

150  ' 

160  '  THE  PROGRAM  USES  THE  SAME  SUBROUTINES  AS  "BKDISP" 

170  ' 

180  ************************************************************ 
1000  OPTION  BASE  1  'ARRAYS  BEGIN  WITH  ELEMENT  1 

1010  GOSUB  9000  'INITIALIZE  THE  INTERFACE 

1020  GOSUB  9200  'LOCAL 

1500  CLS 

5500  INPUT  "ENTER  DATA  FILE  NAME" ; DF$ 

5505  GOSUB  9100  'REMOTE 

5510  OPEN  "0" , #1, DF$ 

5530  FOR  I  =  368  TO  432  'SAMPLES  THE  CENTER  64  POINTS 

5540  0$  =  "AF  "+STR$ ( I ) 

5550  GOSUB  9300 

5560  GOSUB  9400 

5570  DAT  =  VAL (MID$(E$,2,6) )  'DATA  CONVERSION 

5580  WRITE  # 1 , DAT 

5590  NEXT  I 

5600  CLOSE  #1 

5605  GOSUB  9200 

5610  GOTO  1020 

7000  GOSUB  9200 

7010  END 

9000  DEF  SEG  =  &HD000 

9010  CINIT$  -  "SYSCON  MAD=3 ,  CIC=1,  NOB=l,  BA0=&H300" 

9020  CREMOTE?  =  "REMOTE  26" 

9030  CLCCAL$  =  "LOCAL  26" 

9040  COUTPUT$  =  "OUTPUT  26 [ $  1  " 

9050  CENTER$  =  "ENTER  26 [E 

9060  IE488  =  0 :  A%  =  0:  FLG%  =  0:  BRD%  =  0 
9070  CALL  IE488  (CINIT$,  A% ,  FLG% ,  BRD% ) 

9080  IF  FLGsoO  THEN  PRINT  "INSTALLATION  ERROR:  ";HEX$(FI.G  )  :  END 
9090  RETURN 

9100  CALL  IE488  ( CREMOTE$ ,  A% ,  FLG% ,  BRD°s ) 

9110  IF  FLG%<'  >0  THEN  PixINT  "REMOTE  ERROR:  "  ;  HEX$  ( 1T.G  :  )  :  END 
9120  RETURN 
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9200  CALL  IE488  (CLOCAL$,  A% ,  FLG% ,  BRD% ) 

9210  IF  FLG%<>0  THEN  PRINT  "LOCAL  ERROR:  ";HEX$(FLG%) 
9220  RETURN 

9300  OUTPUTS  =  0$+CHR$ ( 10 ) 

9310  CALL  IE488  (COUTPUT$,  OUTPUTS,  FLG% ,  BRD% ) 

9320  IF  FLG%<>0  THEN  PRINT  "OUTPUT  ERROR:  ";HEX$(FLG% 

9330  RETURN 

9400  E$  =  SPACES ( 15 ) 

9410  CALL  IE488  (CENTERS,  E$,  FLG% ,  BRD% ) 

9420  IF  FLG%<>0  THEN  PRINT  "ENTRY  ERROR:  ";HEX$(FLG%) 
9430  RETURN 
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:  END 


)  :  END 


:  END 
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AC  measurements  were  successfully  accomplished  for  a  2 5  Hr 
magnetic  field  by  inspecting  the  modulation  of  the 
interferometer's  output  light  at  25  Hz.  Sensors  composed  of  >■ 
sputtered  metal  coated  fibers  did  not  detect  the  AC  magnetic 
field  because  the  thin  metal  films  did  not  generate  c agr  f  •  : 
to  change  the  length  of  the  optical  fiber  a  detectable  amount. 
Inc  sensor  did  detect  the  AC  magnetic  field  wh-_ n  the  sen.  :  : 
was  composed  or  a  bulk  nickel  rod.  Thermal  annealing  the  n  a:*: 
rod  in  the  presence  of  an  external  magnetic  field  improve  5  t  *> 
sensor's  response.  Thinning  the  optical  fiber  in  the  sc-rrr  r  .= 
did  not  significantly  affect  the  sensor's  perform  an  re  ,  :  .:*■ 
increasing  the  length  of  the  optical  fiber  did  enhance  the 
sensor's  response.  The  sensor's  performance  also  improved  wh.e 
the  AC  magnetic  field  was  superimposed  on  a  DC  bias  magnet  . 
field.  In  the  sensor's  optimum  configuration,  an  AC  magnet u: 
field  flux  density  of  10  mG  was  detected. 


